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Physies—the Profession 


BY HENRY A. BARTON 


Director, American Institute of Physics 


HYSICS, the profession, has come of age, or 

at least has now clearly become economically 
and immediately useful and profitable. ‘Physi- 
cist’’ seems at last to be an enviable title! Now 
those who have earned the right to own it are 
beginning to realize that they have-to protect it. 
The headaches as well as the pride of ownership 
are at hand. We can hardly hope to avoid the 
consequences of worldly recognition and achieve- 
ment with which we have watched our friends 
the chemists and engineers struggle. 

It is time to shape a plan and policy designed 
not only to provide an insistent safeguard for the 
quality of the title ‘‘physicist,’’ but specifically 
further to promote the vocational recruiting of 
able young men, to require high standards in 
their education, to fight for appropriate rewards 
for their work, to eliminate imposters and to 
guide social forces such as labor unions and 
government regulations to a proper understand- 
ing of the profession’s character and needs. 

This is a program 
policy 


the outlining of a plan and 
which cannot be carried out in a day, 
nor without trial and error and the participation 
of many thoughtful men. The Executive Com- 
mittee of the American Institute of Physics so 
understands the problem and is considering how 
to organize a study and assemble opinions to 
obtain—even now when physicists are preoccu- 
pied with war—a guiding draft which may be 
considered point-by-point for adoption. Probably 
a committee will be formed to report a plan. 
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Certainly all physicists should concern them- 
selves with the subject. 

Our problems being different in several ways 
from those of the chemical profession, the engi- 
neering profession and others, we cannot expect 
to find in their experience a readymade procedure 
to answer our needs. It is we who best know our 
profession and surely we have the mental re- 
sources to conceive how its interests are best to 
be served. Nonetheless we can learn quickly and 
find abundant food for thought from studying 
the conclusions of others. For example, the 
American Chemical Society has recently care- 
fully and exhaustively attacked several aspects 
of the “profession’’ problem. The _ resulting 
recommendations concerning employer-employee 
relations constitute a document of the greatest 
interest. It is not and could not be short, but it 
is worth quoting below as an assignment of 
almost ‘‘required reading’’ for physicists. 

In the perusal of this document, certain impli- 
cations must be borne in mind. For example, a 
headquarters organization is implied which has 
resources to serve members with legally as well 
as professionally sound advice, perhaps political 
influence and in other capacities to a degree not 
maintainable with the funds and _ staff now 
available in the American Institute of Physics. 
This organization would have to make a sys- 
tematic continuing study of the excellence of 
institutions for the advanced training of physi- 
cists—a study directed by an authoritative 





_—— 


board empowered to recommend whether or not 
any institution should be “accredited” for such 
training. It would also have to maintain a board 
of altruistic leaders of the profession which, like 
the Committee on Economic Status of the 


chemists, is willing “to act in an advisory 
capacity to both employers and employee groups 

. . and, if deemed helpful, to publicize its find- 
ings.”’ Active and effective Boards and Com- 
mittees, such as these would have to be, require 
much staff work in the headquarters office. 

The recommendations of the American Chemi- 
cal Society which follow are quoted from the 
December 25, 1941, issue of the News Edition of 
that Society: 


Employer-Employee Relationships for Professional Chem- 
ists as Recommended by the American Chemical Society, 
Unanimously Adopted by the Board of Directors, Sep- 
tember 8, 1941. 


Because of efforts to compel chemists and chemical 
engineers to join labor organizations in order to obtain or 
retain employment in certain plants, the Board of Directors 
of the American Chemical Society has given consideration 
to the broad problems of employment in the field of 
chemistry. 

So that the position of the Society may not be misunder- 
stood, the Board of Directors issues the following statement 
for the more complete information of our membership: 

The Society has taken no stand against ‘‘collective 
bargaining” for professional men when such bargaining is 
not controlled by nonprofessional groups and where the 
bargaining unit is composed exclusively of professional men. 

The Society condemns no one of its members for joining 
any noncoercive labor union so long as he does so 
voluntarily. 

The Society, however, is unalterably opposed to the 
forcible inclusion of professional men in bargaining units 
dominated and controlled by nonprofessional employees, 
whether that inclusion be brought about by economic pres- 
sure upon an employer, by intimidation of the professional 
employee, or by operation of either state or federal law. 

The Society will bend every effort to maintain for all its 
members the ‘right to work” and the “right to employment 
and promotion”’ on the basis of worth and merit. 

Accordingly, the Board of Directors goes on record as 
opposed to affiliation of its members with any organization 


‘that conditions promotion primarily on the basis of 


“seniority,” or that insists that they join any labor 
organization where they would be in a minority, with no 
power to protect themselves while paying ‘‘protection”’ 
thereto as an essential to the privilege of earning a liveli- 
hood and with their wage scales negotiated by those whose 
selfish interest would require that benefits be sought for the 
larger number of nonprofessional workers to the detriment 
of the relatively few professional employees. 
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Letters from some members, supported by facts and 
intelligent argument, claiming that conditions during the 
depression and under the emergency have secured for 
tradesmen and laborers, under union leadership, advances 
in income not enjoyed by trained professional men often 
directly associated with the former, have been given 
consideration. 

Already a number of our more progressive employers 
have made surveys of their personnel and have taken, or 
propose to take, corrective measures. 

In order that some logical conclusion may be reached, it 
seems wise to consider the distinction between professional 
and nonprofessional employees engaged in chemical work 
in order properly to differentiate between professional and 
nonprofessional workers. 

At the present time there are a large number of tech- 
nicians employed in the field of chemistry. 

The term “technician” may be applied to routine workers, 
trained in the laboratory but with no special education or 
mental proficiency in chemistry and allied sciences. Nor- 
mally they will be individuals with secondary school 
education only, or perhaps ones who have failed to acquire 
a baccalaureate degree through lack of funds, insufficient 
effort or capacity, or attendance at underequipped educa- 
tional institutions. These are almost always paid wages on 
the hourly basis and are often included in labor union 
agreements. No proper objection can be raised to their 
inclusion in a bargaining unit composed of employees doing 
various kinds of skilled, semiskilled, or unskilled work. 

The two groups of employees which should not be sub- 
jected to forcible inclusion in a heterogeneous bargaining 
unit are those who have received degrees in chemistry or 
chemical engineering and are engaged in those fields, either 
as chemical internes or as professionals. 

The term 


chemical interne’? may be applied to those 
who are essentially in the final qualifying stage for their life 
work in the field of chemistry. They have received their 
baccalaureate degrees with majors in chemistry or chemical 
engineering, have proved their proficiency not only in 
chemistry but also in mathematics, physics, modern lan- 
guages, etc. They are acquiring thereby the necessary 
training and experience to qualify for full professional 
status or standing. Those in industrial corporations are 
engaged in professional work on problems confidential to 
management, whether their work be control, research, or 
development. 

The term ‘‘professional’’ should be applied only to those 
who have the baccalaureate degree, or its equivalent by 
specific accomplishment, and who, having been graduated 
from institutions approved by the American Chemical 
Society, have had at least two years of postgraduate training 
in chemistry or chemical engineering in institutions of like 
grade or have for an equal period obtained experience in 
chemical work. For graduates of other educational insti- 
tutions five years of postgraduate training and/or experi- 
ence subsequent to the baccalaureate degree should be 
required. These are the minimum requirements for full 
professional membership in the American Chemical Society. 

The Board of Directors has heretofore hesitated to make 
any suggestions on the matter of adequate compensation. 
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On the whole our educational institutions and chemical 
corporations have realized the value of brains and have 
compensated employees accordingly. The chemical pro- 
fession as a whole is on a level with any other profession in 
material return. There are today in our profession, as in 
many others, individuals whose compensation is high, 
indeed almost fabulous, and others who are grossly under- 
paid. There are many who, even in the higher brackets, 
earn more than their income. There are some, even in the 
lower brackets, whe are not worth their hire, and who, for 
their own sakes and for the reputation of the profession, 
should seek other more fitting fields for a livelihood. In 
spite of this generalization there has come to the attention 
of the Board of Directors evidence that there are today 
many cases where worthy professional men with years of 
study and training are grossly underpaid and are receiving 
less compensation than men without any specialized educa- 
tion employed in many of the trades and crafts. 

In the opinion of the Board of Directors, the time has 
arrived when it must publicly suggest to all employers of 
chemists and chemical engineers that they review with care 
their employer-employee relations in this line of endeavor. 
The American Chemical Society stands ready through a 
carefully chosen Committee on Economic Status to act in 
an advisory capacity to both employers and employee 
groups of its members and, if deemed helpful, to publicize 
its findings. Enlightened management will, it is believed, 
welcome consultation with such a committee composed of 
carefully chosen individuals. 

As a beginning, the Board of. Directors suggests that all 
worthy chemical internes and professionals when engaged 
in professional work be paid on a salary and not on an 
hourly basis. 


The extent of training, experience, and capacity, as well 
as individual personality and merit, are factors that must 
be considered when employing one to engage in professional 
work. Bearing these factors in mind, as well as the objects of 
this Society, and the requirements for membership therein, 
the board ventures to make the following suggestions to 
employers for their consideration in respect of starting 
salaries of chemical internes and professionals under present 
economic conditions, when such persons are employed to do 
professional work: 

A minimum starting salary not less than $1,500 per 
annum for any chemical interne; and one of not less than 
$1,800 per annum for any chemical interne of better than 
average ability or training, including especially those 
graduated from institutions whose chemical work is of 
particularly high quality such as those approved by the 
American Chemical Society; and a minimum starting 
salary of not less than $2,400 per annum for any individual 
who has attained ‘professional’ grade as hereinbefore 
described, and, according to the extent of training and 
capacity, higher starting salary in proportion to training, 
experience, merit, and individual accomplishment. 

The foregoing suggestions are for minimum starting 
salaries only and are not to be regarded as suggestions for 
maximum salaries or as recommendations for a wage scale. 
They are made in respect of only those who are engaged for 
professional work in the field of chemistry and chemical 
engineering, as the Society cannot concern itself with 
problems of employment in other fields. 

Inasmuch as technicians normally come under nonpro- 
fessional status, the technician grade does not come within 
the purview of this organization. 





Calendar of 


February 

13-14 American Philosophical Society, Philadelphia, Pennsylvania 

16-19 Technical Association of Pulp and Paper Industry, New York, 
New York 

20-21. American Physical Society, Detroit. Michigan 

March 

12-13 Society of Automotive Engineers, New York, New York 

23-25 American Society of Mechanical Engineers, Houston, Texas 

April 

10-11 American Physical Society, Southeastern Section, Oxford, 
Mississippi 

16-17 National Petroleum Association, Cleveland, Ohio 

19 American Ceramic Society, Cincinnati, Ohio 
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Meetings 

20-24 American Chemical Society, Memphis, Tennessee 

23-25 American Philosophical Society, Philadelphia, Pennsylvania 
30-May 1 American Physical Society, Baltimore, Maryland 

May 

11-13 American Institute of Chemical Engineers, Boston, Massa- 


chusetts 
31-June 5 Society of Automotive Engineers, White Sulphur Springs, 


West Virginia 


June 

17-20 American Physical Society, Pacific Coast Meeting, Salt Lake 
City, Utah 

25-27. American Physical Society, State College, Pennsylvania 
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Conference on the Hardening of Metals 


BY FREDERICK SEITZ 


Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 


HE development of the theory of crystalline 

matter during the past ten years has paved 
the way for a deeper understanding of those 
properties of metals that are of most interest to 
engineers and metallurgists. Among these prop- 
erties there are probably none of greater concern 
than those centered about the subject of hard- 
ness. For this reason it is both timely and 
expedient that the physicists whose interest lies 
in the field of solids devote a share of their 
effort to extending the understanding of these 
properties. With a view toward furthering this 
course, the Department of Physics of the Uni- 
versity of Pennsylvania and the Philadelphia 
Chapter of the American Society for Metals 
jointly sponsored a conference last May on 
“Hardening of Metals’ and brought together 
both the physicists and metallurgists whose 
primary field of research deals with the basic 
properties of metals. The following papers were 
presented : 


The General Concept of Metals, Dr. W. 
Bell Telephone Laboratories 

The Basic Problems of the Deformation of Metals, Dr. 
S. L. Hoyt, Battelle Memorial Institute 

The Hardening of Iron, Proressor R. F. MenL, Carnegie 
Institute of Technology 

Precipitation Hardening, Dr. W. L. 
Company of America 

Hardening by Alloying Agents, Dr. S. SIEGEL, Westing- 
house Research Laboratories 

Hardening by Cold Work, Proressor G. Sacus, Case 
School of Applied Science 

X-Rays and Hardening, PRoressor C. S, 
Carnegie Institute of Technology 

The Strength of Metals in High Speed Deformation, 
Dr. A. NapAl, Westinghouse Research Laboratories. 


SHOCKLEY, 


FINK, Aluminum 


BARRETT, 


It was hoped initially that printed versions of 
all eight papers would be made generally avail- 
able in the pages of the Journal of Applied 
Physics. The pressure of other occupations 
proved too great for this, however, and only two 
of the contributors were able to furnish manu- 
scripts. Fortunately, several of the speakers have 
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written summarizing articles on their subjects in 
this journal or in others in the past few years. 
In particular, the interested reader is referred to 
the following: 


S. L. Hoyt, “Plastic deformation in metals,” Metals 
Progress 38, 659 (1940). 

W. Shockley, ‘‘Nature of the metallic state,” 
Phys. 10, 543 (1939). 

A. Nadai and M. J. Manjoine, “High speed tension 

Parts I, IT, III, 


J. App. 


tests at elevated temperatures,” 
Proc. A. S. T. M. (1940-41). 


It is probably not out of place to summarize 
briefly the principal lines along which the theory 
of hardness is developing. It seems to be generally 
agreed at the present time that all pure metals 
would be stronger and harder by two or three 
powers of ten if it were not for the occurrence 
of weak spots. These weak spots seem to aid 
rupture directly by allowing cracks to start for 
lower stresses than would otherwise be the case, 
and seem to aid plastic flow by assisting in the 
formations of a type of lattice defect which can 
be propagated through the crystal under applica- 
tion of very low stresses. The precise nature of 
this type of lattice defect is still a matter of 
speculation; however, a particular model known 
as a “‘dislocation”’ has been employed extensively 
in theoretical work during the past few years 
and appears to possess most of the features that 
are necessary. A summary of this work may be 
found in the February, March, June and July 
issues of this Journal for the past year. This 
tvpe of lattice defect was first proposed by 
Prandtl to explain internal friction and was 
subsequently introduced into the theory of 
plastic flow by Taylor and Orowan. From the 
standpoint of this theory metals are more ductile 
than most other solids of comparable cohesive 
properties because the energy required to form 
dislocations in them and the forces required to 
make the dislocations move are lower than for 
the other solids. From the same standpoint, 
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much of the technology of hardening metals is 
centered about procedures which have the effect 
of inhibiting the flow of dislocations. It is be- 
lieved at the present time that most of the 
hardening practices produce residual stresses on 
an atomic or submicroscopic scale and that the 
dislocations experience difficulty in passing these 
regions because of stresses opposing their motion. 
Thus cold-work and alloying agents harden 
metals because they distort the lattice and 
greatly increase the stress required to make dis- 
locations move. It has recently been shown by 
Mott and Nabarro that alloying agents are most 
effective in producing residual stresses when 


present not in the most highly dissolved form, 
as are the atoms of zinc in brass, but when 
clumped together in the form of finely divided 
precipitate. It is very likely that the age-harden- 
ing properties of Duralumin and the quench- 
hardening properties of steel are related to 
residual stresses produced by such precipitates. 

The type of research that is most urgently 
needed in order to complete the unraveling of 
the theory of hardness is that which the physicist 
can do more readily than any other type of 
scientist. There is probably no field of physics in 
which research will find more use in as short a 
period of time. 





Precipitation Hardening * 


By WM. L. FINK 


Aluminum Company of America, New Kensington, Pennsylvania 


BOUT 1910, Alfred Wilm,! a German metal- 
lurgist, made an epoch-making discovery 

a method of heat treating certain aluminum 
alloys, especially an aluminum-copper alloy 
known as Duralumin. This treatment consisted 
of heating the alloy to approximately 500°C (now 
known as a solution heat treatment), quenching, 
and then allowing the material to stand at room 
temperature for a few days (now known as 
precipitation heat treatment or aging). The 
precipitation treatment causes a_ substantial 
increase in the strength and hardness of the 
alloys. Figure 1 shows the effect of room tempera- 
ture aging on the tensile strength of three 
aluminum-copper alloys. It will be noted that 
the tensile strength increases rapidly during the 
first few hours and then progressively slower, so 
that after about four days it has reached a nearly 
constant value. The upper graph of Fig. 2 shows 
the effect on various properties of a binary 
aluminum-copper alloy which was quenched in 
* Presented at Conference on Hardening of Metals 
held at the University of Pennsylvania on May 17, 1941. 
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cold water and given a precipitation treatment 
at 100°C instead of room temperature: When 
such a temperature is employed for the precipi- 
tation heat treatment, the process is often re- 
ferred to as “artificial aging.’’ Except for some 
inflections at short aging times which will be 
discussed later, the tensile strength, yield 
strength, and hardness increase and the elonga- 
tion decreases in a manner consistent with the 
tensile strength changes shown by Fig. 1. 

At the time Wilm made his discovery, he 
thought the heat treatment of the aluminum 
alloys was substantially the same as the heat 
treatment of steel. He had found, by thermal 
analysis, an arrest point which he called the 
recalescence point because he thought it was 
similar to the recalescence point in steel. It was 
later found that the thermal arrest was really the 
melting point of a complex eutectic and that 
aluminum, in fact, has no allotropic transfor- 
mation, but for several years no one was able to 
furnish a satisfactory explanation of this new 
method of hardening. 
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Then in 1919, Merica, Waltenberg, and Scott? 
of the Bureau of Standards published a paper, 
“Heat treatment of Duralumin,” 


which finally 
gave the explanation. The following two para- 
graphs from the paper summarize the theory: 
“Upon slowly cooling an alloy containing 3 
percent of copper from 500°C, the CuAl. precipi- 
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Fic. 1. Effect of aging three aluminum-copper alloys at 
room temperature. 


tates from solid solution to maintain equilibrium 
along the line b-e [See Fig. 3—(Fig. 12 of 
Merica’s paper) |. The alloy so obtained is soft 
and does not harden upon aging. Rapid cooling 
of the same alloy from 500°C, by quenching, 
partially or wholly suppresses this precipitation 
of CuAl.. If the alloy is held at a low tempera- 
ture, such as that of liquid air (—180°C), no 
further hardening takes place upon aging. The 
alloy is not in equilibrium, but the rate of nuclear 
formation of CuAls is so small that no CuAl. 
precipitates to bring about equilibrium. If, how- 
ever, the temperature of the alloy is raised to 
100°C or even to ordinary room temperature, 
according to the theory which the authors 
propose, the mobility becomes sufficiently great 
- that precipitation of the CuAl, takes place in the 
form of very fine particles of colloidal dispersion. 
To this precipitation is due the hardening during 
the aging of Duralumin.” 

The question of the optimum size of the 
colloidal particles is further considered in the 
second paragraph. “Upon aging a quenched 
specimen, at 200°C, for example, the hardness 
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first increases to a maximum and afterward 
decreases. During that aging there has been first 
a formation of fine nuclei of CuAl., followed by 
coalescence of these particles into ones of larger 
size. There is, therefore, a certain average size of 
particle of CuAl, for which the hardness of the 
material is the maximum; atomic dispersion of 
the solute is not the dispersion that produces the 
maximum hardness, but some intermediate one 
between it and that at which the particles become 
visible by ordinary means.” 

This concise and accurate description of the 
hardening process stimulated further work. It 
soon became apparent that the heat treatment of 
Duralumin was not the unique phenomenon 
which had been believed. Metallurgists found 
that it was possible to develop precipitation- 
hardenable alloys of many metals with special 
combinations of properties for certain 
mercial uses. 


com- 


However, even with this guiding principle 
much experimentation is required to develop new 
alloys. If enough were known about the details of 
precipitation hardening, it might be possible to 
reduce the experimental work. Therefore, let us 
consider this precipitation hardening process in 
some detail in the light of our present-day 
knowledge. 

Precipitation hardening is one example of the 
general problem of the decomposition of one 
phase into two phases. Much experimental and 
theoretical work has been done on this subject. 
Tammann’ and his students have studied the 
crystallization of organic compounds and glasses. 
Zsigmondy* and his students have used the 
ultramicroscope in studying the formation of 
colloidal precipitate in transparent solutions. 
Gibbs,® Einstein,® Vollmer and Weber,’ Stronsky,$ 
Becker,® and Déring"® have all contributed to the 
theory. In this paper no attempt will be made to 
review their work but merely to apply the 
general ideas to the age hardening of alloys. The 
aluminum-copper alloy will be used as an ex- 
ample, but it should, of course, be understood 
that the ideas apply generally and were, in fact, 
developed on other much of the 
pioneering work being done on the precipitation 
of droplets of liquid from a gas. 

As stated previously, the first step in the heat 
treatment of an aluminum-copper alloy consists 


systems, 
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in heating it above the solvus curve in order to 
dissolve the copper constituent ; i.e., atomically to 
disperse the copper in the aluminum. The distri- 
bution of solute atoms at this stage is of im- 
portance and has been expressed by Einstein® in 
the form of a differential equation, but the 
equation cannot be used readily in making 
calculations. Fortunately, at elevated tempera- 
tures, for example just below the solidus, a 
random distribution is probably approached so 
that an approximation can be made by using 
Poisson equation.'! Table I, calculated on this 
basis, shows the number of groups of twelve 
atoms in one cubic millimeter which contained a 
given number of copper atoms. It will be noted 
that the last group in the table (i.e., four copper 
atoms in a group of twelve) has the composition 
of the precipitate, CuAl,. There are a large 
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Fic. 2. Precipitation hardening of an aluminum-copper 

alloy (4.00 percent Cu) at 100°C. Upper graph is for 


material quenched in cold water, and lower graph for 
material quenched in oil at 100°C. 
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number of them—650 trillion per cubic milli- 
meter. It should be realized that at the solution 
heat treating temperature these groups are not 
stable, but are constantly forming and disinte- 
grating by the chance diffusion of copper atoms. 
However, a statistical equilibrium is obtained, 


TABLE I. Probable numbers of groups, containing a 
total of 12 atoms with various numbers of copper atoms, 
in one cubic millimeter of an Al-Cu alloy (4.78% Cu). 


No. of No. of groupst 
Cu Atoms P(s)* in 1 mm* 
0 0.7788 39 X 10"? 
1 0.1947 97 X 1016 
2 0.0243 12x 10% 
3 0.0020 10 10% 
4 


0.00013 65 x 10" 


* Probability of occurrence of J Cu atoms in a group of 12 atoms. 
+ Total number of groups, containing 12 atoms in 1 mm* =50.8 X10!7. 


and the number of groups occurring in a volume 
as large as one cubic millimeter would be 
substantially constant. 

Upon quenching, the conditions which existed 
in the solid solution at the heat treating tempera- 
ture are maintained momentarily. However, the 
solid solution is no longer stable but tends to 
decompose into two phases—a more dilute solid 
solution, and that crystal form of CuAl, known 
as 6’. Assuming that this decomposition occurs in 
the normal manner (i.e., according to the general 
method which has been found by the investi- 
gators mentioned above), the 6’ phase will 
develop from nuclei formed by the chance 
association of solute atoms. 

In crystalline materials (both metallic and non- 
metallic) the precipitate particles are oriented in 
a definite manner with respect to the parent solid 
solution. This so-called Widmanstiatten relation- 
ship has been studied in considerable detail by 
Dr. Mehl and his collaborators.” The basic 
requirement is an orientation such that a good 
matching of atoms exists at the interface between 
the two phases. The aluminum-copper system is 
no exception. Preston has established the 
Widmanstiatten relationship shown in Fig. 4. 
The figure at the left shows the aluminum-copper 
solid solution with a chance grouping of four 
copper atoms (the full circles) and a region 
blocked out around them corresponding to one 
unit cell of 6’. The figure on the right shows one 
unit cell of 6’ with the dimensions which were 
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observed in the case of large precipitate particles. 
It will be noted that basal plane of the 6’ lattice 
has exactly the same arrangement of atoms as 
the cube plane of the aluminum-copper solid 
solution. In accordance with the finding of Dr. 
Mehl and his collaborators, one would expect a 
nucleus to grow preferentially in the plane in 
which the matching is good. That is, in this case, 
the precipitate 6’ would be expected to grow as 
plates parallel to the cube planes of the solid 
solution. 
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Fic. 3. Aluminum end of aluminum-copper equilibrium 
diagram according to early work of Merica, Waltenberg 
and Scott. 


A few years ago, Guinier™ and his co-workers, 
and Preston™ independently, found streaks on 
Laue patterns of single crystals of aged aluminum- 
copper alloys, of the type which might be ex- 
pected from very thin plates of @’ parallel to the 
cube plane of the aluminum lattice. According to 
these investigators, the platelets formed in a 
short time at room temperature are only about 
one unit cell thick, and less than 10 unit cells in 
diameter, but the size increases with increasing 
temperatures and time. Guinier and Preston did 
not commit themselves in regard to the structure 
of these platelets prior to the time when they 
became so thick that they acted as three-dimen- 
sional x-ray gratings and were definitely identified 
as 6’. However, these streaks are consistent with 
the theory that precipitation in the aluminum- 
copper alloy occurs by the normal process of the 
formation and growth of nuclei of 6’ in ac- 
cordance with the Widmanstiatten mechanism. 

Barrett, Geisler, and Mehl'® have found the 
same types of streaks in the case of an aluminum- 
silver alloy. It is very significant that in this 
alloy the precipitate (y’) forms plates on the 


78 





octahedral planes of the aluminum, in accordance 
with the Widmanstatten relation. The formation 
of precipitate particles on different planes of the 
aluminum lattice, depending upon the structure of 
the precipitate, is strong evidence that the pre- 
cipitate platelets grew from nuclei which had the 
structure of the precipitate (i.e., 6’ or y’). The 
weight of the evidence will further increase if 
additional systems all show the same behavior. 

Another important factor in precipitation 
hardening is the nucleus size. Both theoretical 
and experimental work have indicated that the 
higher the temperature at which aging occurs, the 
larger must be the nucleus for stability and 
growth. Although no quantitative determination 
has been made, it has been suggested'*® that in the 
aluminum-copper alloys a nucleus of twelve 
atoms (i.c., one unit cell of 6’) is stable at room 
temperature and below, but that by the time the 
aging temperature reaches 100°C, the smallest 
stable nucleus probably contains about 48 atoms 
or about 4 unit cells. More will be said about 
this in connection with electrical conductivity 
changes. 

In view of some of the theories to be discussed 
later in this paper, it is well to point out that the 
growth of nuclei into larger precipitate particles 
involves only the normal type of diffusion; that 
is, the attachment of solute atoms to the growing 
crystal diminishes the concentration of solid 
solution immediately adjacent to the precipitate. 


= NGA - t= 
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Fic. 4. Widmanstitten relationship between 6’ precipi- 
tate and aluminum-copper solid solution, as well as 
crystal structure of 6’. 


Other solute atoms then diffuse from points 
of high concentration to this region of low 
concentration. 

Such is the mechanism which would be ex- 
pected to apply to the age hardening of aluminum 
alloys in the light of our general knowledge of the 
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Fic. 5. Changes in electrical resistivity of a quenched aluminum-copper alloy (4 percent Cu) 
at room temperature. 


decomposition of solid solutions. The question 
then naturally arises, ‘“Are the property changes 
actually observed during precipitation hardening 
those which would be expected?” It was pointed 
out by Jeffries and Archer" that small particles 
of precipitate might act as keys to obstruct slip 
and thus strengthen the metal. This mechanical 
conception is probably an oversimplification, and 
possibly an explanation in terms of “dislocations” 
and interatomic forces would be better. However, 
this ‘‘slip-interference theory” did help in visu- 
alizing the process and has had a big influence in 
bringing about a general agreement that precipi- 
tate particles of the right size will produce higher 
vield strengths, tensile strengths, and hardness 
than the same amount of alloying element in 
solid solution. Regardless of the theory, experi- 
mental proof of precipitation hardening lies in 
the fact that in certain alloys some of the 
precipitate particles are large enough to be 
resolved microscopically at the stage of maximum 
hardness or earlier. 

In regard to volume changes which occur 
during the early stages of age hardening, a 
semiquantitative calculation can be made if 
the crystal structures of the phases and the 
Widmanstatten relation are known. Enough data 
are available to make such calculations in the 
case of the aluminum-copper system. It was 
found!’ that the calculated volume change was 
in the right direction and of the proper order of 
magnitude. As the aging progresses, particles 
become larger and the calculations become very 
complex. Changes in the concentration of the 
solid solution by diffusion; complete disregistry 
of the lattices of the precipitate and the matrix; 
the complex internal strains are all important 
factors which are difficult to evaluate. However, 
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no facts are known which are not in qualitative 
agreement with the mechanism just outlined. 

In regard to electrical conductivity changes, 
accompanying precipitation hardening, wave 
mechanics indicate that particles substantially 
larger than single solute atoms would offer 
maximum resistance to the flow of electrons. 
Professor Mott!’ in England has estimated that 
the probable diameter for maximum resistivity is 
four or five interatomic distances. It is of con- 
siderable interest that a physicist in this country 
who is an expert in the field of wave mechanics 
but was unacquainted with the phenomenon of 
age hardening predicted qualitatively the changes 
in electrical conductivity which are actually 
observed. When an aluminum-copper alloy is 
aged at room temperature, the stable nuclei are 
smaller than the particles which offer maximum 
resistance to the flow of current. Therefore, the 
resistivity first increases to a maximum and then 
decreases as shown in Fig. 5. 

If an aluminum-copper alloy were aged first at 
room temperature and then at an elevated 
temperature (say 150°C or 200°C), many of the 
small precipitate particles would be expected to 
dissolve rapidly immediately after heating to the 
higher temperature, because many of the particles 
would be smaller than the nuclei which are stable 
at the higher temperature, and because the solid 
solution immediately adjacent to the particles 
would be undersaturated at the new aging 
temperature. This re-solution of the small parti- 
cles would be expected to soften the material and 
decrease the resistivity. These changes do, in 
fact, occur as illustrated in Figs. 6 and 7. 

In brief, there are no experimental facts 
pertaining to precipitation hardening alloys 
which are inconsistent with the assumption that 
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precipitation hardening like the decomposition of 
any other crystalline solid solution occurs by the 
formation and growth of nuclei of a second 
phage in accordance with the Widmanstitten 
mechanism. 


This would be a good place to end this paper if 
the development of the theory of age hardening 
had been straightforward and universally ac- 
cepted. Such was not the case, and there are 
metallurgists even today who still adhere to 
widely different theories. Consequently the 
history will be considered very briefly. 

Prior to publication of the paper by Merica, 
Waltenberg, and Scott,? metallurgists had gener- 
ally agreed that a given alloy possessed maximum 
hardness and electrical resistivity when in the 
solid solution state. One can search the literature 
in vain for any data on which to justify such an 
assumption. Perhaps the most plausible reason 
for coming to such a conclusion is that in the 
microscopic range smaller constituent particles 
produce greater hardness and resistivity than the 
same weight of large particles. One might there- 
fore jump to the conclusion that the finest 
possible particles (the atoms) would produce 
maximum hardness and resistivity. Obviously 
this is an unwarranted extrapolation and one 
which a student of colloidal chemistry would not 
be apt to make. Unfortunately metallurgists had 
generally accepted this hypothesis. 

When Merica, Waltenberg, and Scott proposed 
their theory of age hardening, there were many 
metallurgists who thought it could not be true 
because everyone knew that an alloy possessed 
maximum hardness when in the solid solution 
state. In fact, several years after the publication, 
a famous metallurgist still contended that the 
solid solution state was the hardest and that 
precipitation hardening probably did not exist, or 
if it did, it was indeed anomalous. Finally, 
however, the facts were so conclusively proved 
that all metallurgists had to agree that precipi- 
tate particles of the right size offer more resist- 
ance to deformation than does random atomic 
dispersion. 

One would naturally suppose that rejection of 
the hypothesis that the solid solution state 
possessed maximum hardness would have led to 
the rejection of the twin hypothesis that the 
solid solution state possessed maximum electrical 
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resistivity, inasmuch as both hypotheses rested 
on exactly the same untenable line of reasoning. 
This, however, was not the case. When it was 
found that room temperature aging of Duralumin 
caused an increase in resistivity, it was immedi- 
ately assumed by some metallurgists that there 
could be no precipitation during room tempera- 
ture aging. It was argued that, since the solid 
solution state possessed maximum resistivity, an 
increase in resistivity could not be consistent 
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Fic. 6. Schematic representation of hardness changes 
in a quenched aluminum-copper alloy which is aged first 
at a low temperature (e.g., room temperature) and subse- 
quently at elevated temperatures. 


with precipitation but must indicate some other 
type of solid solution with still greater resistivity. 
The metallurgists who reasoned thus divided age- 
hardening alloys into two groups. The first group 
in which the resistivity decreased during aging 
were “normal alloys’’ and the hardening was 
admittedly caused by the precipitation of fine 
particles. The second group in which the re- 
sistivity increased were ‘‘anomolous,”’ and the 
hardening of these alloys was presumed to be 
caused by some sort of a pre-precipitation solid- 
solution reaction. 

Once this view was expressed, others came 
forward with various arguments to support it, 
none of them carrying very much weight by 
themselves, but all taken together were effective 
in guiding popular opinion away from the theory 
of precipitation hardening. For example, it was 
frequently repeated that no precipitation could 
be observed under the microscope during the 
early stages. This certainly was a feeble argu- 
ment since it would be expected that particles of 
colloidal size (i.e., submicroscopic particles) 
would harden the alloys most effectively. 
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Finally, however, it was found that no change 
in lattice parameter could be measured by the 
sensitive x-ray back reflection method during 
room temperature aging of aluminum-copper 
alloys. This was generally considered to be a fatal 
blow to the simple theory of precipitation 
hardening. It was known that parameters could 
be measured with considerable precision. It was 
also known that the parameter of solid solution 
varied with concentration. It was argued, there- 
fore, that with no measurable change in parame- 
ter there could be no significant precipitation. 

Consequently concoction of pre-precipitation 
theories was the order of the day. It is impossible 
to review all of these theories in this paper. 
However, here is one theory for the age hardening 
of aluminum-copper alloys proposed by a promi- 
nent metallurgist as late as 1937: It was tacitly 
assumed that solute atoms in these alloys are 
distributed uniformly instead of having the 
random distribution which characterizes most 
solutions. This necessitated the 
mechanism for grouping copper 


invention of a 
atoms. It was 
assumed that the first stage of aging was the 
diffusion of copper atoms to certain planes of 
the solid-solution lattice to form a kind of 
imperfect ordered arrangement which, by defini- 
tion, had a high electrical resistivity (probably 
resulting, it was believed, from work hardening 
caused by the diffusion). The next step in the 
aging was then supposed to be the formation of 
CuAl, molecules. Just why copper atoms sur- 
rounded by aluminum atoms should wait for 
hours or days to form molecules is not clear. 
However, these molecules were supposed to form 
gradually so that two particular aluminum atoms 
became firmly attached to each copper 
These molecules 


atom. 
or small groups of molecules 
were then supposed to diffuse to certain regions 
to form large groups at sites selected for future 
precipitation. This involved two difficulties 
First, it is difficult to visualize diffusion of such 
large particles through the aluminum lattice 
with the necessary velocity ; and second, it meant 
necessarily that these molecules were diffusing 
from regions of low concentrations to regions of 
high concentrations (a new inverse or up-hill type 
of diffusion which must be assumed in most 
of these pre-precipitation hypotheses). The 
breaking-up of the ordered structure to form 
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these groups of molecules would decrease the 
resistivity inasmuch as the ordered structure was 
defined as being a type which increased the 
resistivity. When sufficiently large agglomerates 
of molecules had been formed, they were trans- 
formed into crystallites, which were conceived as 
being in solid solution! Finally, these crystals 
become so large that the parent lattice could not 
stand the stress and they were then rejected from 
solution. The theory also included an interesting 
alternative method of precipitation. If large 
groups of molecules were adjacent, then the 
parent lattice between them could not stand the 
stress set up and was rejected from solution. 
Such a precipitate would have a composition 
between that of the phase to be precipitated (6’) 
and that of the solid solution. This theory has 
been reviewed because it is one of the most 
recent and because other pre-precipitation theo- 
ries have one or more of the faults illustrated by 
this theory. 

As I said before, people were forced to such 
theories because they believed that x-ray back 
reflection patterns had conclusively demonstrated 
the absence of significant precipitation. However, 
several years ago it was found” that a great deal 
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Fic. 7. Changes in electrical resistivity of a quenched 
aluminum-copper alloy which was aged alternately at low 
and high temperatures. [This figure was taken from an 
article by Willy Hartnagel entitled, ‘‘Die Kalt-Aushartung 
von Aluminium-Kupfer-Legierungen,” in Zeits. f. Metall- 
kunde 30, 84 (1938).] 


of precipitation could occur without any meas- 
urable change in lattice parameter. 
shows 


Figure 8 
a cast aluminum-magnesium alloy speci- 
men (10 percent Mg) which had been given a 
solution heat treatment at 425°C quenched in 
water and aged four hours at 300°C. No change in 
lattice parameter could be detected. Obviously 
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there are many precipitate particles clearly re- 
solved at this magnification—500 diameters. 

A similar situation the case of 
aluminum-copper alloys. The photomicrograph 
of Fig. 9, for example, shows precipitation along 
slip planes of 4 percent copper alloy which had 
been given a solution heat treatment and then 
aged for five minutes at 100°C. The precipitate 
particles are too small to be resolved at this stage, 
but the which precipitation has 
occurred are clearly revealed. No change could be 
determined in the lattice parameter at this aging 
temperature until the specimen had been aged 
for 32 hours. 


exists in 


regions in 


Once it had been discovered experimentally 
that lattice parameter measurements furnish no 
criterion for the absence of precipitation, it 
seemed surprising that they should ever have 
been considered as evidence. 
During the early stages of precipitation, most of 
the solid solution is unaffected. The initial 
precipitation involves those groups of copper 
atoms which existed prior to the quench. Further 
aging in alloys like the aluminum-copper alloys 
depletes the solid solution adjacent to the 
precipitate particle but does not change the 
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Fic. 8. Structure of a cast aluminum-magnesium alloy 
specimen (10 percent Mg) which had been aged for four 
hours at 300°C after solution heat treatment. 


concentration or parameter of the bulk of the 
solid solution until much later when diffusion, 
which is slow at these temperatures, has equalized 
the concentration. These regions of depleted solid 
solution are very small, have a variable compo- 
sition and are undoubtedly strained. Under such 
conditions several percent of the volume of the 


82 





alloy would have to be depleted before any 
change could be detected in a back reflection 
pattern. In other words, the parameter measure- 

















Fic. 9. Structure of a wrought aluminum-copper alloy 
specimen (4.00 percent Cu) which had been aged for 5 
minutes at 100°C after solution heat treatment. 


with 
theory of precipitation hardening. 


ments are not inconsistent the simple 

After it was shown that previous reasons for 
the assumption of pre-precipitation theories were 
untenable, an additional argument was brought 
forward. Briefly this argument was that time 
hardness curves obtained during the aging of 
binary alloys sometimes exhibited two peaks. 
The explanation suggested was that the first peak 
represented pre-precipitation hardening of some 
sort, whereas the second peak represented true 
precipitation hardening. Referring again to Fig. 2 
showing the aging of an aluminum-copper alloy, 
two peaks are clearly shown in the tensile 
strength and yield strength curve of this 4 
percent copper alloy quenched in cold water and 
aged at 100°C. However, at the time of the first 
peak, precipitation could be observed micro- 
scopically in the plastically deformed regions of 
the specimens; that is, along slip planes and 
grain boundaries. It seemed evident, therefore, 
that the two peaks are caused by different rates 
of precipitation in different portions of the 
specimen, the precipitation being much more 
rapid in the plastically deformed portions. This 
conclusion confirmed by quenching the 
sample in oil at 100°C so that little, if any, plastic 
deformation occurred. The aging curve obtained 
on these specimens quenched in hot oil is shown 
in the lower graph. In other words, double aging 


was 
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peaks offer no evidence against 
precipitation hardening theory.*! 


the simple 


It seems that the sooner we relegate to history 
all of these pre-precipitation theories (which 
were, after all, based on invalid assumptions) and 
return to the theory of precipitation hardening, 
the more rapid will be our progress in this field. 

Up to the present time we are justified in 
summarizing our knowledge of precipitation 
hardening as follows: 

(1) Precipitation hardening of alloys results 
from the controlled decomposition of a super- 
saturated solid solution, as originally postulated 
by Merica, Waltenberg, and Scott. 


(2) This decomposition is accomplished, as in 
other systems, by the formation of nuclei and 
their growth. 

(3) The nuclei form by chance grouping of 
solute atoms, as in other precipitation reactions. 

(4) Precipitation occurs in accordance with 
the Widmanstatten other 
crystalline systems. 

(5) Growth of the precipitated particles occurs 
by the usual process of crystal growth which 
involves the normal down-hill type of diffusion. 

(6) The changes in physical and mechanical 
properties are accounted for by this normal 
precipitation process. 


mechanism, as in 
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The Hardness of Metals as Affected by Alloying Agents* 


BY SIDNEY SIEGEL 


Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 


HE subject of hardening by alloying agents 

may be divided broadly into two sections. 
In one, the alloys exist as a single homogeneous 
phase, a solid solution. The solute atoms are 
distributed through the lattice on an atomic 
scale. In the other, the alloys exist as a con- 
glomerate of several phases. The added atoms are 
dispersed on only a microscopic scale in the alloy. 
Several of the most important aspects of this 
part of the subject, precipitation hardening, and 
the hardening of steels, are discussed in other 
papers of this series. 

As examples of these two possibilities, we have 
the Zn-Cd and Zn-Sn systems. Figure 1 shows 
the variations of the stress at the yield point for 
specimens of these two alloys systems.' Since Cd 
is soluble in Zn to about 1 percent at room tem- 
perature, while Sn is soluble to less than 0.05 
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Fic. 1. The variation in the critical shearing stress of 


Zn-Cd and Zn-Sn alloys with alloy content. 


percent, the Zn-Cd specimens are solid solutions, 
the Zn-Sn specimens are heterogeneous alloys, 
conglomerates of the Zn-Sn eutectic in the nearly 


* Presented at Conference on Hardening of Metals 
held at the University of Pennsylvania on May 17, 1941. 
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pure Zn matrix. The soluble Cd has a much 
larger hardening effect than does the insoluble 
Sn, dispersed as a coarse eutectic. 

Solid solutions are not, however, generally 
harder than two-phase systems of the same 
degree of alloying. In the Al-Cu system, an 
alloy containing 4 percent Cu can be quenched 
from the homogeneous field to retain the solid 
solution at room temperature for some time. If 
the alloy is aged at 200°C, its hardness first 
increases considerably as precipitation of the 6’ 
phase occurs, but aging for longer periods results 
in a decrease in hardness, as the precipitated 
phase agglomerates into large particles of the 
stable @ phase. A phenomenon similar to this is 
the decrease in hardness accompanying the 
“spheroidization” of steel, during which Fe3C 
dispersed as pearlite agglomerates into coarser 
aggregates. There appears to be a critical degree 
of dispersion in heterogeneous alloys at which 
maximum hardness is attained. If this degree of 
dispersion is far exceeded, then the hardness of 
the alloy will be determined primarily by the 
hardness of the solid solution matrix. 

This point is well exhibited in the Ag-Cu 
system which has two limited solid solution 
fields, with a eutectic in the two-phase region. 
Figure 2 shows the Brinell hardness of Ag-Cu 
alloys in the quenched and annealed states.* The 
hardness of the solid solutions is nearly the same 
for both heat treatments. At intermediate com- 
positions, however, the quenched alloy contain- 
ing finely divided eutectic is considerably harder 
than the annealed alloy. The quenched alloy 
behaves as a mixture of solid solution and 
eutectic, with a hardness about the average of 
those of the eutectic and the solid solution. The 
eutectic can here be considered to behave as one 
component of the system. The annealed alloy, in 
which the eutectic has been so coarsened that 
the alloy must be considered a mixture of two 
solid solutions, exhibits a hardness which is 
simply the average of the solid solution hardness 
values. 
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Before going further into a discussion of solu- 
tion hardening, it is important to consider the 
effects of the methods of measurement upon the 
data. In the examples cited above, two measures 
of hardness were given. In the work of Rosbaud 
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Fic. 2. The variation of the Brinell hardness of Cu-Ag 
alloys with Ag content, in the annealed and quenched 
conditions. 


the hardness is measured in terms of the critical 
shearing stress, the stress required to produce 
the first measurable departure from purely elastic 
behavior. This in itself is not a well-defined 
physical property of the metal, since it depends 
on the sensitivity of the apparatus, but at any 
rate it does limit the plastic deformation to a 
small quantity. In the work of Kurnakow, 
Brinell hardness was measured. Most of the 
data on hardness is obtained by the more usual 
tests, such as Brinell, Vickers, Rockwell, etc. 
In these tests a load is applied, and the amount 
of plastic deformation produced by the indenter 
is measured after the load is removed. The 
hardness is then given either on an arbitrary 
scale, or in terms of the ratio of load to area of 
indent. In such measurements, not only the 
initial ability of the material to resist plastic 
deformation, but also its work-hardening ca- 
pacity determine the size of the indent and hence 
the hardness values. This point is brought out in 
Fig. 3, which shows the stress-strain curves for 
an alloy having the composition Cu;Au, after 
different heat treatments.* The critical shearing 
stress for Curve A is higher, and on this criterion 
the quenched alloy is harder, but if the conditions 
of the test are such that considerable plastic 
deformation occurs, the crossing point of the 
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two curves may have been passed, and such a 
test would show the annealed material to be 
harder; not because of its greater initial hard- 
ness, however, but merely because of its greater 
work-hardening capacity. 

For this reason it is desirable in seeking to 
increase our understanding of the problem of 
solution hardening, to use a simple test such as 
critical shearing stress, rather than a complex 
test such as indentation hardness, in comparing 
different alloys. The effect of work-hardening on 
the ball-indenter hardness values may be taken 
into account if measurements are made at differ- 
ent loads and the data subjected to Meyer 
analysis,‘ but this is rarely done in most of the 
technical work. 

From the technical point of view, it is just 
this work-hardening capacity which makes homo- 
geneous alloys more useful than pure metals. 
This may be seen from Fig. 4, which shows the 
Vickers hardness of a series of Cu-Zn alloys in the 
annealed and cold-rolled states.' In the annealed 
state the hardness of 70:30 brass is 55 as 
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Fic. 3. The strain-hardening curves for single crystals of 
CusAu in the ordered and disordered states. 


compared with 40 for pure Cu. The solution 
hardening effect is thus rather slight. In the 
cold-rolled state, Cu has a hardness of 110, 
while 70:30 brass has the value 230. Thus 
70 : 30 brass can be hardened by a factor of 4, 
while pure Cu only by a factor of 2.5, and it is 
this greater work-hardening capacity of the 
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alloy which is important, not the small increase 
in intrinsic hardness of the annealed solid solu- 
tion. This subject is discussed in detail in one 
of the accompanying papers, and the work- 
hardening capacity of alloys has been mentioned 
here merely because of its important bearing on 
the measurement of the hardness of alloys. 


SOLUTION HARDENING 


lor the practicing metallurgist the most im- 
portant means of hardening metals are those 
which depend on transformation, as in the steels; 
on precipitation, as in Al alloys; and on work- 
hardening capacity, as in the brasses. From the 
physicist’s standpoint, however, a clearer under- 
standing of solid solution hardening is needed 
before we can develop our ideas concerning 
heterogeneous alloys. 

In the case of a solid solution we are con- 
cerned with a single-phase alloy and the changes 
in hardness observed are due to the matrix 
alone, not to the presence or degree of dispersion 
of another phase. Solution hardening is thus a 
true property of the lattice. 

The Ag-Au system is an example of an ideal 
solid solution. Figure 5 shows the critical shearing 
stress of a series of these alloys. The 50-50 
atomic percent solution is about ten times as 
hard as the pure metals. This curve is more or 
less typical of other complete solid solutions, 
but not all are as symmetrical as this. The 
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Fic. 4. The Vickers hardness of a-brasses in the annealed 
and cold-worked conditions. 
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Fic. 5. The critical shearing stress of Au-Ag single crystals 
as a function of silver content. 


number of systems which exhibit complete solid 
solubility is rather small, but a large number of 
alloys which have limited solid solubility exist. 
In general it is found that for small concentra- 
tions of the solute the increase in hardness varies 
nearly linearly with the amount of alloying 
element. 

Figure 6 shows the increase in vield strength 
of Fe-Si and Fe-Al alloys with increasing amounts 
of Si or Al.’ An initial linear relationship is 
observed for both alloys, but the hardening 
effect per atom is considerably higher for the 
case of Si than for Al. 

Norbury carried out an extended investigation 
of the Brinell hardness of a@ solid solutions of 
Cu.* He used the Meyer analysis on his data to 
correct for the effect of work-hardening. A nearly 
linear relationship between the initial hardness 
and atomic percent of alloying element was 
found for all the solid solutions. The slope of 
the curve varied from one alloy to another, 
however, and in Fig. 7 the increase in hardness 
per atomic percent solute element is plotted 
against the difference in atomic size of the 
solute element and Cu. The atomic size in each 
alloy was obtained from the measured density of 
the alloy and the fact that all have the face- 
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centered cubic structure. It is observed that the 
larger the difference in atomic size, the larger 
the hardening effect per atom. 


TABLE I. Lead a solid solutions. 


Added Increase in hardness Solubility 
metal per atomic % limit 
Bi 0.5 30% 
Sn 1.8 10 
Cd 4.8 5 
Sb 10.0 2m 
Mg 11.0 0.5 


Ludwick measured the hardness of various Pb 
base solid solutions, and observed a correlation 
between the hardening effect per atom and the 
extent of solubility. His data® are summarized 
in Table I. 

In the solid solutions discussed so far the 
solute atoms are distributed at random on the 
sites of the lattice of the solvent pure metal. 
In certain solid solutions it is possible to have an 
ordered arrangement of the atoms on the lattice 
sites,'!° and in fact by appropriate heat treat- 
ment it is possible to control the degree of 
ordering and thus control the nature of the 
distribution of atoms in the lattice sites. This 
order-disorder transformation occurs in the Cu- 
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Fic. 6. The variation of the yield-point stress of Fe-Si 
and Fe-Al alloys with Si or Al content. 
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Au alloys near the concentrations Cu;Au and 
CuAu. Above 425° the alloys of all compositions 
are face-centered cubic. When in equilibrium, 
CuAu becomes ordered below 425° and also 
becomes slightly tetragonal in structure, while 
Cu3Au is ordered below 390°. but remains cubic. 

The Brinell hardness of these alloys has been 
measured at room temperature, in the annealed 
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Fic. 7. The change in Brinell hardness for one atomic 
percent added element plotted as a function of the differ- 
ence in size of the solute and solvent atoms. 


and quenched states.'! The annealed alloys are 
ordered, the quenched are random. The data 
appear in Fig. 8, from which it is to be con- 
cluded that ordering Cu3Au decreases its hard- 
ness, while the converse is true for CuAu. It is 
somewhat difficult to interpret these curves 
unambiguously, for they may be complicated by 
strain hardening effects. A more detailed exami- 
nation of Cu3Au is shown in Fig. 3. The critical 
shearing stress for the ordered state of the alloy 
is 2.1 kg/mm?’, and is raised to 4.4 kg/mm? if 
the alloy is disordered. The difference in strain 
hardening capacity of the two materials is 
evident from this curve, as was pointed out 
above, and it is this which makes the data of 
Fig. 8 uncertain. 

Thus far only a@ solid solutions have been 
mentioned. Many alloy systems show homo- 
geneous intermediate phases over a finite range 
of compositions, and may also form intermetallic 
compounds at definite compositions. 

The 8-phase, having a body-centered cubic 
structure is found in many alloy systems. The 
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alloys in the annealed and quenched states, plotted against 
Au content. 


critical shearing stress of 8-brass® is 10 kg/mm”, 
while 70:30 a-brass single crystals have a 
critical shearing stress of 1.5 kg/mm*, consider- 
ably lower. In general, the 8-phases are harder 
than the a and other inter- 
mediate phases with more complex structures, 


solid solutions, 


such as the y-phase in brass, are still harder. 
The hardness-composition curve for the Cu-Zn 
system appears in Fig. 9. 

In other systems definite chemical compounds 
are found, the most common example of which 
is FesC. In general, the hardness of these inter- 
metallic compounds is very high, and they do 
not have any appreciable plastic range, but show 
brittle fracture. 


SUMMARY 


In heterogeneous alloys, the hardness is de- 
termined mainly by the structure. The nature 
and degree of dispersion of the hardening phase 
determine the hardness of the alloy as a whole. 
For a solid solutions, several generalities can be 
stated. First, the solutions are always harder 
than the pure metals. Second, in dilute solutions 
the hardness increases linearly with the number 
of solute atoms. Third, the hardening effect per 
atom is greater the more unlike the solute and 
solvent atoms are. Fourth, a disordered solid 
solution is harder than the same alloy in the 
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ordered state, if no change in crystal structure 
occurs. In intermediate solid solutions the crys- 
tal structure is generally more complex and the 
alloys are harder, with little or no plasticity. 
Intermetallic compounds are extremely hard and 
brittle as a rule. 


The metals deform plastically by slip on 
definite crystallographic planes, in definite crys- 
tallographic directions. It has been realized for 
some time that the observed critical shearing 
stresses could not be explained in terms of a slip 
mechanism in which one part of a crystal slipped 
as a unit past another part on the other side of 
the slip plane. The present view is that a certain 
type of imperfections, called “dislocations,” exist 
in the lattice and slip is produced by the motion 
of these dislocations through the crystal." These 
dislocations move under the action of the ex- 
ternally applied stress, and any property of the 
lattice which impedes their motion makes the 
lattice harder. 

On the basis of such a picture, Taylor has been 
able to develop a theory of strain hardening, 
and it appears that at least the fundamental 
ideas of the theory are correct. In a qualitative 
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Fic. 9. The Brinell hardness of the Cu-Zn system. 
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way we may also describe solution hardening 
effects in terms of this theory. The solute atoms 
interpose barriers to the motion of the disloca- 
tions, and thus a larger stress is required to 
enable this motion to continue through lattice, 
which is necessary in order for slip to occur. 
In sufficiently dilute solutions, the number of 
stranger atoms which interfere with the motion 
of a given dislocation increases nearly linearly 
with the number of dissolved atoms, so we 
would expect the hardening to be proportional 
to the atomic percent of solute. Also, the barrier 
interposed to the motion of a dislocation is 
determined by the distortion produced in the 
lattice by the dissolved atom, and so would 


depend on the size difference between the solute 
and solventatom. The difference between ordered 
and disordered solution is due to the fact that 
it is the deviations from perfect periodicity in 
the lattice which impede the motion of disloca- 
tions through it, and in the disordered state 
these deviations are greater than in the ordered 
state. 

In the intermediate phases, slip planes are not 
as well defined as they are in more simple struc- 
tures, and so it is more difficult for slip to occur. 
The intermetallic compounds in general do not 
behave as metals at all since the type of binding 
is different, and in general they exhibit brittle 
failure. 
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Chapter VII. Viscosity of Colloidal 
Solutions. Application to 
Synthetic Polymers 
Abstract 


OLLOIDAL particles present in solutions 
and suspensions increase the effective vis- 
cosity of the solvent. Einstein has presented the 
theory of this effect valid for spherical particles, 
and given an equation, deviation between which 
and experiments being due to either solvation or 
electric charge of the particles. Elongated par- 
ticles require a further modification of the 
formula. The theory is being successfully applied 
to the study of synthetic fiber molecules, and 
valuable information as to the shape of the latter 
can be derived by viscometric measurements on 
colloidal solutions of such polymers. 


28. Introduction. Einstein’s Equation 


According to the principle of this monograph, 
of proceeding step by step from simple to more 
complicated systems, colloidal solutions and 
suspensions should be dealt with next, homo- 
geneous liquids having been considered in the 
two preceding chapters. The addition of col- 
loidal or suspended particles to the liquid phase 
causes an increase in viscosity. In this chapter 
the quantitative relations controlling this in- 
crease will be presented, followed by an applica- 
tion of these relations to certain problems of 
scientific and technical interest. Since the 
mathematical treatments are too involved to be 
included, only some approximate considerations 
and the results of computations will be given. 

Colloidal systems can be divided into sols and 
gels. Sols are those systems in which the disperse 
phase is subdivided into small, separate particles, 
floating in the dispersing medium. Gels are 
in which the disperse phase forms a 
coherent network, soaked in the liquid dispersing 
medium. Gels evidently are more complicated 
systems than sols. The present chapter will be 
concerned with sols only, gels being treated in 
the following chapter. 


systems 
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The name sol is used here in the widest sense 
given by the definition above, and independently 
of the size of the floating particles. Actually, 
three further groups can be differentiated within 
the main heading. 

First, there is a dispersion of the smallest 
particles, up to about 10~® cm, invisible even by 
the ultramicroscope (but now within the range 
of the electron microscope'). This type of sol is 
called hydrophilic, or more generally lyophilic 
sol, indicating a strong attraction between par- 
ticle and solvent. Recently they are also called 
macromolecular systems, implying that the par- 
ticles are really nothing other than molecules, 
only of considerable size, and that the systems 
are real solutions in the sense of the kinetic 
theory. Typical examples are dilute solutions of 
gelatin in water, or polystyrene in benzene. 

Second, there are the next larger particles, 
visible by the ultramicroscope, generally termed 
hydrophobic or lyophobic sols, and having a 
size from about 10~® to 10~* cm. They consist of 
minute molecular aggregates, mostly crystal- 
lites. A typical example is a colloidal gold sol. 
There is no intimate connection between the 
particles and the solvent, and they do not form 
real solutions. 

The third group contains the largest particle 
size, above 10~* cm. They are called suspensions, 
if the disperse phase is solid, emulsions, if it is 
liquid. Examples are a kaolin suspension and an 
oil emulsion in water. Like the second group. 
they are two-phase systems with one of the 
phases subdivided. 

The following deductions on viscosity are 
strictly valid for such two-phase systems as 
groups 2 and 3. The application of the same rules 
to macromolecular solutions (group 1) is to a 
certain extent an extrapolation but, judging 
from its experimental verifications, appears to 
be a safe working hypothesis. 

The question now arises, what is the reason 
for the observed increase of viscosity in colloids 
as compared with that of the dispersing medium 
alone? This is accounted for as follows: The 
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presence of the particles will set up certain new 
boundaries within the liquid, thus altering the 
distribution of the stream lines. It is generally 
assumed that the slip coefficient [Eq. (22) ] is 
zero; thus each single particle will create new 
boundary conditions, according to which the 
new system of flow lines can be computed. This 
is a purely hydrodynamical problem. 

Thus it can be seen that the viscosity coef- 
ficient » of the liquid phase itself remains unal- 
tered. It is the distribution of flow lines, and 
thus the components of the deformation velocity 
[ Eq. (1) ], that change, even though the external 
stresses at the boundary of the total system 
remain unaltered. As a result of these changes 
the ratio mean velocity : mean stress will be 
reduced at the boundaries of the whole system, 
as will be shown in a simple case. 








Fic. 44. Influence of a single suspended particle on the 
effective viscosity of a liquid. 


Let us assume (Fig. 44) a parallel flow u in the 
x direction with a positive gradient (du /dy) =a 
in the y direction. A small spherical particle of 
radius r is situated in the plane of symmetry 
(y=0 and u=0). In Section 4 this case was dis- 
cussed; it was mentioned that the resulting 
torque on any volume element is zero and thus 
that the rotation of volume elements must occur 
at a constant angular velocity. It can be shown 
that in case of a uniform laminar flow each 
volume element rotates with an angular velocity 
u/2 and evidently the same rotation will be 
exhibited by the particles. 

This means that the velocity of the liquid at 
the point A will be r#/2 as the liquid adheres to 
the surface of the particle. Since the original 
velocity (without particle) at this point was ra, 
it is now reduced to half of this value. This 
retardation of the motion evidently cannot be 
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restricted to one point only, but will extend 
further in the positive y axis. The same holds, of 
course, also for the negative y axis, where the 
direction of the motion is reversed. As a result, 
the velocity even at the boundary, B, (y=h/2), 
will be less than it was originally. The value at 
the boundary will vary with the distance x from 
B. In computing average values along the whole 
boundary, it will be found that, while the average 
of u is reduced compared with the original value, 


namely to 
hi 107 
u= “(1- —"rn) 
2 3 


(where is the number of particles per unit 
volume), the average stress at the boundary is 
unaltered. 

It is now possible to consider the system as 
having an equivalent viscosity 7’ of a homo- 
geneous liquid. Using the fundamental Eq. (6) 
between stress and gradient, one has for the 
stress in the original liquid na. In the equivalent 
homogeneous liquid the average stress is given 
by »/a(1—(107/3)r'n). Since both stresses, as 
mentioned above, are equal, it follows that 


107 
f=n(1+—"rn), 
3 


since the second term is small compared to. unity. 
Denoting the volume of one particle by v, Eq. 
(83) becomes 


(82) 


(83) 


n’ = n(1+2.5nv) 


as originally derived by Einstein. Another fre- 
quently used form of the equation is obtained by 
denoting the fraction (n’—n)/n as n.»=specific 
viscosity, and writing c for the concentration of 
the particles in g/cc: 


(84) 


Nsp = 2.5¢/p, (85) 


p being the density of the particle material. 
The validity of the Einstein equation has been 
proved for a number of lyophobic sols and sus- 
pensions of various kinds. It is safe to assume 
that deviations from the theoretical form 
originate in some sort of departure from the 
original assumptions. Two of the most frequent 
departures, namely solvation and electric charge, 
will be considered in this section; a third one, 
namely the departure from the spherical in the 
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shape of the particles, in the two following 
sections. 


It is often found that »,, varies with the 
solvent, a fact not accounted for by the theory. 
At the same time, results obtained from certain 
sols, particularly of the lyophilic type, often 
indicate a coefficient in Eq. (85) larger than 2.5. 
The explanation of these facts is that a certain 
amount of solvent is closely associated with the 
colloid particles, thus causing its apparent volume 
v to be larger than that computed on the basis 
of concentration and density of the solid. This 
solvation process easily explains the deviations, 
and it is even possible to use the equation in 
order to evaluate the degree of solvation. Figure 
45 indicates the trend of »,, for agar sols if (for 
a given concentration of colloid) the solvent is 
gradually changed from water to alcohol.? The 
explanation above is further strengthened by the 
fact that, while the originally aqueous system 
with a high specific viscosity is stable, it becomes 
rather unstable at higher alcohol concentrations 
and lower specific viscosity. This is because 
solvation is one of the chief stabilizing factors in 
colloidal systems. 

There are other cases in which the behavior of 
Nesp cannot be explained on the above basis alone. 
An explanation can be given in considering the 
Smoluchowski extension of Eq. (85) in the form: 


2.5¢ ; se" 
Nep= 1+ ( -) , (85a) 
p or*n\2r 


where o and e=electrical conductivity and dielec- 
tric constant of the medium, ¢=potential dif- 
ference at the boundary particle-solvent. Ac- 
cording to this extension, an electric charge (¢) 
increases the viscosity effect of suspended par- 
ticles. Figure 46 shows the behavior of an agar 
sol with increasing concentration of added 
electrolytes.” It is known that the electric charge 
of negative sols decreases with increasing con- 
centration and valency of the cations present; 


therefore, according to Eq. (85a) the specific 


viscosity must also decrease. 

Among the large number of experimental 
papers on Eq. (85) a recent one by Mardles* 
should be mentioned. He took viscosity data of 
kaolin and graphite suspensions in various 
organic liquids. In certain cases agreement with 
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Fic. 45. Variation of the specific viscosity of an agar 
sol with varying mixture ratioof the solvent (water-alcohol), 
after Bungenberg de Jong. 


Eq. (85) was obtained but, in others, viscosities 
four times higher than those indicated theoreti- 
cally were measured. Discrepancies were found 
to be due to interaction of particles, not con- 
sidered in the derivation. 


29. Elongated Particles 


It will be of interest to consider how the 
fundamental Eq. (85) is altered if the particles 
in suspension are elongated rather than spherical. 
This problem has been treated by Jeffery’ and 
Eisenschitz® among others, and a good review 
has been written by Burgers.® Elongated par- 
ticles, which can be treated either as ellipsoids 
of large axis ratio or as cylindrical bodies, are of 
considerable importance for the study of a 
number of organic plastic materials consisting 
of thread-like molecules, to be discussed in the 
next section. 

Accurate treatment of such particles is fairly 
involved, in that the influence upon the effective 
viscosity n’ depends upon their orientation with 
respect to the flow lines. Because of the rotatory 
motion of volume clements in the case of a 
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streamlined flow with finite gradient, the long 
particles should assume a direction parallel to 
the flow, as soon as the stationary condition is 
reached. However, the particles will perform a 
Brownian motion, characteristic of colloidal 
particles, which is the exact analogue to the 
thermal motion of molecules, expressed by Eq. 
(17) in Section 8. This motion will tend to keep 
the particles in a state of random orientation, 
distributed almost equally among all the possible 
angles in space. 

Since 7’ and thus 7,, depend upon the state of 
orientation, it is important to have an idea of 
this latter. It will ultimately depend upon the 
degree of predominance of the shearing motion 
as compared with the Brownian motion. A 
complete analysis shows that the situation will 
be controlled by the condition 


uR/2kT=1, (86) 


where @ is the shear gradient and R the viscous 
resistance of the liquid against rotation of the 
particle. For cylindrical particles we have 
arnl? 
R=—————_ - (87) 
3(log 1/r—0.80) 

(/=length, and r=radius of cylinder). If the left 
side of (86) is large compared to unity (shear 
prevailing), the particles will assume positions 
slightly deviating from the direction of flow. If, 
on the other hand, the expression in (86) is 
small compared with unity (Brownian motion 
prevailing), the orientation is essentially random. 

The theory shows that the specific viscosity in 
these cases is given by an equation similar to 
Eq. (85), with the difference that the factor 2.5 
has to be replaced by a more general expression, 
namely : 

fP/4r? 
an emremrnnn (88) 
log (l/2r) —0.8 
where f is a dimensionless factor, depending upon 
the degree of orientation. 

If, according to (86), Brownian motion still 
prevails, and the orientation is more or less 
random, f very nearly assumes the value 1/15. 
If, on the other hand, according to (86), shearing 
motion is prevalent and the orientation of the 
particles becomes gradually parallel, the factor 
becomes proportional to (2k7T/a#R)'; hence, 
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because of Eq. (87), approximately proportional 
to t/l. In the extreme case of very high shear, 
the particles are all parallel to the direction of 
the flow and, according to the expression for f, 
Nsp approaches the value zero. The viscosity is 
then not affected by the presence of the particles. 

It was mentioned in the introduction to this 
chapter that its contents refer to sols and not to 
gels. The underlying idea is the assumption that 
in sols the individual particles are independent 
of each other, whereas in gels the mutual influ- 
ence has reached the extreme stage of forming a 
network. Equation (85) and also the subsequent 
ones are deduced on the basis of a complete 


independence of individual particles, as is 


° 
° 





wo 
°o 


70 











++++ 


Complex Pt 


SPECIFIC VISCOSITY (SOLUTION WITHOUT SALTS+I00) 
@ 
° 


60 














! 2 3 4 
NORMALITY OF ELECTROLYTE 


Fic. 46. Variation of the specific viscosity of an agar 
sol with concentration and valency of added cations, 
after Bungenberg de Jong. 


evidenced by our having introduced the number 
of particles n simply as a factor in the second 
term of Eq. (82), originally deduced from con- 
sideration based on one single particle. This re- 
striction can be formulated by saying that all 
these equations refer, strictly speaking, to low 
concentrations such that 


nv<l, (89) 
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which means that the volume of a particle should 
be small as compared with the space at its 
disposal (1/). The less this condition holds, the 
more it is likely that neighboring particles will 
influence each other by means of collisions and 
also attraction forces. This will necessitate the 
introduction of other terms in the equations. 


Finally, when the concentration becomes large 
enough for the particles to form a_ network, 
radical changes in the viscous behavior of col- 
loidal systems occur, to be discussed in Chapter 


VIII. 


30, Application of the Theory of Viscosity 
of Colloids to Polymers of Fiber 
Structure 


In recent years the considerations just dis- 
cussed have gained in importance, because of 
their application to the study of a number of 
synthetic plastics. This field in chemistry, de- 
veloped by Carothers,’ Staudinger,® and others, 
has also considerable industrial significance. It 
is concerned with materials obtained synthet- 
ically from simple compounds by means of 
polymerization. The original material is called a 
monomer, in contradistinction to the polymer, 
the final product. The monomeric molecules are 
usually joined together in thread-like, so-called 
“fiber” structures, and a number of their physical 
properties can be explained on this basis. The 
number of monomers present in one polymer or 
macromolecule is an important characteristic of 
the material, since it determines some of the 
physical and chemical properties. This figure is 
sometimes denoted as the degree of polymeriza- 
tion. If it is designated by p, the molecular 
weight of the monomer by Mo and that of the 
polymer by /, we evidently have the relation: 


M = pMo. (90) 


Thus a determination of the molecular weight 
M is equivalent to a determination of p, since 
Mo is known in each case. 

There are two methods available for the 
measurement of M, namely the osmotic pressure, 
or cryoscopic method, and the ultracentrifuge 
process. Both are fairly reliable, but have their 
limitations. The cryoscopic method, for instance, 
is used chiefly for products of relatively low 
molecular weight. Consequently it was thought 


04 


desirable to have available a third method, to be 
used without restriction in the whole range of MV. 
Such a method was found by Staudinger in the 
viscometry of solutions of polymers. 

The polymers are generally soluble in organic 
solvents like benzene, chloroform, and the like, 
and the viscosities of such solutions are, of 
course, increased as compared with the value for 
the solvent. It is usual to operate with the 
specific viscosity 7.» as defined in Section 28. 
Staudinger has found that, by using different 
polymers (varying 7) of the same monomer, the 
specific viscosity for a given concentration 
increased with the molecular weight J/. If a 
definite relation between the two could be 
established, 1/ could be computed from the vis- 
cosity of the solutions of the polymer. 

Let us first see what relation between »,,/c¢ 
(specific viscosity per unit concentration) and M/ 
might be expected on the basis of expression 
(88) of the last section. 


TABLE VI. Viscometric determination of molecular weights 
of nitrocellulose samples, after Staudinger. 


CALCULATED FROM 7 OBSERVED 

Msp. Ea. (93) CRYOSCOPICALLY 
59 54,000 50,000 
82 75,000 82,000 
210 190,000 178,000 
450 410,000 443,000 


TABLE VII. Viscometric determination of molecular weights 
of cellulose triacetates, after Staudinger. 





M 

CALCULATED FROM OBSERVED 

nsp/c Ea. (93) CRYOSCOPICALLY 
4.7 47,000 45,000 
5.5 55,000 54,000 
6.2 62,000 60,000 
12.0 120,000 112,000 
17.0 170,000 155,000 
23.3 230,000 275,000 


If the molecules are supposed to be rod-shaped, 
then r can be taken as being independent of the 
length and thus as a constant. The length / is 
then proportional to both p and AM. It can be 
seen that for strong Brownian motion the factor 
in (85) becomes proportional to ?/r*, since the 
variation of the logarithmic term in the de- 
nominator is slow. With r=constant, it follows 
that approximately 


nep/¢» M2. (91) 
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With, more pronounced shear this relation 
must change since f diminishes and becomes 
proportional to 1//. Consequently the approxi- 
mate relation 

Nsp/ 6 M (92) 
should hold. 

Now, let us consider the experimental findings. 
Staudinger finds an empirical relation of the 
form: 


Nep/C=KM. (93) 


In Eq. (85) c was the concentration in g/cm', but 
Staudinger prefers to use the molarity (in terms 
of the monomer) per liter. This difference entails 
only a different value of the constant K. 

In Table VI, Staudinger’s data on_nitro- 
cellulose are shown, the solutions being in 
acetone. Four different degrees of polymerization 
were selected and the experimental figures for 
N«p € are shown in the first column. The ratio 
between 7s, and ¢ is not strictly linear on the 
whole concentration range; consequently, the 
ratio given is the limiting value for small concen- 
trations. In selecting the value K=11X10~ the 
molecular weights .7 are computed, using Eq. 
(93). In the third column are the molecular 
weights obtained from osmotic measurements, 
and the agreement can be seen to be satisfactory. 

Another series of data is shown in Table VII 
for cellulose triacetate with chloroform as solvent, 
the degree of agreement being about the same. 
K is here 1.0X10™. 

In looking at these data, it seems that Eq. (93) 
should hold fairly well. Theoretically, one could 
conclude that the data refer to the range of pre- 
dominating shear, characterized by Eq. (92). 
Unfortunately such a picture would show the 
situation over-simplified; it has very definite 
limitations, viz., the Staudinger equation (93) 
does not hold in all cases; in fact, the deviations 
of observed from computed values are often 
large. This is not surprising, since the shear 
values used in the experiments are in reality not 
so large as to cause shear to prevail. Thus either 
Eq. (91) or, still better, its original form (88) 
should be used, since the variation of the logar- 
ithmic term cannot well be neglected. 

Thus, while it seems that the simple and 
attractive Staudinger relation had to be dis- 
carded (generally, at least), one can make use of 
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the experimental data for n,, and see what con- 
clusions can be drawn if (88) is applied. It will 
be seen that the conclusions are interesting, and 
they do lead us one step further in the under- 
standing of the structure of polymers. 


TABLE VIII. Evaluation of viscometric measurements on 
polystyrene solutions, according to Burgers. 


M nep/C l/2r LinA 2r,1InA 
519 0.38 6.4 40 6.2 
2,460 0.73 12.7 105 8.3 
5,220 1.13 17.9 169 9.5 
7,600 1.41 20.9 213 10.2 
13,000 2.288 29.2 318 10.9 
30,000 te 32.6 452 13.9 
80,000 5.6 51.8 855 16.5 
270,000 24.0 18.6 


120.0 2270 


A computation of this type, as carried out by 
Burgers,® follows. It is concerned with polysty- 
rene, the polymer of monostyrene CH(CsH5)CHg; 
the polymer being soluble in benzene, tetralin, 
etc. (Addition of mineral oil to these solvents 
causes a transition lyophilic sol suspension.°) 
Determinations of the molecular weight were 
carried out by Staudinger cryoscopically, and by 
Signer!’ by means of the ultracentrifuge. These 
data are shown in the first column of Table VIII. 
The second column indicates values of ».)/c, as 
measured by Staudinger. It can be seen that Eq. 
(93) does not hold for this case, M increasing in 
the ratio 520 :1, and »,,/c only in the ratio 
63 : 1, in other words »,p/¢~ M$. 

By applying (88) with f=1/15, corresponding 
values of //2r have been computed, shown in the 
third column. Now from M we have the weight 
of one single molecule, given by M/(6X 10") and 
in dividing the latter by 1.08, the density of 
polystyrene, we have the volume of one molecule. 
Thus we have the relation M/(1.08X6 X10”) 
=rir*, and 1 and r then can be computed, the 
values being shown in the last two columns of 
Table VIII. 

Thus, the striking result is obtained that the 
thickness of the molecule is not constant, but 
increases with the molecular weight appreciably. 
An explanation of this behavior is supplied by a 
theory of Kuhn" on the statistical bundle-shape 
of polymer molecules. The theory is based on the 
assumption that the unit molecules are loosely 
joined together; and that this loose connection 
allows the shape of the polymer to be twisted as 
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a result of the independent motions and rotations 
of the single units within the whole molecule. 

It is assumed that the angle between two 
neighboring units, the so-called valency angle, 
has a given value (110°) but, apart from this 
restriction, free rotation of the units is possible. 








< | eff 
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Fig. 47. Projection on a plane of a bundle-shaped 
polystyrene molecule, according to Kuhn, 


While the actual position of the single members 
is unknown, the laws of probability can be applied 
to molecules containing a large number of units. 
The projection upon a plane of the three-dimen- 
sional bundles is shown in Fig. 47. The length of 
each unit is J», and the total, or effective length 
loc is evidently less than plo, if p is the degree of 
polymerization. The total length can be found 
by applying the laws of probability, and it is 
found that 

lese= (8 W) iby Pp, (94) 
where ) depends upon the angle of valency. For 
110° 

b= (2) 4p. 


The effective length is proportional to the square 
root of the degree of polymerization. The effective 
cross-sectional area of the bundle then increases 
with p. These conclusions are born out (quali- 
tatively, at least) by the data presented in Table 
VIII, showing that the increase of / is less than 
proportional to the increase of M. 

‘Kuhn further extends his theory to explain 
certain optical and mechanical properties of 


. solutions of polymers. In using the conception 


for solids, a sound explanation of the elastic 
properties of rubber is arrived at. 

The theoretical evaluation of viscometric 
studies on polymers has thus proved rather 
fruitful. However, as stated before, the simple 
linear relation suggested by Staudinger, cannot, 
at least generally, be accepted as valid. Kuhn 


96 





applies his theory to the hydrodynamical theory, 
the results from which are presented in Section 
29, and concludes that »,.,/¢ ought to be propor- 
tional to a power of M between 0.5 and 0.9. The 
data in Tables VI and VII and in Table VIII 
are in good agreement with this conclusion. 

A different empirical relation between the vis- 
cosity and chain length of polyesters was found 
by Flory” who claims that this relation can be 
used for accurate determinations of the molecular 
weight. The relation is a linear one between the 
logarithms of viscosity and the average chain 
length if a mixture of several chain lengths is 
present. 

The relation between the molecular weight and 
shape (/,r) of molecules on the one hand and the 
viscosity of their solutions on the other has also 
been used for the study of proteins. Polson! con- 
siders ellipsoidal protein particles of an axial 
ratio varying from 2.8 to 16.6. A good review of 
this matter is presented by Mark and Simha,!"' 
Burgers'® includes oblate rotational ellipsoids in 
his treatment, and also systems of rigidly con- 
nected spheres. 

The Staudinger equation is used also for the 
study of rubber. Kemp and Peters!® compare 
molecular weights of crepe rubber fractions 
obtained by means of viscometric and cryoscopic 
measurements. The two sets of data are not in 
complete agreement. 
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The Radiation from Carbon Arcs* 


By H. G. MACPHERSON 


Research Laboratories, National Carbon Company, Inc., Cleveland, Ohio 


EW people aside from those engaged in the 

industry realize the present commercial 
importance of the carbon arc. In fact most of us 
base our opinion of the arc on memories of the 
faltering light of the obsolete arc street lights and 
therefore think of it as an outmoded method of 
illumination. This general ignorance concerning 
the present status of the carbon arc is perhaps a 
tribute to its smoothness and dependability, for 
we all frequently watch the light from an arc for 
several hours at a time and are aware only of 
having seen a motion picture. Every theater uses 
the carbon arc as a light source with which to 
project the motion pictures on the screen and, 
as a matter of fact, the optics of the projection 
system are such that the spot of illumination on 
the screen is a section of a greatly enlarged and 
slightly out of focus image of the arc crater. The 
high uniformity of illumination and absence of 
flicker obtained by most motion picture houses 
is thus a direct indication of the smoothness of 
burning of modern carbon arcs. 

Since the carbon arc is the brightest com- 
mercially available light source, it finds usage 
today in any field where a very intense source of 
light is needed. The brightness of the arc and the 
white color of the light make it the most satis- 
factory source for the theatrical projection of 
motion pictures. Other typical fields where the 
carbon arc is important include photo-engraving, 
the irradiation of milk, the accelerated testing 
of the resistance of paints and other materials 
to photochemical action, and the _highest- 
powered searchlights of the Army and the Navy. 

The carbon arcs used today are of three 
general types. In the order of their natural his- 
toric development, they are the low intensity 
arc, the flame arc and the high intensity arc. All 
classes of arcs can be operated either on alter- 
nating or direct current. However, since in their 
most important uses the low intensity and the 


* Presented at the Symposium on Light and Lighting of 
the Ohio Section of the American Physical Society, Cleve- 
land, Ohio, October 10-11, 1941. 
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high intensity arcs are operated on direct current, 
we will speak here about direct-current arcs, 
except as otherwise noted. 

The so-called low intensity arc is an improved 
descendant of the original carbon arc discovered 
by Sir Humphrey Davy in 1801. This is an arc 
between plain carbon electrodes or electrodes to 
which certain salts have been added to a soft 
carbon core for purposes of stabilizing the arc. 
This “arc supporter’ does not affect the light 
emitted, so that the principal source of light of 
the low intensity arc is incandescent solid carbon 
at or near its sublimation temperature. 





Fic. 1. Low intensity arc. 


Figure 1 shows a side view of a 25-ampere 
55-volt low intensity arc. The positive carbon is 
the larger one and has a diameter in this case of 
about one-half inch. The positive carbon is 
always used as the light source in a direct current 
arc because it is the brighter electrode and has 
the larger incandescent area. 

Since the light from the low intensity arc 
comes by virtue of the incandescence of the 
anode surface, to understand its light emission 
we must know the reasons why the anode 
becomes heated and what determines its tem- 
perature. The positive carbon tip is heated to 
its high temperature by the discharge of the 
electrical current there. Part of this is accom- 
plished by absorption of energy from the elec- 
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trons entering it and part by conduction of heat 
from the adjacent are stream. The gas in the arc 
stream is extremely hot, with a temperature in 
the neighborhood of 6000°K, and has a relatively 
high ion density, enabling it to conduct the 
current easily. The current is carried through the 
arc stream largely by the electrons, since they 
move much faster than positive ions because of 





Fic. 2. Alternating current flame arc. 


their small mass and consequent high mobility. 
However, there is an equal density of positive 
ions interspersed throughout the arc stream so 
that it has no net space charge. The only re- 
sistance to the motion of the electrons is that 
supplied by the inert atoms and molecules with 
which they frequently collide. Thus, a voltage 
gradient of approximately 20 volts per centimeter 
of arc length is sufficient to force the electrons 
along. Incidentally, it might be noticed that the 
product of this arc stream voltage gradient and 
the arc current give the wattage of the heat 
input to the arc stream necessary to maintain it 
in its ionized condition. 

Near the anode surface, the conditions are not 
so favorable for the conduction of current. The 
surface of the positive carbon is about 2000° 
cooler than the arc stream, and the gas immedi- 
ately in front of it consists largely of freshly 
evaporated carbon vapor in temperature equi- 
librium with the surface. At 4000°K, the carbon 
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vapor is a very poor conductor of electricity. 
Consequently, there are few positive ions in the 
narrow layer of gas close to the anode surface 
and it is therefore impossible to force a high 
density of electrons into this layer without 
developing a large unbalanced negative space 
charge. Since only a small number of electrons 
can exist close to the anode surface, in order to 
carry the arc current they must be forced through 
the layer at a very high speed. It therefore takes 
a high voltage drop concentrated across this 
narrow anode layer to force the electrons carrying 
the arc current across it and into the anode. In 
a pure carbon arc, this anode drop of voltage is 
about 35 volts. 

Most of the heat developed by this drop of 
voltage in the anode layer is transferred to the 
surface of the positive carbon, part of it by the 
impact of the highly accelerated electrons and 
part of it by thermal conduction through this 
narrow layer of gas. Conditions are favorable for 
a high heat transfer between the arc gas and the 
anode surface because of the extremely large 
temperature gradient between the two. When 
the electrons finally reach the anode surface, they 
provide an additional heating effect in giving up 
their heat of condensation. 

All of these factors tend to increase the tem- 
perature of the anode surface. Just how high 
its temperature becomes depends upon a balance 
between this heat input and the heat lost from 
the positive carbon by radiation, evaporation, 
and conduction. In general, this heat balance 
favors the highest possible temperature for the 
anode surface. The limiting temperature is 
the sublimation temperature of carbon, about 
3950°K. This temperature limits the maximum 
brightness of the low intensity arc to a value of 
about 180 candlepower per square millimeter. 
Most low intensity arcs operate with a maximum 
brightness slightly lower than this, from 160 to 
175 candlepower per square millimeter. An 
increase in current in the low intensity arc 
cannot, therefore, appreciably increase the 
maximum brightness but can only increase the 
area of the crater surface. 

It is interesting to observe why carbon is the 
ideal material for use as an electrode in such an 
arc. It remains solid at a higher temperature than 
any other material of suitable electrical and 
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thermal conductivity and its property of vola- 
tilizing directly from a solid to a gaseous state 
provides for a convenient disposal of the con- 
sumed portion. 

The flame arc was a natural development from 
the low intensity arc obtained by enlarging the 
core in the electrodes and replacing part of the 
carbon by chemical compounds capable of 
radiating efficiently in a highly heated gaseous 
form. These compounds are vaporized along with 
the carbon and diffuse throughout the arc 
stream, producing a richly colored arc flame, the 
color of the flame depending upon the compounds 
used. Typical materials of importance in provid- 
ing efficient radiation in various parts of the 
spectrum are iron for the ultraviolet, rare earths 
of the cerium group for white light, calcium 
compounds for yellow, and strontium compounds 
for red. 

Figure 2 is a picture of an alternating current 
flame arc. 

The chemical compounds added to the core 





Fic. 3. High intensity arc. 
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have a considerably lower ionization potential 
than carbon. The greater ease of ionization of 
these materials reduces the temperature of the 
anode layer necessary for the conduction of 
current into the anode, and therefore reduces the 
area and brilliance of the anode spot so that its 
contribution to the total light becomes unim- 
portant. The radiation emitted by the flame then 
consists chiefly of the characteristic line spectra 
of the elements in the flame material and the 
band spectra of the compounds formed. The 
concentration of flame elements in the arc 
stream is not very high; so that the flame arc, 
while brighter than many other light sources, is 
considerably less bright than either the low or 
the high intensity arc. Since the whole arc 
flame is made luminous, however, the light 
source is one of large area; and high radiating 
efficiencies of up to 80 lumens per watt are 
obtained. 

We have shown that the transition from the 
low intensity to the flame arc is accomplished by 
incorporating in the carbon electrodes elements 
or compounds to produce the particular band of 
radiation desired. In turn the high intensity arc 
may be obtained from the flame arc by increasing 
the current density in the electrodes to a point 
where the anode spot spreads over the entire tip 
of the carbon, resulting in a rapid evaporation 
of flame material as well as carbon from the core. 
The principal source of light is now the anode 
surface and the gaseous region immediately in 
front of it. Since the flame material is more 
easily ionized than the carbon, its presence in 
the anode layer results in a lower anode drop at 
the core area than at the shell of the carbon. 
This tends to concentrate the current at the core 
surface and results in the hollowing out of the 
crater as the current is increased. 

Figure 3 shows a typical high intensity are at 
a current of 125 amperes. 

The rapid evaporation of flame material 
produces a high concentration of this efficiently 
radiating gas in the crater and immediately in 
front of it. This gas, of course, radiates in all 
directions, even back toward the crater surface 
and consequently tends to serve as a blanket 
preventing the radiative cooling of the crater 
face. The heat dissipated at the crater face must 
then be entirely through evaporation of flame 
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material and carbon and through conduction 
back along the positive carbon. This, of course, 
tends to increase the evaporation of material 
within the crater and aids in the tendency for 
crater formation. Thus, in a high intensity arc 
there is a close correlation between the crater 
depth and the brilliancy of the arc gas within and 
immediately in front of the crater; for a given 
type of positive carbon there is a linear relation- 
ship between the crater depth and the excess 
brightness over that of a low intensity arc. 

An increase of current in a high intensity arc 
increases the crater area only slightly but pro- 
duces a marked increase in crater depth and 
brilliancy. The maximum brilliancy of the crater 
obtained in various types of direct current high 
intensity arcs used in common commerical lamps 
ranges from 350 to 1200 candlepower per square 
millimeter. Experimental carbons have been 
produced with brilliancies of the order of 2000 
candlepower per square millimeter. 

The problem of understanding the mechanism 
of light production in the high intensity arc is an 
intriguing one. The reason for this lies in the 
fact that a large proportion of the light emitted 
from the high intensity crater is continuous 
radiation and not, as one might expect, the line 
radiation of the elements. True, one would 
expect some continuous radiation from the incan- 
descent carbon at the base of the crater. How- 
ever, the carbon surface cannot exceed the 
sublimation temperature of carbon, and there- 
fore the brilliancy of this source of continuous 
radiation cannot exceed 180 candlepower per 
square millimeter. The observed brightness of 
the continuous spectrum is several times this. 
Furthermore, it appears to originate within the 
crater gas itself rather than at the bottom surface 
of the crater. Thus, we have the problem of 
explaining the emission of an extremely strong 
continuous radiation by a hot gas. 

We ordinarily assume that continuous radi- 


_ation is emitted by incandescent solids and that 


gases radiate the characteristic line spectra of 
the elements and the band spectra of compounds. 
This simple situation exists in the low intensity 
are and in the flame arc. However, in the high 
intensity arc these sources are no longer clearly 
segregated, since the gas emits such a strong con- 
tinuous spectrum. 
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It is true that many other cases of emission of 
continuous radiation by gases are known and 
some of these are understood. For example, there 
is the hydrogen discharge which gives off a 
strong ultraviolet continuous spectrum which is 
explained as a molecular continuum. Also, fairly 
strong continuous spectra have been obtained 
and recognized as electron recombination spectra. 
However, the continuous spectrum appearing in 
the high intensity arc is not recognizable as 
belonging to any of these standard types. 

Several explanations have been put forward 
in the literature for the strong continuous radi- 
ation from the high intensity arc. P. R. Bassett! 
in 1921 gave a theory that solid incandescent 
particles of cerium carbide are formed in the 
crater of a high intensity arc and that the con- 
tinuous radiation is obtained from them. There 
are a number of indications that discredit this 
theory. Perhaps the most important is that a 
strong continuous spectrum is obtained using a 
wide variety of coring materials, not all of which 
could form extremely high boiling compounds in 
the arc crater. 

Another explanation was introduced by Fin- 
kelnburg? in 1935. He thought that the con- 
tinuous spectrum was an electron radiation. In 
the arc discharge, the fast moving electrons 
frequently collide with positive ions. When an 
electron comes close to a positive ion, it is 
attracted toward it because of its opposite 
charge. Sometimes the electron will be captured 
by the positive ion to form a single neutral atom 
and sometimes it will be only deflected and 
retarded in its progress. In both cases, the elec- 
tron loses a certain amount of energy due to its 
loss in velocity. The lost energy is radiated. 

When the electron is captured by the positive 
ion, it not only loses its kinetic energy of motion, 
but it also loses the énergy that it takes to get 
the electron out of such an atom; that is, its 
ionization energy. The ionization energy can 
take certain discrete values, depending upon the 
energy state of the atom into which the electron 
falls. The energy radiated in the case of capture 
of the electron is the total energy lost by it, or 
the sum of the ionization and kinetic energies. 
Since the wave-length is inversely proportional 


1P. R. Bassett, Trans. A. E. S. 44, 153 (1923). 
2 W. Finkelnburg, Zeits. f. Physik 93, 201 (1935). 
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to the energy radiated, and since the energy 
radiated may take on any value greater than the 
ionization energy, the wave-lengths radiated 
may be of any wave-length shorter than the 
particular wave-length corresponding to the 
ionization energy. The spectrum emitted by the 
recombination of an electron with a positive ion 
to form a neutral atom is therefore continuous 
in small bands, each band being associated with 


important. Furthermore, it is sure to exist in 
some measure in the arc, since there are always 
electrons and positive ions present. However, 
there are also things to be said against explaining 
the major part of the observed continuous 
spectrum by the electron radiation theory. No 
evidence is found of any bands such as would 
occur with electron recombination. Further, the 
wrong kind of atoms give the best continuous 
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the energy state of the atom into which the 
electron may fall. Experimentally it is found that 
these electron recombination bands are rather 
wide and diffuse without any sharp boundary. 

When an electron collides with a positive ion, 
but is not captured, it loses a certain amount 
of energy by radiation corresponding to its loss 
in velocity during the collision. The energy lost 
in the collision is the difference in kinetic energy 
before and after the collision. Since this dif- 
ference can vary continuously, any amount of 
energy can be radiated, and the radiation will be 
continuous. Finkelnburg gave this type of radi- 
ation the rather catchy title ‘‘free-free’’ radiation, 
because the electron is free from the influence of 
any other particles before and after its encounter 
with the positive ion. While free-free radiation 
is intriguing as a source of light, calculations 
show that it is probably negligible in comparison 
with the recombination continuous spectrum 
emitted by electrons. 

The theory of electron radiation has some 
factors in its favor. The principal factor is a 
negative one: since not too much is known about 
electron radiation experimentally, a chance 
always remains that it might turn out to be 
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spectrum. One would expect the strongest con- 
tinuous radiation to come from arcs cored with 
alkali metals if the continuous spectrum were due 
to electron radiation, because these atoms are 
most easily ionized and produce the most 
copious supply of electrons and positive ions. It 
is found, however, that these alkali metals are 
among the poorest light emitters. 

An examination of the elements which produce 
the strongest continuous radiation in a high 
intensity arc provides the key to what is probably 
the explanation of its origin. We find that the 
continuous radiation is strongest with the flame 
constituents that provide the largest number of 
strong spectrum lines, and that the continuous 
spectrum is most intense in the wave-length 
regions occupied by the greatest density of these 
lines. Therefore it seems probable that the source 
of continuous radiation is the broadening of 
spectrum lines by their transmission through the 
are vapor. 

Spectrum lines are ordinarily thought of as 
being sharp with very little breadth, because the 
ordinary spectroscope or spectrograph is incap- 
able of resolving them. However, it is well known 
that spectrum lines do have some finite breadth. 
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produced by several factors. In the first place 
they have a natural breadth because the atoms 
exist in their excited states for only a limited 
length of time. This natural breadth of lines is 
negligibly small in the are source. Another factor 
broadening lines is the Doppler effect whereby 
there are slight changes in the frequency of the 
light, depending upon the speed at which the 
emitting atom is traveling towards or away from 
the observer. The collisions of the atoms during 
their radiating processes also produce slight 
changes in the frequency of the light emitted 
and tend to broaden the spectrum lines. These 
latter two factors give a half-breadth of line of 
about 0.03A as initially emitted in the arc gas. 
Thus the lines as initially emitted are still quite 
sharp. 

However, these lines are broadened further by 
the process of absorption and re-emission as they 
are transmitted through the crater gas. The 
dense group of flame atoms in the hot crater of 
a high intensity arc quickly build the density of 
radiation in the line up to the point where the 
center of the line is as strong as it can be. This 
point is reached when the radiation density is so 
great that the atoms absorb energy as fast as they 
emit it. The maximum radiation intensity under 
these conditions is the radiation intensity of a 
blackbody at the same temperature as the arc 
gas. The outskirts of the line have the same limit- 
ing intensity as the center, and if the radiating 
depth of the gas is great enough, this slightly 
off-key radiation will also build up to a blackbody 
intensity. Even at some distance from the center 
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Fic. 5. Spectrum of mag- 
nesium resonance lines broad- 
ened in a high intensity arc 
crater. The three strongest 
lines are Mg II lines at 2795.5 
and 2802.7A and a Mg I line 
at 2852.1A. An image of the 
positive crater face of the arc 
was focused on the slit as 
shown in the drawing on the 
left. The spectrogram of any 
part of the arc lies at the same 


horizontal level as that at 
which it is shown in the draw- 
ing. 


of the line, an appreciable intensity of radiation 
is built up. The calculated intensity distribution 
of such a broadened line is shown in Fig. 4. The 
assumptions used for the calculation of this 
curve were chosen to correspond to conditions in 
the high intensity arc crater, and to a flame 
element such as cerium that has a multitude of 
lines of low excitation potential. It can be seen 
that the intensity of radiation does not fall off 
rapidly even at a relatively long wave-length 
interval from the center of the line. 

The broadening in the case of an element with 
strong resonance lines is shown in the spectrum 
reproduced in Fig. 5. This is the spectrum of a 
magnesium cored high intensity arc, and shows 
the great broadening of the resonance lines within 
the are crater. With flame elements producing a 
large number of weak but closely spaced lines, 
the broadening of the individual lines is much 
less extensive, but the overlapping tails of these 
broadened lines combine to give an intense con- 
tinuous spectrum. 

The sources of radiation of the carbon arcs, 
then, are the continuous spectrum of an incan- 
descent solid in the low intensity arc, the line 
and band spectra of elements and compounds in 
the flame arc, and a combination of these with a 
continuous spectrum originating from the trans- 
mission broadening of spectrum lines in the 
high intensity arc. All of these types of arcs have 
been under continuous development since their 
discovery, and at present they have reached a tech- 
nical excellence that gives them a preeminent place 
wherever an intense source of radiation is needed. 
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The Mechanism of the Electric Spark 


By Lreonarp B. LOEB AND JOHN M. MEEK. Pp. 
188+-xiii, Figs. 43, 16233 cm. Stanford University 
Press, Stanford University, 1941. Price $3.50. 


While disruptive electrical discharges through gases have 
long been the subject of intensive study, still, owing to the 
extremely short time involved in spark and lightning dis- 
charges, only in recent years has reliable experimental 
evidence been obtained regarding the important initial 
stages of these discharges. This evidence has disclosed new 
significant features which now serve as guides to a better 
understanding of the somewhat involved processes opera- 
tive in these discharges. Another advance rests on the 
general acceptance of the hypothesis that in the initiation 
of sparks in gases at atmospheric pressure a prominent part 
is played by photo-ionization of molecules of the gas by 
photons, liberated in the gas as a result of electron impact. 
Based on the vantage ground of this new knowledge, there 
has now been rounded out a comprehensive theory of the 
whole mechanism of spark discharges, which is described in 
the volume under review, by authors who themselves have 
played an important role in its development. 

In the opening chapter, a brief account is given of the 
Townsend theory of the spark. While this theory accounts 
adequately for spark formation when the product of the gas 
pressure by the spark distance is below a certain value, 
several reasons are given why the theory cannot be valid 
when larger values of the above product prevail. Among 
these objections perhaps the most crucial is the one that 
according to this theory the time required for the develop- 
ment of a spark is about a hundred times larger than that 
actually found by experiment. 

The second and main chapter deals with the new streamer 
theory of spark discharge. Briefly stated, the principal 
feature of this theory, as applied to short sparks, is that a 
Townsend electron avalanche as it approaches the anode 
may under certain conditions become converted into a 
positive streamer moving in the opposite direction. An 
electron avalanche leaves a trail of positive ions in its wake 
which rapidly increases in density as the avalanche pro- 
ceeds on its way. The space charge arising from these 
comparatively immobile ions produces a radial field about 
the avalanche, reduces the applied field ahead of it, and 
enhances that behind it. This rearward field may become so 
strong that any electrons, produced in it by photons or 
drawn into it from the outside, will cause intense ionization 
there. The new supply of positive ions rapidly extends the 
intense field backward toward the cathode, and when the 
streamer reaches that terminal a large flow of electrons 
moves to the anode through the now plasma-like con- 
ducting channel, and this flow constitutes the main dis- 
charge current. 

The next step was to determine quantitatively the con- 
ditions under which these positive streamers can form. The 
stand is taken that the crucial condition to be satisfied is 
that the radial field of the avalanche be strong enough to be 
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able to draw into the head of the avalanche a large number 
of the photoelectrons from the surrounding space. Meek 
made the arbitrary assumption that to do this effectively 
the radial field X, at the outer edge of the avalanche must 
be equal to the impressed field X. So an equation for X,, 
involving X, is derived on the bases of certain simplified 
conditions and when this value is set equal to X, this 
quantity becomes the field required for the production of a 
spark. The evaluation of X requires in the end the making 
use, by a trial and error method, of values obtained from a 
curve giving the experimental relation of a/p as a f(X/p). 
Notwithstanding the approximations made and the inaccu- 
racy of some of the data used, the values of X thus obtained 
are found to agree uncannily with experimental values of 
the fields required to produce sparks. Misgiving may arise, 
at first, regarding this close agreement when the possible 
errors introduced by some of the assumptions are con- 
sidered. Thus the number of positive ions contributing to 
X1 is taken to be the number produced in a uniform field X 
by one electron in its motion from the vicinity of the 
cathode. The field ahead of the avalanche is however 
admittedly decreased by the positive charge of the ava- 
lanche and so the all important ionization coefficient a 
must also diminish as the avalanche progresses. No positive 
ions are supposed to be added to the avalanche by the 
action of the other electrons drawn into the channel by X, 
itself. The ions spread, as they are, along the avalanche are 
supposed to occupy a sphere of radius equal to the distance 
electrons diffuse, on the average, in the time the avalanche 
has been advancing. Evidently a substantial number of 
them is at a greater distance from the axis of the avalanche. 
Perhaps the most serious difficulty of all arises from the fact 
that the data available for the relation of the coefficient of 
diffusion to the mobility of electrons, which is involved in 
the equation for X, may well be in error by a factor of two 
or three or more. 

However in the end it turns out that somehow errors of 
this magnitude do not much matter because it is found that 
the criterion X,=X is not a rigid one and that, owing to the 
presence of an exponential term in the formula, almost the 
same value is obtained for X when, instead of the above, 
the relation used is X¥,=0.1X. 

The third and final chapter of the book is devoted to a 
calculation of breakdown potentials in various types of 
gaps. The book considers not only short sparks where a 
uniform applied field obtains but also long sparks where the 
fields are not uniform, lightning discharges in their various 
forms, sphere gaps, coaxial cylinders, and corona discharges. 
Attempts are made to explain all of these by some form of 
the streamer theory. If the explanations at times seem 
somewhat forced, it must be remembered that the phe- 
nomena under discussion are very complex and that in all 
likelihood the theory has not as yet reached its final form in 
all details. 

Those interested in disruptive discharges will be glad to 
have available this useful compilation of recent papers on 
this latest theory of spark formation. 

JOHN ZELENY 
Yale University 
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National Lectureships of Sigma Xi 


National lectureships on scientific subjects for 1942 have 
been announced by Sigma Xi, Science Service reports. 
Lecturers and their topics will be: 


Energy Production in Stars, H. A. Berne, Department of Physics, 
Cornell University 

Some Recent Work in the Field of High Pressures, P. W. BRIDGMAN, 
Department of Physics, Harvard University 

Structure of Liquids, Joun G. KirKwoop, Department of Chemistry, 
Cornell University 

Modern Power Generation, Lionet S. MARKS, Department of Engi- 
neering, Harvard University 

Recent Results from Studies on the Anterior Pituitary, H. M. Evans, 
Department of Biology, University of California. 


* 


Necrology 


Rollin Landis Charles, Professor and Head of the De- 
partment of Physics at Franklin and Marshall College, died 
on December 13 at the age of 56. 

Walter Mathew Dunagan, Associate Professor of Theo- 
retical and Applied Mechanics at Iowa State College, died 
on November 24 at the age of 47. 


* 


Bulletin of Mathematical Biophysics 


The March, 1942, issue of the Bulletin of Mathematical 
Biophysics, edited by N. Rashevsky and published by the 
University of Chicago Press, will contain the following 
articles: 

Diffusion as a Function of Aggregation in Colloidal Media, HERMAN 
BRANSON 

Theory of Steady-State Activity in Nerve-Fiber Networks: IV. N Cir- 
cuits with a Common Synapse, ALSTON S. HOUSEHOLDER 

Kinetic Theory of Diffusion Forces in Metabolizing Systems, H. D. 
LANDAHL 

Suggestions for a Mathematical Biophysics of Auditory Perception 
with Special Reference to the Theory of Aesthetic Ratings of Combina- 
tions of Musical Tones, N. RASHEVSKY 

Non-Linear Excitation Theory: Non-Accommodative, Sub-Threshold 
Effects, ALVIN M. WEINBERG. 


* 


Recent Awards 


The John Fritz Medal for 1942, sponsored by the four 
leading national engineering societies, will be awarded to 
Dr. Everette Lee DeGolyer, consulting petroleum engineer 
of Dallas, Texas, in recognition of his work in the application 
of geophysical exploration to the search for oil fields, 
according to a recent issue of Science. The award will be 
presented at a dinner of the American Institute of Mining 
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and Metallurgical Engineers at the Waldorf-Astoria Hotel, 
New York City, on January 14. 


Dr. William David Coolidge, Vice President of the 
General Electric Company in Charge of Research and 
Director of the GE Research Laboratory at Schenectady, 
has been awarded the Duddell Medal of the Physical 
Society in London. This medal is given ‘“‘to persons who 
have contributed to the advancement of knowledge by the 
invention or design of scientific instruments or by the dis- 
covery of materials used in their construction.”’ It has been 
awarded seventeen times previously. Two other Americans, 
the late Professor A. A. Michelson of the University of 
Chicago and Professor W. G. Cady of Wesleyan University, 
have received it. In making the award the Society cites Dr. 
Coolidge’s achievement in producing ductile tungsten from 
the powdered metal and the development of the heated- 
filament cathode x-ray tube which is known by his name. 


The Edison Medal for 1941 has been awarded by the 
American Institute of Electrical Engineers to Dr. John 
Boswell Whitehead, Director of the School of Engineering 
of the Johns Hopkins University, ‘‘for his contributions to 
the field of electrical engineering, his pioneering and 
development in the field of dielectric research, and his 
achievements in the advancement of engineering educa- 
tion,” a recent issue of Science states. 


* 


New Defense Training Courses 


The Graduate Electrical Engineering Department of the 
Polytechnic Institute of Brooklyn is offering, under the 
Engineering Science and Management Defense Training 
Program of the United States Office of Education, a group 
of four courses and one seminar in ultra-short waves for 
graduate students of physics or electrical engineering, as 
follows: 


Ultra-Short Electromagnetic Waves, PRoressorR E. WEBER 

Electrical Measurements at Ultra-High Frequencies (300 to 3000 
Megacycles per Second), Proressor G. B. HOADLEY 

Fundamentals of Television Engineering 

Ultra-High Frequency Electronics, Dr. D. O. Nortu, RCA Radio- 
tron, Harrison, New Jersey 

Seminar on Selected Topics, Proressors G. B. HOADLEY, E. WEBER, 
Dr. D. O. NoRTH, AND OTHERS. 


The program is to start the week of January 26. Prefer- 
ence is given to applicants engaged in important defense 
projects. 


* 


American Society for X-Ray and Electron Diffraction 


An invitation extended to the newly-formed American 
Society for X-Ray and Electron Diffraction to become an 
Associated Society of the American Institute of Physics 
was accepted by that organization at a business meeting 
held December 31, 1941, at Cambridge. 
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Anisotropy of Ferromagnetic Powder Particles 


E.tiot T. BENEDIKT 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received October 3, 1941) 


Measurements of the magnetic moment of sediments of ferromagnetic particles deposited 
in a magnetic field show that the particles are magnetically anisotropic. The dependence of 
this anisotropy on the magnetic field is investigated for sediments prepared under different 


conditions, and its probable origin discussed. 


INTRODUCTION 


HE present work deals with some magnetic 
properties of sediments deposited in a 
magnetic field from suspensions in water of fer- 
romagnetic powders. The experiments were 
performed on powdered iron (reduced by 
hydrogen), ferromagnetic FexO; and Fe;QO,. The 
size of the particles was determined micro- 
scopically, and by Stokes’ law. The iron particles 
were of an irregular cubical shape, and practically 
all of them had a side of about 8X10~ cm; the 
oxide particles, of an irregular spherical shape, 
had a diameter of the order of 10-* cm, though 
some of them were considerably smaller. (FesO; 
particles having a diameter of ~4X10~ were 
observed. ) 
EXPERIMENTAL 


1. Preparation of the Specimens 


For the preparation of the sediments, the fol- 
lowing apparatus was used (see Fig. 1). A glass 
cylinder ¢c was placed inside a pair of Helmholtz 
coils hh, which generated a field of 20.8 oersteds 
per amp., and whose axis was in the direction of 
the horizontal component of the terrestrial 
field, found to be 0.275 oersted in our laboratory. 
Near the bottom of the cylinder, with its center 
coincident with the center of the system of coils 
a glass disk d, ?’”’ in diameter and of known mass 
with a straight line marking a direction on its 
surface was resting horizontally on the brass 
support ss, and kept in place by two brass springs 
s’s’. The support could be rotated around its 
vertical axis so that the direction marked on the 
surface of the disk formed with the axis of the 
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coils a given angle # measured on the goni- 
ometer g. 

The sediments were prepared by fixing the 
direction of the disk, applying a given field in 
the coils, and introducing from the top a con- 
centrated suspension of the ferromagnetic powder 
into the cylinder, already filled with water. 
After a short time (5 min. for iron, 15 min. for 
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Fic. 1. Settling chamber. 


the oxide particles) practically all the particles 
were deposited. The current in the coils was 
then interrupted and the liquid very slowly 
evacuated through the bottom of the cylinder. 
After about 3 hours the main part of the water 
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mixed with the deposited powder having been 
evaporated, the glass disk, covered with a thin 
layer of sufficiently compact sediment, was 
removed from the support and placed in a des- 
iccator where it was completely dried in about 
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Fic. 2. Torque in dyne cm acting on a sediment of iron 
powder particles deposited in a field of 10.7 oersteds as 
function of the magnetic field in oersteds, acting at 45° 
with the direction of sedimentation. 


12 hours. This procedure prevented oxidation in 
the case of iron particles. The mass of the sedi- 
ment was later determined by difference of 
weight; it was generally of the order of 4-5 
centigrams for the iron, 2-3 decigrams for the 
oxide sediments. 


2. Measurement 


The measurements of the magnetic moment 
of these specimens was accomplished by means 
of a magnetometric device, consisting of another 
pair of Helmholtz coils giving a field of 24.7 
oersteds per amp. A light glass support, sus- 


. pended on a phosphor bronze strip whose torsion 


constant was 3.18 dyne cm/rad. carried a glass 
disk exactly similar to the one used for the 
specimen with its plane horizontal, its center 
coincident with the center of the system of coils, 
and the direction marked on it coincident with 
the axis of the coils when the suspended system 
was at rest. A concave spherical mirror was 
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attached to the upper part of the support and 
the deflection of the suspension measured, as 
usual, by means of a beam of light reflected on 
a plane scale 1 m distant from it. The disk carry- 
ing the sediment was placed on the suspension 
with the direction marked on it coincident with 
the one marked on the support. Thus, the field 
of the magnetometer formed an angle equal to 
3% with the direction in which the sedimentation 
field was acting with respect to the specimen. 
The deflections were then measured. These were 
obtained by increasing the field 7 of the mag- 
netometer at regular intervals of 2.47 oersteds 
(corresponding to intervals of 0.1 amp. of the 
current flowing through the coils) from 0 to a 
maximum, generally 29.7 oersteds, then decreas- 
ing through 0 to an equal and opposite value, and 
increasing again to the original maximum value. 
From the deflection the torque acting on the unit 
mass of ferromagnetic material was then calcu- 
lated from the formula Q=(k/200m)6 dyne cm 
(where k is the specific torsion of the suspension, 
m the mass of the sediment, 6 the deflection” in 
cm read on the scale, and corrected for cur- 
vature). Figure 2 represents the function 
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Fic. 3. Torques in dyne cm acting on sediments of iron 
powder particles deposited in a field of 10.7 oersteds (con- 
tinuous lines) and in the magnetic field of the earth, 0.275 
oersted (dotted line), as function of the magnetic field in 
oersteds acting at different angles with the direction of 
magnetization. 
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Q=(H) obtained with a specimen deposited in 
a field of 10.7 oersteds for 8=45°. The dotted 
line represents the ‘‘virgin’’ curve obtained when 
increasing the field from 0 to the maximum for 
the first time, the superior continuous curve 
when decreasing it, the inferior when increasing 
it again. By repeating the cycle the last two 
curves are obtained again. The function Q(/) 
was of this type in all the cases. Other results are 
shown in Fig. 3. The continuous curves represent 
the part of the Q(/7) function obtained decreasing 
the field from its maximum to its minimum on 
specimens deposited in a field of 10.7 oersteds 
with different values of 3. For #=90° the above 
function becomes linear. The dotted line repre- 
sents measurements obtained with a sediment 
deposited under the action of the field of the 
earth alone, for }=45°. The same type of results 











Fic. 4. Anisotropy function of a sediment of iron powder 
particles deposited in a field of 5.27 oersteds. 


were obtained with sediments of ferromagnetic 
FesO3; and Fe;QO,. The above results can be 
reproduced only within about 10 percent, 
because of inevitable displacements after the 
sedimentation of the particles at the surface of 
the specimen, due to unavoidable little mechan- 
ical shocks, flow of the liquid during its evacu- 
ation, and the magnetic torque acting on them 
in the course of the measurements. 


THEORY 


The experimental results can be explained by 
assuming the sediments to be magnetically aniso- 
tropic. Their magnetic moment per unit mass 
can be written in the form: 


Y= 2X0t+214+ =X, 
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Fic. 5. Anisotropy function of a sediment of iron powder 
particles deposited in a field of 16.1 oersteds. 


where Xp is the constant remanent moment with 
which the specimen is left by the sedimentation 
field h, in a direction paralell to it, £1, Z_2 the 
moments induced by the applied (magneto- 
metric) field H' in directions respectively 
parallel and perpendicular to Xo. These directions 
are also assumed to be the principal directions 
of magnetization of the specimen. 

It is logical to consider 2; and Y2 functions of 
the projections H,=H cos 3’, H,.=H sind of H 
along their directions: 


ys 
a 


1=2,(M), Y2= 22(H2). (1) 


The component of the magnetization normal to 


1The demagnetizing factor was found, assuming the 
sediment to have a thickness of the order of a few tenths 
of a millimeter and considering it an oblate spheroid to be 
at most of the order of 0.2 [see, for instance, L. P. Tarasov, 
‘“‘Demagnetization coefficient of oblate spheroids,’’ Phys. 
Rev. 53, 105 (1938) ]. The susceptibility (per unit volume) 
of the powders for the fields used was found by means of 
ballistic measurements performed on tubes containing the 
substance to be about 1.5 for iron, 0.06 for oxide powders. 
The demagnetizing field was therefore for iron sediments 
at most 15 percent, for oxide sediments (which were some- 
what thicker) 1.5 percent of the applied field. It was 
neglected, since difficult to establish, the main difficulty 
lying in the determination of the thickness of the specimen. 
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H will be: 
y= Zo sin +>; sin )—L.2 cos ¥ (2) 


and the torque acting on the unit mass of the 
specimen : 


Q=,/T =o sin 3- H+ 3(xi1—x2) sin 2077, 


where x:=2:/H, x2=22/H. The latter equation 
accounts for the parabola-like effect of the Q(/7) 
curves, since x1, x2 vary slowly with /7. 

If #=45°, we get from (1), (2): 


V22n = LotX1(Ho) —22(Mo), 
where //)=/1,;=H2.=(1/v2)H. The function: 
A(Ho) = 2,(1To) — Y2(Ho) 


may be deduced from the experiments, and is a 
measure of the anisotropy of the sediment; we 
shall call it anisotropy function. Figures 4 and 5 
show this function obtained with sediments of 
iron particles deposited in different fields. Figure 
6 shows the function A obtained with sediments 
of ferromagnetic Fe2O3, and Fe;O,4. It exhibits 
the same features of the one obtained with iron 
powder, but is considerably smaller. The dotted 
lines represent the “virgin” function and are 
obtained when increasing the field from 0 to its 
maximum for the first time, the upper and lower 
full lines, respectively, when decreasing and in- 
creasing it. These curves have the shape of 
unsymmetrical magnetization curves due to the 
fact that they represent the difference of two 
“‘displaced’’ magnetization functions, i.e., ob- 
tained from already magnetized material (such 
functions are unsymmetrical). 

The anisotropy of the sediment must be 
obviously due to the anisotropy of the particles, 
or group of particles by which it is formed. Let 
@, be the component of the magnetization in the 
direction of highest susceptibility, 2 in a direc- 
tion perpendicular to it.* When settling in a 
magnetic field h the particles will experience a 
torque rotating @, toward its direction. The 
energy required to turn o; in a direction per- 





? All the results are expressed in absolute electromagnetic 
units. 

* This does not imply that the directions in a plane 
perpendicular to @; are magnetically equivalent. 
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is the mass of the particles, (which was calculated 
from their size); ¢:—o2 is at least equal to A.‘ 

Thus, for iron particles € is at least of the order 
of magnitude of 10~-'°; for oxide particles of the 
order of 10~°, and only exceptionally as low as 
10-". Therefore, since ¢ is much greater than the 
energy of the thermal agitation kT =4-10-" erg 
at room temperature, all the particles are com- 
pletely oriented in the direction of h.5 After the 
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Fic. 6. Anisotropy function of sediments of ferromag- 
netic FesO; (curve C) and Fe;O, (curve C’) powder parti- 
cles, deposited in a field of 21.0 oersteds. 


sedimentation, they will have a remanent mo- 
ment @o in the direction of highest magnetization. 
From the preceding discussion we see that we 
must have: 


Lo=Go, Li(Ho)=G1(o), Y2(FTo) =62* (Ho), 


where the bars denote average values, and o2* 
is the projection of #2 on the plane of the speci- 
men. The function A, therefore, is also a measure 
of the anisotropy of the particles. 

The magnetic anisotropy of the particles is 





‘ The anisotropy of the sediment can be only lower than 
that of the particles, if we assume only a partial orientation 
of these. 

5 The magnetic field due to the other particles is neg- 
lected; for its order of magnitude see note (1). 
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probably due to the fact that they are formed of 
a small number of anisotropic single crystals, and 
partly to an anisotropic demagnetization deter- 
mined by their shape. The last effect is, however, 
almost completely eliminated due to the fact 
that the particles of the sediment form a con- 
tinuous layer, and consequently the demag- 
netizing actions, which are due to the free 


magnetism appearing on the boundaries of each 
particle, are almost canceled by the effect of the 
free magnetization appearing on the neighboring 
particles. 

I want to express here my deepest gratitude 
to Professor F. Bitter of the Massachusetts 
Institute of Technology, who made the present 
work possible. 





Indexing Powder Diffraction Patterns of Isomorphous Substances 


Lupo K. FREVEL 
The Dow Chemical Company, Midland, Michigan 
(Received October 29, 1941) 


Two methods of indexing powder diffraction patterns of isomorphous substances are pre- 
sented: the first method depends on matching the logd values of two or more isomorphs 
(d being the interplanar spacing of a given reflection); the second method relates the spacing 
shifts Ad/d of two isomorphs to the symmetry of the particular crystal system. These methods 
are made practical by the publication of tabulated powder diffraction data. 


HE problem of indexing the powder photo- 

graphs of two or more isomorphous sub- 
stances resolves itself into two cases; namely, 
the case for which the lattice constants of one of 
the isomorphs are known, and the other for 
which no lattice parameters have been deter- 
mined from single crystal data for any member 
of the isomorphous group. 


CASE I. LATTICE CONSTANTS KNOWN FOR ONE 
MEMBER OF ISOMORPHOUS GROUP 


Since, in general, an approximate proportion- 
ality exists between corresponding lattice con- 
stants of two isomorphous substances, it becomes 
expedient to compare their log d values, where d 
is the interplanar spacing of a given reflection. 
The simplest case to be considered is isomorphism 
in the cubic system. Here the log d values of 
corresponding reflections of any two cubic 
isomorphs! differ by an additive constant. 
Accordingly, one log d strip will suffice for each 
cubic isomorphous series. In view of the fact 
that the (hkO) reflections in the tetragonal 
system or the hexagonal system involve only 


‘The log d values of corresponding reflections of any 
two cubic substances will differ by a constant. 


VOLUME 13, FEBRUARY, 1942 


one parameter, it is readily seen that the argu- 
ments for the cubic system apply essentially to 
the (hkO) reflections of isomorphs in these two 











TABLE I 
 — 

SnO2 (rutile) gq 444 
(hkl) d d 2 Set dp? : dg? : dr-?: 
110 3.34A 3.24A 0.030 a (least) 2:4:5:8:10 
101 2.64 2.49 0.059 
200 2.36 2.29 0.030 a 
111 2.30 2.19 0.049 
120 2.11 2.05 0.029 a 
121 1.75 1.69 0.035 
220 1.67 1.62 0.030 a 
002 1.58 1.485 0.062 b (largest) 
130 1.492 1.449 0.029 a 


Filtered MoKa radiation was used to obtain the powder 

diffraction patterns of stannic oxide and rutile. (ki) are 

the indices of the various reflections. d=interplanar 

: : Ad ’ 

spacing measured in A. — the mean fractional change 
d+ D 

of corresponding interplanar spacings. The ratio test 

(d,* : d,* : d,* :) identifies the tetragonal system. 


systems. Figure 1 gives a graphical illustration 
of this matching of log d;xo for two or more 
isomorphs in the tetragonal system. 

In the case of the orthorhombic system, the 
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Fic. 1. The respective log d values of the four isomorphous groups are matched by appro- 
priate translations along the abscissa. For the cubic system (In2O; and Y.O;), all reflections 
with the same indices (kl) match completely. The matching of (hk0) reflections in the tetrag- 
onal system is indicated by dotted lines. In the case of the orthorhombic system (BaSQ,, etc.) 
the matching of corresponding reflections is only approximate. 


matching of corresponding reflections on a 
sliding logarithmic scale is facilitated either by 
selecting as a comparison standard such an 
isomorph that the more prominent reflections 
have about the same relative intensities for the 
two isomorphous substances or by approxi- 


110 


mating the axial ratio of the un-indexed sub- 
stance by forming solid solutions of the standard. 

The utility of the comparison method of 
indexing isomorphous substances depends on the 
availability of correctly indexed powder diffrac- 
tion data for the various known crystal struc- 
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_Ad Ad Ad 

Hgols Hg2Br2 Hg2Cle Ad Ad Ad 
(hkl) a "7 - Set (a 2 ), Br (7:3), Cl (733), Cl dp? : dqg7? : dre7?: 
101 4.544 4.30A 4.16A 0.054 0.087 0.033 
110 3.49 3.29 3.17 a (largest) 0.059 0.096 0.037 1:2:4 
004 2.90 2.77 ye b (least) 0.046 0.064 0.018 
200 2.46 2.32 2.24 a 0.059 0.094 0.035 
114 2.23 2.12 2.06 0.050 0.079 0.029 
105 2.10 2.01 1.97 (0.044) 0.064 0.020 
204 1.88 1.78 1.73 0.055 0.083 0.028 
220 1.74 1.64 1.58 a 0.059 0.096 0.037 
215 1.60 1.52 1.478 0.051 0.079 0.027 


Filtered MoKe radiation was used to obtain the powder diffraction patterns of the mercurous halides. (hkl) are the 





indices of the various reflections. d=interplanar spacing measured in A. Ad =the mean fractional change of corre- 
d+ 

sponding interplanar spacings. The ratio test (d,?:d,?:d,*:) identifies the tetragonal system. (See International 

Tables for the Determination of Crystal Structures (1935), Vol. II, 455-457, 460-464.) Of the two choices of indices for 


the prism reflections; namely, n(100, 110, 200) or n(110, 200, 220), the latter set for »=1 permits of satisfactory indexing 
of the remaining reflections; the second set for m=2 leads to a larger but equivalent lattice. 
































rasce III. 
( Ad ) Ad Ad Ad ‘ Ad Ad 

BaSeO. BaCrOs BaMnO« BaSOx ad ( Ad ( 3) ( ai) ( Ad ( ad 
(hkl) d d d d Set c+ 2 / Se, Cr é+ 2 / Se, Mn é+ 2/Se,S d+ 2/Cr, Mn é+ 2/Ce.$8 é+ 2/Mn,S 
200 — A 4.550A — A 4.435A —_ — —_— — 0.0256 _— 
oll 4.430 4.396 4.323 — —_ — 0.0077 0.0244 0.0167 
111 4.031 3.983 3.955 3.885 0.0120 +0.0035 0.0190+0.0035 0.0369+0.0035 0.0071+0.0035 0.0249+0.0035 0.0179+0.0035 
201 — - - 3.779 —- —_ = oa —- — 
002 3.680 3.668 3.662 3.573 a 0.0033 0.0049 0.0295 0.0016 (least) 0.0262 0.0246 (largest) 
210 3.535 3.527 3.503 3.436 0.0023 0.0091 0.0284 0.0068 0.0261 0.0193 
102 3.405 3.403 3.390 3.309 * 0.0006* 0.0044* 0.0286 0.0038 0.0280* 0.0242* 
211 3.187 3.174 3.157 3.100 0.0041 0.0095 0.0276 0.0054 0.0235 0.0182 
112 2.923 2.898 2.883 2.830 0.0086 0.0138 0.0323 0.0052 0.0237 0.0186 
020 2.846 2.769 2.739 2.720 b 0.0274 (largest) 0.0383 (largest) 0.0453 (largest) 0.0109 (largest) 0.0179 (least) 0.0070 (least) 
212 2.549 2.535 2.526 2.480 0.0055 0.0091 0.0274 0.0036 0.0219 0.0182 
220 2.366 2.344 2.317 —_— — _— 0.0093 0.0209 0.0116- 
221 2.286 2.256 2.234 2.205 0.0132 0.0230 0.0361 0.0098 0.0229 0.0131 

BaSeOs BaCrO.s BaMnO,s BaSO.« Filtered MoKea radiation was used to obtain the powder diffraction 
a 9017+40.008A 9.125+0.030A 9.0752+0.025A 8.885+0.022A patterns of the above isomorphous barium salts. Camera radius =20.73 
b 5.692 +0.006 5.535 0.010 5.483 +0.013 5.433 +0.009 cm. (hkl) are the indices of the various reflections. d =interplanar spacing 
ce 7.357+0.005 7.333 +0.010 7.312 +0.012 7.141 +0.010 Ad 


measured in A. =the mean fractional change of corresponding 
d 
+ 2 
interplanar spacings. 





tures. The published tabulated powder data?‘ CASE II. LATTICE CONSTANTS UNKNOWN 
are quite useful in establishing isomorphism; 
however, a refinement of the data is needed to 
permit the unambiguous indexing of the ten 
innermost reflections of non-cubic substances. 
In general it is advisable to compare more than 
two isomorphs inasmuch as a series of com- 
parisons frequently reveals the ambiguity of a 
reflection. 


The first method of indexing is simple and 
direct and is convenient to use when suitable 
comparison standards are available. However, 
it is also possible to utilize isomorphic com- 
parisons to aid materially in the indexing of 
non-cubic powder patterns for which only 
isomorphism has been established. In principle 
sesakocce the method is the same as that of relating the 

2. I. Mikheiev, V. N. Dubinina, and G. A. Kovalev, temperature shifts> of a powder pattern to the 
Annales de I’Institut des Mines 4 Léningrad, 11, IIT (1938). anisotropy of a crystalline substance. The 


3]. D. Hanawalt, H. W. Rinn, and L. K. Frevel, J. I. : 
E. é 10, 457 (1938). aie temperature shifts are merely replaced by the 


_* A, K. Boldyrev, V. I. Mikheiev, G. A. Kovalev, and 
V. N. Dubinina, Annales de I’Institut des Mines 4 Lénin- 





grad, 13, I (1939). 
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5 L. K. Frevel, J. App. Physics 8, 553 (1937). 
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interplanar spacing shifts of the two isomorphs. 
Tables I and II illustrate the practicability of 
the method as applied to the tetragonal system. 

The orthorhombic system presents three 
degrees of complexity : 

(1) Three pinacoid reflections identifiable. By 
selecting a sufficiently large isomorphous series 
(five or more members) it is possible to identify 
the planes (100), (00) and (00/) and the problem 
of indexing the various reflections reduces 
essentially to establishing the correct orders for 
the pinacoid reflections. 

(2) Two pinacoid reflections identifiable. This 
case is indexed in the same manner as outlined 
in the procedure utilizing the temperature shifts 


of the various reflections® of an orthorhombic 


crystal. Table III illustrates the method as 
applied to four isomorphous barium salts. By 
inspection of the spacing shifts it is seen that the 
salts are non-cubic. The tetragonal and hexagonal 
systems are ruled out by the fact that the 
permissible axial ratios, kdox0/l/doo:,8 do not lead 
to a common set of indices for the four patterns 
of the isomorphous barium salts. The next 
system to try is the orthorhombic system. On a 
tetragonal Hull-Davey chart containing only 
(00/) and (O&l) reflections, one adjusts a log d 
strip so that one of the identified pinacoid 
reflections falls on a permissible straight line; 
e.g., (ORO) with R=1V2V3V---, and moves the 
strip parallel to the abscissa until a match is 
found for the other identified pinacoid reflection ; 
e.g., (00/) with /=1V2V3V---. For the proper 
selection of k and / the bc prism reflections, if 
present, will also match. In the case of the four 
barium salts it is found that reflections of type 
(Okl) are restricted to 011, (002 and 020), for 
the first ten observed reflections. The identifica- 
tion of the (AOl) and (hkO) reflections is carried 
out in a comparable manner. The condition 
that the same indices have to fit the four iso- 
morphous patterns is rather stringent and assures 


. a reasonably unique identification of the tabu- 


lated reflections. This conclusion is borne out in 
part by the fact that the starred set of Table III 
does not lead to a common set of indices for the 
four isomorphous salts, although on the basis 


$],k=1V2V3V---. 
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of spacing shifts alone the starred set has to be 
considered as pertaining to a possible pinacoid 
reflection. 

It is also possible to index an orthorhombic 
pattern with only two lattice parameters known 
by solving the following Diophantine equation: 


h od ~?— h id, 
= (h2*k ? —h ke")b ~2 (h ol ? —h 112?)c-, 


where d,; and dz are any independent pair of 
interplanar spacings ; 


hy, he, ki, Re, li, lz=OVIV2V3V---; 


b and ¢ are the identified lattice parameters. 

(3) One pinacoid reflection identifiable; e.g., 
(0k0). The procedure for identifying the prism 
reflections (hkO) and (O&/) in the present case is 
analogous to that for identifying the (hk0) 
reflections of case 2. An interesting alternative 
to this method of picking out prism reflections 
among the various powder reflections is the use 
of a log d-'(Ab/b—Ad/d)' strip’ in conjunction 
with the semi-logarithmic chart containing only 
pinacoid and prism reflections. If the spacing 
shifts Ad/d are sufficiently accurate, it is theo- 
retically possible to identify the indices h and | 
for all of the reflections at one time. The accuracy 
of the interplanar spacing of Table III is only 
0.18 percent and does not permit the use of the 
latter method. 

The proposed methods of indexing isomor- 
phous diffraction patterns have a direct applica- 
tion to the identification of crystalline solids by 
diffraction methods. For each well-established 
crystal structure; e.g., spinel, rutile, wurtzite, 
barite, etc., a weighted log d strip is prepared. 
This collection of strips then serves as a rapid 
reference file for the identification of patterns, 
especially solid solution patterns, that are not 
listed in the published tabulated powder diffrac- 
tion data.2-4 


Ab Ad\} 1/Ab Aa\} 
Tne d-1 a’ tee 8 
log d (; 7) log AG ~*) 


1 ee | 

stoi ee Se 
+3 log t + (ab /b—Aa/a) 
Ab Adoxo | 1/Ab_ Aa\* . o 
where = 7 and log =($ ) is a constant for 
all reflections. 
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Arc Cathodes of Low Current Density at High Amperage 


J. SLEPIAN, W. E. BeRKEy, AND M. J. Koroip 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 
(Received November 12, 1941) 


From photographs, and from the marking on the electrodes, Slepian and Haverstick found 
that short duration arcs of 25 amperes in gases at a pressure of a few centimeters had apparent 
current densities at the cathode of less than 100 amperes per cm?. This is in sharp contradiction 
to the density of several thousand amperes per cm* which seems to be required by generally 
accepted theories of the cold cathode arc. Similar results have been found up to currents of 
5000 amperes, burning for 1/120 sec. with electrodes of various metals in air at less than 
10-cm pressure. Apparent cathode current densities calculated from the observed marking on 
the electrodes were less than 1000 amperes per cm?. The markings on the electrodes were 
apparently produced by oxidation only, the energy density developed at the electrodes being 
insufficient to bring the electrode surface to the melting point in the 1/120 second duration 
of the arc. The circuit interrupting capacity of these short low current density arcs are found 
to be not less than that of short arcs at atmospheric pressure. 


HE fact that arcs with apparently low 

current density cathodes could be formed 
between initially cold electrodes in a low density 
gas was first noted more than 10 years ago by 
Slepian and Haverstick.! In the early work d.c. 
arcs of 25 amperes in Ne, A, and H: at a few 
cm Hg pressure were found by inspection to 
have markings on the cathodes which indicated 
current densities of less than 100 amperes per 
cm*. It was pointed out that with such low 
current densities, the electric field at the cathode 
produced by the space charge of positive ions 
would be insufficient for auto-electronic emission, 
thus making Langmuir’s theory of the cold 
cathode arc untenable.* 

This paper reports experiments in which short 
arcs with currents up to 5000 amperes, lasting 
1/120th of a second are found to have apparent 
current densities less than 1000 amperes per cm? 
in air at pressures below 15 cm Hg. The circuit 
interrupting ability of such low pressure short 
arcs in an alternating-current circuit was found 
by experiment to be as great as at atmospheric 
pressure. 

In one set of experiments arcs were drawn by 
separating current carrying contacts in low 
pressure air. In others, the arc was started by a 
fine fuse wire bridging the already separated 
electrodes. In still others the arc was started 
by an impulsive high voltage breakdown between 


J. Slepian and E. J. Haverstick, Phys. Rev. 33,52 (1929). 
Irving Langmuir, Zeits. f. Physik 46, 282 (1927). 
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the separated electrodes by means of a previously 
charged condenser. A contactor was mounted in 
a bell jar provided with pumps and manometer, 
as shown in Fig. 1. The electrodes were separated 
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Fic. 1. Test apparatus and circuit. 





a short distance, usually one millimeter. On 
110-v a.c. and 220-v a.c. circuits the arc always 
extinguished itself at the first normal current 
zero following the opening. A magnetic oscillo- 
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36 cm Hg 


15 cm Hg 





5 cm Hg 


Fic. 2. Tinned copper. Electrode markings vs. air pressure. One test only 


1 cm Hg 


0.2 cm Hg 
1000 amp., 


for about 0.01 sec. 


graph pictured the voltage and current relations 
during a test. The apparent current density was 
taken as the area which was observably affected 
by the arc as determined by the residual mark- 
ings on the electrodes divided into the maximum 
current, as determined by the oscillograph. 

, The photograph of Fig. 
markings obtained with currents of approxi- 
mately 1000 amperes as a function of the air 
pressure, the arcs being 


2 shows the electrode 


drawn by separating 





contacts. It was found that below approximately 
15 cm the arc markings were very diffuse. The 
cathodes appear to be only slightly etched while 
the anodes do not appear to have lost any metal 
during the arcing. A small molten area in the 
center of several electrodes tested at low pressure 
is attributed to the very high heat 
concentration which explodes away the last 
point of contact of the separating electrodes. 


ohmic 


This type of burning was greatly reduced by 
coating the copper 
with a lead alloy coat- 
ing. This molten spot 
was not found when 
the arc was initiated 
by a fuse wire or a high 
voltage condenser dis- 
charge. A  compari- 
son of electrode burn- 








ing with bare copper, 
tinned copper and sil- 
the 1000- 
was in- 


ver when 


2 aan ampere arc 

Fic. 3. Typical 
markings of 1000- 
amp. arcat low pres- 
sure, (above) arc ini- 
tiated by fuse wire, 
(left) arc started by 
contact separation. 


itiated by contact 
separation and by a 
fine fuse wire is shown 
in Fig. 3. 

Where low apparent 
current densities were 
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Fic. 4. Electrode erosion by repeated tests. (a) After 1000 tests with 1000 amp. —0.61 sec. 
(b) After 8 tests with 5000 amp.—0.01 sec. 


observed, remarkably little melting of the 
electrodes was noted. The markings seemed to 
be an oxidation, with perhaps a slight sputtering 
away of material. This result is consistent with 
calculations of the temperature rise of the 
electrode surface which might be expected in 
the short time from the heat developed by the 
diffuse are. 

Fig.4(a) shows a cathode which has carried arcs 
of 0.01-sec. duration 1000 times with currents of 
1000 amperes r.m.s. In this test the electrodes 
were separated one millimeter and the are was 
initiated by a timed condenser discharge. This 
electrode is 2”. diam.X1/20” thick and only a 
slight erosion can be noted on the cathode 
surface. Fig. 4(b) shows 
electrodes which have car- 
ried 0.01-sec. duration arcs 
eight times with a cur- 
rent of 5000 amperes. The 
are was initiated by a fuse 
wire through the small 
hole in the electrodes. 
These arcs interrupted an 
a.c. inductive circuit of 
180 volts r.m.s. or 254 


Line & 
278 /i™- 


Aaac Tt 
2100 A./lx, 





Vow vrs 


volts crest. 

Typical oscillograms of 
tests are presented in Fig. 
5. Test L.P. 80 shows volt- 
age and current record 
of a 1000-ampere current 
test. L.P. 89 shows an 
interruption test at 220 
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Fic. 5b. Typical oscillogram 
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volts where the current is limited to 530 amperes 
crest by an air core reactor. 

That the immediately regained dielectric 
strength of an arc space at a normal current 
zero in a short a.c. arc was practically equal to 
the normal cathode drop of a glow has been 
pointed out by Slepian.* The normal cathode 
voltage of a glow is very nearly equal to the 
minimum breakdown potential between plane 
electrodes of the gas; and both are independent 
of current and gas pressure for the same gas and 
electrode material. Thus, theory predicts then 
that the a.c. arc extinction voltage for a short 
arc should not be lowered by reducing the gas 
pressure. The expected manner of recovery of 
dielectric strength after current zero in a short 
gap arc in low pressure is as shown in Fig. 6. 
After the rapid initial recovery the rate of gain 
will be faster up to the breakdown in a low 
pressure arc space due to faster ion diffusion to 
the electrodes. This expectation has _ been 
confirmed by experiments. 

The well-known variation of breakdown volt- 
age with the gas density for a given gap spacing 
suggests that by reducing the breakdown voltage 
to somewhere near the minimum by reducing 
the gas density, a low ratio of breakdown to a.c. 
are extinction voltage would be obtained. Such 
a gap could be used for overvoltage protection 
of electrical apparatus where the longer time 


3 J. Slepian, Trans. A. I. E. E. 47, 1398 (1928). 
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Fic. 6. Reignition voltage vs. time for short arc on 
a.c. test; A—atmospheric pressure, B—reduced air pressure 
(in order of 0.1 to 15 cm Hg). 


lag of breakdown obtained in reduced gas 
density is not objectionable. 

It has been found by experiment that: 

(1) The electrode surface destruction by an 
arc with a.c. currents of several thousand 
amperes lasting 1/120 seconds is greatly reduced 
at low air pressures as compared with normal 
atmospheric pressure. 

(2) The a.c. are extinction voltage for a short 
gap at low air pressure is not less than that of 
the same gap at atmospheric pressure, as 
predicted by theory. 

(3) The actual ratio of breakdown to arc 
extinction voltage can be made close to unity in 
a low pressure gap and utilized as an overvoltage 
protection device. 





Erratum: Determination of the Static and Dynamic 
Constants by Means of Response Curves 


[J. App. Phys. 12, 866 (1941) ] 


R. K. BERNHARD 


The Pennsylvania State College, State College, Pennsylvania 


4 I ‘HE equation which appears at the bottom of page 868, 


(4mc — 2B err — on =0, 


should read instead as follows: 


(4mc—2k Jw 
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Operational Analysis of Non-Linear Dynamical Systems 


Louis A. PIPEs 
Graduate School of Engineering, Harvard University, Cambridge, Massachusetts 


(Received November 24, 1941) 


An operational method of determining the free and forced vibrations of non-linear dynamical 
systems is presented. The method is an operational adaptation of the one developed by Lindsted 
and Liapounoff. The advantage of the operational procedure is that by using a table of opera- 
tional equivalents, one may write the various successive approximations almost at a glance. 
The procedure is illustrated by applications to non-linear problems in mechanics and electricity. 


I. INTRODUCTION 


HE beautiful simplicity of the mathematical 

theory of electrical circuits and small oscil- 
lations of mechanical systems is a consequence 
of the fact that the Lagrangian equations of 
motion of these systems are linear differential 
equations with constant coefficients. When, how- 
ever, the mechanical system under consideration 
is undergoing vibrations of large amplitudes or 
the electrical circuit contains an iron-cored in- 
ductance coil, the mathematical formulation of 
the problem leads to non-linear differential 
equations. 

The solution of these non-linear equations is an 
extremely complicated task. The case of free 
vibrations may be solved in terms of elliptic 
functions.' Usually the differential equation can 
only be integrated graphically or numerically. 
These methods are laborious and give the solu- 
tion of only the particular system under con- 
sideration, and the solution cannot be generalized 
to other similar systems having different nu- 
merical constants. 

Methods of successive approximations have 
been developed which are of great power in 
studying systems deviating but slightly from the 
ideal linear case.2~* The method presented in this 
paper is an operational adaptation of the one 
developed by Lindsted and Liapounoff*? and 
involves a series expansion in a parameter. 


' George Duffing, Erzwungene Schwingungen bei verdnder- 
licher Eigenfrequenz (Braunschweig, 1918). 

2A. M. Kriloff, Bulletin of the Russian Academy of 
Sciences, 1933, No. 1, p. 1. 

Lord Rayleigh, Theory of Sound (1894), Vol. 1, p. 77. 

*C. Shaefer, Ann. d. Physik 33, 1216 (1910). 

°H. Poincaré, Mecanique Celeste, Vol. I, p. 58. 

°L. A. Pipes, “An operational treatment of non-linear 
dynamical systems,”’ J. Acous. Soc. Am. 10, 29 (193 8). 
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Il. THE LAPLACIAN TRANSFORMATION 


Before proceeding, let us review briefly the 
theory on which the modern operational calculus 
is based. The modern development of operational 
calculus is based upon the Fourier-Mellin integral 
equations. This set of equations state that if: 


g(p) =p e~?'*h(t)dt, (1) 
then : 


c—jx p 


h(t) =—— 


1 “+i” g(p)e”'dp 
f - where j=\/-—1 (2) 
2rj 


and ¢ is a real constant of such magnitude that 
the poles of the integrand lie entirely to the left 
of the path of integration. For brevity, let us 
use the notation: 


Lh(t) =g(p) (3) 
and 


L-'g(p) =h(t) (4) 


to express symbolically the relations (1) and (2). 

Many interesting and important theorems 
have been established’ which make these integral 
equations powerful tools in the solution of 
dynamical problems. We turn now to a con- 
sideration of the manner in which this powerful 
method of analysis may be used to solve certain 
types of non-linear differential equations. 


III. FREE VIBRATIONS OF A NON-LINEAR 
OSCILLATOR 


The method may be illustrated by a considera- 
tion of a mechanical oscillating system consisting 
of a mass attached to a spring. The equation of 


7L. A. Pipes, ‘The operational calculus,”’ J. App. Phys. 
10, 172, 258, 301 (1939). 
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the free vibration of such a system is 


md*x /dt*+ F(x) =0, (5) 


where md*x/df is the inertia force of the mass, 
F(x) is the spring force, and x is measured from 
the position of equilibrium of the mass when the 
spring is not stressed. Let us consider the sym- 
metrical case where 


F(x) =kx+ bx’. (6) 


Hence (5) becomes 


md*x dt? +kx+bx* =0 (7) 
or 
d*x/ dt? +w*x+ax*=0, (8) 
where 
w=(k/m)', (9) 
a=b/m. (10) 


Equation (8) occurs in the theory of non-linear 
vibrating systems and certain types of non-linear 
electrical systems and serves to illustrate the 
general method of analysis. Let us consider the 
solution of the Eq. (8) subject to the initial 
conditions : 


(11) 


That is, the mass is displaced a distance a 
and allowed to oscillate freely. We are interested 
in studying the subsequent behavior of the 
motion. 

Let us now multiply the Eq. (8) by pe-?‘dt 
and integrate from 0 to «. We thus obtain: 





rf e?'[ d*x/ dl? +w*x+ax' |\dt=0. (12) 


0 


Let us now use the notation of Section I, and 
write: 


Lx(t) = y(p) (13) 


Now by an integration by parts, it may easily 
be shown that: 


Ld?x dt? = — pio— p*xo+ p’y, (14) 


where #» and x» are the initial velocity and dis- 
placement of the particle at ‘=0. 

In view of the initial conditions (11) and the 
transform (14), the Eq. (12) may be written in 
the form: 


(p?+w*)y = p*a—aLx’. (15) 

Now let 
X=Xop tax; +a*x.+a*x3+---, (16) 
w?=wo’ +Ciatcoa*+c30'+---+, (17) 
Y=Yotayita®yetaryst+, (18) 
yi(p)=Lx(t). (19) 


In these expressions the quantities x,(t) are 
functions of time to be determined, wo is the 
frequency which will be determined later. The c 
quantities are constants which are chosen to 
eliminate resonance conditions in a manner that 
will become clear as we proceed. The y,(p) 
quantities are the Laplace transforms of the x,(t) 
quantities. 


In most non-linear dynamical systems, the quantity a is small compared to w* and the series (16) 
may be shown to converge. In the following discussion, let us limit our calculations by omitting all 
the terms containing a to a power higher than the third. Substituting the above expressions into (15), 


we obtain 


bP? (yotayita*y.+a*ys) + (wo? +61:a+ C207 +630°) (yo tayita*y2+a*ys) 


(p* vot woo) Fal p*yi two? Vi tCivotLxo*) + a?(p2y2+wo7V2+CoVo + Civ + L3x07x1) 


= p*a —_ aL (xp tax;t+a*xe+a*x;)*. (20) 
If we now neglect all terms containing a to powers higher than the third, we obtain 
+03(p*y3t+worvstcsvotcovitcivetL[ 3x02x2+3x0x1? ]) = pra. (21) 


This equation must hold for any value of the quantity a. This means that each factor for each of 
the three powers of a must be zero. Hence the Eq. (21) splits up into the following system 
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of equations: 


b?yy + wo2Vo = pra, (22) 
p*yitwo"*yv1 = —Ciyo— Lxy*, (23) 
Pp? Vet wo*V2= —CoVo— C11 — L3x0*x1, (24) 
p?va+wo2V3 = —C3Vo— C2V1 — C12 — L[ 3x02x2+ 3x 0x1? |. (25) 

Using the notation, ‘ 
T(¢) =1/(p?+¢°), (26) 


then the Eq. (22) may be written in the form: 
Vo= p?aT (wo) = Lxo. (27) 
From the table of transforms, see Appendix I we have 
L-'p*aT (wo) =a cos(wot) = Xo. (28) 


This represents the first approximation to the solution of the Eq. (8) subject to the initial condi- 
tions (11). The transform of the second approximation as given by (23) may be written in the form: 


Vi = —C1Vo0T (wo) — T (wo) Lxo’. (29) 


From the table of transforms, we have: 
a? 
Lx,* = La’ cos? (wot) =— [ 3p? T (wo) + p*T( 3w) ]. ( 30) 
4 


Substituting (27) and (30) into (29), we obtain 
3a*] a’ ; 
= = pTa0)|c0 + == | - p?T (wo) T( 3wo). (31 ) 
4 4 


Now from the table of transforms it is seen that 
} & \p?T?(wo) = sin (wot) 2w. (32) 


Hence the first term of the right member of (31) corresponds to a condition of resonance. We may 
eliminate this condition of resonance by placing the coefficient of this term equal to zero. Then 


[c.a+3a*/4]=0 (33) 


or 


¢,= —3a?/4. (34) 
This determines the constant c;. With the resonance condition eliminated, (31) reduces to 
a’ 


u=— {P*T (wn) T (30). (35) 


Making use of the table of transforms, we obtain the second approximation 


a* a’ 


x= L~'y, = L~'— pT (wo) T (30) =- [ cos(3wot) —cos(wot) }. (36) 
4 2 


Wo 
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If we limit our calculations to the second approximation, we obtain from (16) 


aa* 
x =a cos(wol)+ | cos(3wol) —Cos(wol) |. (37) 
32wy" 


The angular frequency is obtained by substituting the value of c; given by (34) into (17). This 
gives 


wo? = w?+ fa°a. (38) 


From this we see that the presence of the term x* in the equation introduces a higher harmonic 
term cos (3wof) and the fundamental frequency is not constant but depends on the amplitude a, 
and it increases with a provided that the quantity a@ is positive. 

The third approximation is obtained by substituting the above values of yo, v1, xo, and x; into (24). 
This gives: 


c,a* 
V3™= —Cop*a T?(wo)+ ~-T (wo) p? T (wo) T(3wo) = T (wo) L(3x07x1). (39) 
4 
We must now compute 
3a°® | 
L(3x07x1) = L: COS*(wot) Lcos(3wot) —cos(wol) |}. (40) 
| 32w 9" 


By using the table of transforms, we easily obtain 


3a°*p* 
L(3x07x1) => -LT(Swo) +7 3) —2T (wo) |. (41) 


2wy" 
Substituting this value of L(3xo?x,) and making use of the identity 


1 
T(a)T(b) = -[T(a) —T(b)], (42) 


(b?—a?) 
we write (39) in the form 
c,a* 3a® 


C\a : 3a® 
y= PTHu)| — c+ +— |+e7@) 7600] - —- 
32w 9" 64a" 32w 9" 4.320" 


3a® 
a p*T(an)T(Su0)| ——] (43) 
4.329 


To eliminate the condition of resonance, we equate the coefficient of the p?7*(wo) term to zero, 
substituting the value of c, into the coefficient. We thus obtain 


Co= 3a‘ 128’. (44) 
On substituting the value of c; into the second member of (43), we see that this term vanishes and 
we have 
3a°® : , 
y= —p*7 (wo) T(5wo). (45) 
1289" 


Using the table of transforms to obtain the inverse transform of y2, we have the third approximation 


a® 
X2 = —————[_cos(5wof) — cos(wof) ]. (46) 
10249‘ 
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From (16) we thus have the third approximation 


aa* a*a® . 
x =a COs(wot) +——[_cos(3w!) — cos(wo!) ]+ —[cos(5w!) —cos(wo!) ], (47) 


Wo 1024 





where now the fundamental frequency is given by (17) to be 
wo”? = w+ 2a?a— 3ata?/128a9?. (48) 
The fourth approximation is obtained in the same manner from (25). We compute 


Tp2 


L(3x07?X2+ 3x 0x1") aaaipnaidiits - [2 T (7) + T(5w) —_ 3 zt 3wo) | (49) 
4.1024 ' 


by using the table of transforms. Substituting the values of the quantities yo, y1, and y2 as given above 
into Eq. (25) and making use of the relation (49), we obtain: 


CoT (wo)a* p* C11 (wo)a*p? 
ys = —c3p*aT*(wo) a nneemnunecemmee [ T( 3wo) ~~ T (wo) |- Les —[T(Swo) as T(wo) | 
\ 32w" 1024w9' 
3a’ p* T (wo) 
——————[2T(7w0) + T(5wo) —3T(3wo) ]. (50) 
4-1024w' 
The condition for no resonance leads to 
C3 = Co0* / 32w9? + 30° / 10249". (51) 


Substituting the values of c; and cz given by (34) and (44) into (51) we obtain 


cz; = 0. (52) 
Suppressing the resonance terms, the Eq. (50) reduces to 
3a? 6a’ , 
V3= p?T (wo) r(3w.)| |+p-7%e0 T( ia»)| — cence | (53) 
2048 be 4096w4 


Computing the inverse transform by the use of the tables of transforms, we obtain: 


a‘ 
X3 —[5 cos(wot) — 18 cos(3wot) +cos(7wol) }. (54) 


~ 32.7680" 


Substituting xo, x1, X2, X¥3 into (16) we obtain the fourth approximation 


aa* aa? 
XxX =a COS(wol)+ —[cos(3wot) —cos(wot) |+ —[cos( Swot) —cos(wot) |] 
32w 9" 1024w! 
ata? 
+————5 cos(wot) —18 cos(3wot)+cos(7wot) ]. (55) 
32,768w9° 


To this approximation, the fundamental frequency wo is given by (17) and is 
3 3 a‘ 
wo? = w*? +-—aa? — —a*—-. (56) 
128 wy” 


Now in all the calculations, terms to higher power than the third have been omitted. We may 
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simplify Eq. (56) by substituting on the right side the value of wo given by (38). We thus obtain : 


3 a‘ 
of | (57) 
128 Lw?+?aa?. 


Expanding the bracket term in powers of a and retaining only powers of a up to the third, we have 


3 


Wo” = w* +—aa* — 


4 


3 


wo”? = w* +-aa* — 
4 


3a* a‘ 


128 w? 


Jara® 
—., (58) 
512w? 


Further approximations may be carried out by the same general procedure. 


IV. THE VIBRATIONS OF A PENDULUM 


The above analysis may be applied to the 
study of a theoretical pendulum. The equation of 
motion of such a pendulum is 


54° sind =0, (59) 
where / is the length of the pendulum and g is 
the gravitational constant. If we develop sin @ in 
a power series in @ and retain only the first two 
terms of the series, we obtain: 


U 


Tr 
6+-6——#=0. (60) 
l 


6/ 

If we stipulate the initial condition that at 
t=0, 6=6@), we may make use of the preceding 
analysis by letting 


6,=a, (61) 
w=g/l, (62) 
a=-—g/6l. (63) 


Using the first approximation for the angular 
frequency as given by (38), we have 
wo = 2g l—g0,° Sl. (64) 


To this approximation, the period of the 
oscillations is given by : 


rae") 


For small amplitude 4) the radical may be 
expanded in powers of #) and we may write 


ly} 6" 
T+ 2r( ) (+ ). 
g 16 


(65) 


(66) 


This formula gives excellent results for small 
amplitudes of oscillation. 


V. NON-LINEAR ELECTRICAL CIRCUIT 


As another example of the method, let us 
consider the free oscillations of an electric circuit 
consisting of an inductance and a condenser in 
series. The equation of such a circuit may be 
written in the form 


do@‘dt+Q/c=0, (67) 


where @ is the total flux linking the circuit, 
Q is the charge on the condenser, and c is the 
capacitance of the system. For simplicity we 
assume that the circuit is devoid of resistance. 
If the inductance of the system consists of a coil 
of wire wound on an iron core, we have in 
general 


i= F(9), (68) 


where 7 is the current flowing in the circuit. 
The functional dependence between the current 
in the circuit and the flux @ depends on the 
functional variation of the flux density, B, of 
the material of the inductance coil with the 
magnetizing force, which is proportional to the 
current 7. This function is many-valued in 
general. However, it may be taken to be single- 
valued at large magnetizing forces for Nicalloy, 
Permalloy, and low-loss steels. 

Several analytical expressions have been sug- 
gested to express approximately the relation 
between the magnetizing current 7 and the total 
flux @ in a ferromagnetic material. Some of the 
expressions are 


m 


t=) a.¢", (69) 
n=1 
1=A sinh(¢), (70) 
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D x 

i=Agt+>. A, cos(nd)+)>. B, sin(nd) (71) 
n 1 n=1 

and other forms may be derived. In ordinary 

circuit theory, the linear relation 


1=—o (72) 


is assumed. The coefficient L is then the in- 
ductance of the system. In our analysis, we will 
assume that the relation between i and ¢ is 
given by 

1 1 


1i=—¢+-¢', (73) 
Lb 


where the coefficient 1/b is a measure of the non- 
linearity of the system. This expression may be 
regarded as the first two terms of the power 
series of the relation (70), and by a proper 
choice of the coefficients may be made to approxi- 
mate the actual i—@ curve of a material quite 
accurately. If we differentiate Eq. (67), we 
obtain 

d*/dt?+1/c=0 (74) 


in view of the relation 
1=dQ/dt. (75) 


Substituting the value of 7 given by (73) into 
(74) gives 
d*o@ 1 | 
ae a ae (76) 
dt? Le cb 


If we now let 


2=1/Le, (77) 
a=1/be, (78) 
p=x. (79) 


Then the Eq. (76) is identical with Eq. (8). 
Let us now consider the case in which at ¢=0 
there is no current flowing in the system and 
the condenser has an initial charge of Q». In 
view of Eqs. (67) and (73) we therefore have the 
following initial conditions of the independent 
variable x. 

x=0, 


| at t=0. (80) 
dx/dt= —Qo C=Vo0 


In this case, therefore, we have to solve the 
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Eq. (8) subject to the initial conditions (80). 
As before we let y be the transform of x. The 
transformed equation is in this case 


(p?+w*)y = pop —aLx’. (81) 


To treat this case, we again establish the 
sequences (16)—(19) and by equating like powers 
of a, we obtain the set of equations 


p*yot wo? vo= pro, (82) 
p*yvitwoyi= —Ciyo— Lx’, (83) 
p*yv2+wo’V2= — CoVo— C191 — L(3x07x1), (84) 
p*vstwo'vs = —C€3Vo— CoVi-—Ciye2 


— L[3x0?x2+3x9x1? |. (85) 


This set of equations is the same as the set 
(22)-(25) with the exception of the first one 
because of the different boundary conditions. 
The first approximation is obtained from (82) by 
writing 

Vo= pvoT (wo) = Lx». (86) 


From the table of transforms this gives the 
first approximation 


Xo =V06 sin (wot) Wo. (87) 
We now write (83) in the form 
\i= T (wo) =-— fae Lxo* |. : (88) 
From the table of transforms, we have: 
3vo* 
Lx,* = -LpT (wo) — pT(3wo) |. (89) 


Wo 


Substituting (86) and (89) into (88), we 
obtain 


39° 
i= pT] =— ie | 


4u” 


3v9* 
oo pT (wo) T( 3wo) = Dx). (90) 
4w9" 


The term p7?(wo) leads to resonance, as may 
be seen from the table of transforms; hence its 
coefficient must be set equal to zero. This yields 


C1= — 3u9" 49”. (91) 


By means of the table of transforms (Ap- 
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pendix I), we then obtain 


35° 
mn=Ly,= ——[sin (wot) —sin(3wof) J. (92) 
32w° 


To this approximation, we obtain the angular 


frequency, wo, from (17) to be 


wo” = w* + 3r97a/ 4”. (93) 
The second approximation is therefore : 
Vo SIN(wol)  3u9%a sin (3wot) 
x= + Sin (wot) — —j. (94) 
Wo 32w° 


The variation of the charge on the condenser 


may be obtained from (67) in the form 


dx 
O= —c— = Qy cos(wol) 
dt 


30," 


+ [ cos(wo!) —cos(3wof) | 
32wo'bc* 


}. (95) 
In view of (93) the natural frequency of the 
system may be written in the form 


lpi 3Q°L?/ 
fau—1—+ | (96) 
2milc 4 be 
to the second approximation. We see that the 
frequency varies with the initial charge on the 
condenser, Qo. 


VI. OSCILLATOR NON-LINEARLY DAMPED 


As another application of the method in the 
study of free oscillations, let us apply it to the 
study of a freely vibrating system whose damping 
force is proportional to the square of the velocity. 
The equation of motion of this system is given by 

¥+w*xFar*=0. (97) 
taken when the 
velocity is in the direction of the negative x 
axis and the plus sign when the velocity is in 
the direction of the positive x axis. Let us assume 
the initial conditions 


The minus sign must be 


x=a 
at ?=0 


(98) 
z=) 
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We thus have to consider the equation 


x+w*x — at*=0 (99) 
for the first half of the oscillation. As before, we 
set up the sequences (16)—(19), and if we limit 
our calculations to the terms containing a’, we 
obtain the following set of equations: 


p* vot wo Vo = p*a, (100) 
p*vitwo7y = —CiWot Lio’, (101) 
p*y2t+wo' V2 =—¢1¥i—- Covot L( 2% 5X1) " ( 102) 
From (100) we obtain 
Yo= p*’aT (wo) = La cos(wot)= Lx. (103) 
Substituting this into (101), we get 
Vi=--— C\p*a T (wo) 
wa" 
+ [ T (wo) — p?T (wo) T( 2wo) |. (104) 
2 
In order to eliminate resonance we must 
choose 
c,=0. (105) 
Hence 
wrt Pn oe ll, 
Vi= sai T (wo) — T(wo)+ T (2a) ° (106) 
2 30" Sw" a 


By means of the table of transforms we 
obtain 
a? 2 a* - 
x1 =L-'y,;=— ——a?* cos(wolt) +— cos(2wol). (107) 
3 6 


Substituting the various transforms involved 
into (102) gives 


2a Sao” 
y2= pT(w0| —¢cat+— 
6 
2a 4ey9” p* 


idscseneeiianintal -T (wo) T (39) 
6 


4a?" 4a*w°" 


+—— p?T (wo) T(2wo) —-———T (wo). (108) 
6 6 


The condition that there be no resonance in 
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this case gives 
Co= a7 wo", 3. 


(109) 


From the table of transforms, we obtain 

2 a* . 
x2=L~'y.= —-a*+—[61 cos(wot) 
3 72 


— 16 cos(2wot) +3 cos(3wot) ]. (110) 


Substituting the above values of xo, x1, and x» 
into (16), we obtain 


aa? 
x =a COS(wot) +—(3 —4 cos wot+cos 2wof) 
6 


aa* 


—(48 — 61 cos wot +16 cos 2wol 
72 


—3 cos 3wot), (111) 


Substituting the above values of c; and ¢2 into 
(17) gives 








wa" 
2 = wo? --+-—____ (112) 
3 
or 
WwW 
oo a (113) 
(1+4a%a?)! 


If we let 7,/2 be the time required for the 
oscillation to reach its maximum displacement to 
the left, we have 

T, w w(it+ta’a?)! 
— =— =-——______ + (14 }a’a’). 


2 Wo Ww Ww 


(114) 


To determine the maximum negative displace- 
ment, +a;, we place wo =z in (111) and we 
obtain 

4 16 
+a,= —a+-aa*——a’a’*. 
3 9 


(115) 


Now to determine the motion for the next 
half-cycle we must solve the equation 


¥+w*x+azr=0 (116) 
subject to the initial conditions 
x=—dy 
at (=0. (117) 
z=0 


However, if we make the charge in variable 


x=-Yy, (118) 
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we see that we must then solve the equation 


ytw*y—ay?=0 (119) 
subject to 
yr 
t=0. 
y=0 (120) 


However, this is the same equation whose 
solution we have just obtained. The only differ- 
ence is that we have a, instead of a in the initial 
conditions. We thus obtain 


T2 Tv 
5 Tt bata) 


ie) 


(121) 


for the time required for the second half-cycle, 
and 
4 16 


y= —A,+—-aa;* ——a’*a;' 


9 


(122) 


for the displacement of the system at the end 
of the second half-cycle. It is thus evident that 
the oscillation has a gradually decreasing ampli- 
tude. Its amplitude after the (n+1)th half cycle 
is given by 

4 16 


An4i= —d,+—-ad,?— 
3 





a*a,,*. 


(123) 


VII. FORCED NON-LINEAR VIBRATIONS 


The above method may be applied to obtain 
the steady-state solution of a non-linear forced 
system. As an example, consider a forced oscil- 
lating system whose equation of motion is 
given by 


¥+w*x+ax*=q cos(bt). (124) 


This is the equation of motion of a non-linear 
system that has a disturbing force proportional 
to cos (bt) impressed upon it. This is a non- 
linear equation whose general solution is not 
known. Let us consider the steady forced vibra- 
tions and disregard the free vibrations that 
depend on the initial conditions. It will be 
assumed that a@ is small and that the resulting 
steady forced vibrations may be represented as a 
periodic function of the time having a funda- 
mental frequency of 6/2 and is therefore repre- 
sentable as a Fourier series in multiples of b. 

To carry out the operational solution, we again 
set up the sequences (16)—(19) and make use of 
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the transform 


Lq cos(bt) = qp*/ (p?+?) (125) 
and that the steady forced vibration is built up 
of an amplitude a in phase with the disturbing 
force. Substituting the above sequences and (125) 
into the Eq. (124) and equating the coefficients 
of the various powers of a in the resulting ex- 
pression, we obtain the following set of equa- 
tions: 


(p?+ wo") vo= p’at+qp?/ (p?+b*), (126) 
(p?+wo")y1 = —Ciyo— Lx’, a8) 
(p?+wo*)¥2= —Co¥o—Ciyi—L3xo*x1, (128) 


where we have assumed that at t=0, x=a, and 
we limit our calculations to terms in a®. From 
(126) we obtain 


Vo= p?aT (wo) +qp*T (wo) T(d). (129) 


From the table of transforms, the inverse of 
this expression is 


q 
xo= L yo = @= cos Wot 
wo”? — b? 


q 
+. cos bt. 


(wo? — b?) 


(130) 


In order that the vibration have a fundamental 
frequency of b/2x, we must place 


a=q/(wo?—b*). (131) 


Hence we obtain 


x»9=a cos bt. (132) 


Substituting this expression into (127) and 
solving for y;, we obtain after some reductions: 


p | 
N= [7(b) — T(wo) JL —c1a — 3a*} 


(wo" — b?) 


a® p* 
4 (wo — 9b?) 


[7(3b)—T (wo) }. (133) 


In order that the vibration have a fundamental 
frequency of b/2z, the coefficient of the term 
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p> T' (wo) must be equated to zero. This leads to 


ca+ fa a* 
-—- =: () (134) 
(wo" — b?) 4(w»? — 9b?) 
or 
3 a” (wo? — Db?) 
c¢y= —-a*+ (135) 


4 4 (9b? — wo?) 


Substituting this value of c,; into (133), we 


obtain 
a* 
Vi= 


[ p?T(b) — p2(3b) }. (136) 


4(w” — 9?) 
Hence 
a’ 
m=L-'y,= [cos bt —cos 3bt }. 
4(wy" —_ Ob?) 


(137) 


Hence to this degree of approximation we have 
x=Xy+ax,=a cos bt 
aa* 


[cos bt—cos 3bt }. 
4( Wo" ae 9b? ) 


(138) 
To determine the amplitude a we make use 
of the equation 
w* = wo? +c1a (139) 
and substitute the relation (131) in the form 
wo =g/atb? (140) 


and the value of c; as given by (135). We thus 
obtain 


aa" 1 


q 3 
a + 4 (1—867a/q) 


Multiplying both members of (141) by a, we 
obtain 


(141) 


1 
v*ataa?=9+b'a+ taal 1 ~ ), (142) 
1—8b*a/q 


Now in most vibrating systems the condition 


b*a>q (143) 
is satisfied for large amplitudes a and (142) may 
be written in the form 

w*a+aa*=q+b*a+aa*/4 (144) 
or 


taa*+(w?—b*)a—gq=0. (145) 
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FIG. 


This is a cubic equation in a. It may be solved 
graphically by plotting 


yi(a) = faa’, 


(146) 


yv2(a) =q+(b?—w*)a, (147) 


as two functions of a as in Fig. 1. 

The intersections of the curves y; and ye give 
the possible values of the amplitude a. It is 
seen from the figure that for small values of } 
there is only one point of intersection and we 
have only one amplitude of the forced vibration. 
However, as b increases, @ increases until for a 
certain critical frequency }» there are three inter- 
sections. It is therefore possible for the vibration 
to undergo abrupt changes in amplitude. The 
significance of the three different solutions was 
studied experimentally by O. Martienssen,* using 
a non-linear electrical circuit, and the stability 
of the three possible solutions was discussed 
theoretically by E. V. Appleton.® 


VIII. CONCLUSION 


A consideration of the above examples shows 
the practical utility of the operational method in 


5Q. Martienssen, Physik. Zeits. 11, 448 (1910). 
®E. V. Appleton, ‘The motion of a vibration galva- 
nometer,” Phil. Mag. [6] 47, 609 (1924). 
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the analysis of non-linear systems. By using the 
table of transforms and the shorthand notation, 
the tedious effort of obtaining the various suc- 
cessive approximations is reduced to a minimum 
and the chance of error is vastly diminished. 


APPENDIX I 
Table of Transforms 


Let 7(w) = > 
p? +? 

1. L sin(wt) =wpT(w) 
2. Lsin*(wt) = 3[1—p?T(2w) ] 
3. L sin?(wt) = }[3pwT (w) —3wpT (3w) ] 
4. L cos(wt) = p?T (w) 
5. L cos*(wt) = [1+ p?7(2w) J 
6. L cos*(wt) = }[3p*T(w) + p*T (3) J 
7. L sin(At) cos(Bt) = }[(A+B)pT(A+B) 


+(A—B)pT(A—B)] 


8. L cos(At) cos(Bt) =" [T(A+B)+T(A—B)] 

9. L sin(At) sin(Bt) =" [T(A—B)—T(A+B)] 

10. L sin(At) sin(Bt) sin(ct) ="[(A+B—0)T(A+B—0) 
+(B+c—A)T(B+c—A)+(c+A—B)T(c+A—B) 
—(A+B+c)T(A+B+c)] 

11. L sin(At) cos(Bt) cos(ct) AG +B—c)T(A+B-—c) 


—(B+c—A)T(B+c—A)+(c+A —B)T(c+A—B) 
+(A+B+c)T(A+B+c)] 
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»” 


F[T(B+e—A) 


12. L:sin(Al) sin (Bt)cos(ct) = 
—T(A+B—c)+T(c+A—B)—T(A+B-+c)] 
13. L cos(At) cos(Bt) cos(ct) =" (ri +B—c) 


+7(B+c—A)+T(c+A—B)+T(A+B+<c) ] 
cos(at) —cos(bt) 


14. L 'p*T (a)T(b) = (b?—a?) 





- t sin(wt) 
15. L'p*T*(w) =—— 
2 


16.1 pT (us) = net) _ # cos (ut) 





2w 2u:? 
- _ . 1—cos(wt) 
17. L“T(w)=——— 
18. T(a)T(b) = ————_[T(a) — T(b)]. 


(b?—a?) 





Note on Temperature Correction Methods in Calorimetry 


R. S. Jessup 
National Bureau of Standards, Washington, D. C. 
(Received November 17, 1941) 


The problem of correcting calorimetric data for thermal leakage has been discussed by 
King and Grover in a recent paper in this Journal. From their analysis of the problem these 
authors have concluded that large errors are involved in the usual methods of correcting the 
data of bomb-calorimetric measurements, and of experiments by the method of mixtures. 
In the present note it is shown, in agreement with W. P. White’s treatment of calorimetric 
lag, that in the case of bomb-calorimetric measurements the errors in question are practically 
completely eliminated by calibrating the calorimeter experimentally. Since experimental 
calibration of the calorimeter is the usual practice in modern bomb-calorimetric measurements, 
the results of such measurements are not affected by the errors discussed by King and Grover. 
In the case of measurements by the method of mixtures these errors can be avoided by so 
conducting the experiment that the final. temperature of the calorimeter will be very near 


the convergence temperature. 


: a recent paper in this Journal, King and 
Grover' have presented a mathematical 
treatment of the problem of correcting calori- 
metric data for thermal leakage. They conclude 
from their analysis of the problem that large 
errors are introduced into the results of bomb- 
calorimetric measurements, and of measurements 
by the method of mixtures, by the usual methods 
of calculating the correction for thermal leakage. 
For a particular hypothetical experiment they 
report that the usual methods of calculation lead 
to errors of the order of 0.2 percent. As the 
precision of modern bomb-calorimetric measure- 
ments is 0.01 or 0.02 percent, an error of 0.2 
percent in the results of such measurements 
would be very serious. In the present note it 
will be shown that in the case of bomb-calori- 
metric measurements, the errors discussed by 
King and Grover are practically completely 
eliminated if the calorimeter is calibrated by 
supplying a measured quantity of energy and 

‘Allen King and Horace Grover, ‘“Temperature cor- 


rection methods in calorimetry,”’ J. App. Phys. 12, 557 
(1941). 
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observing the resulting temperature rise, and if 
the usual methods of correcting for thermal 
leakage are used for both the calibration experi- 
ments and the measurements of heats of com- 
bustion. It will be shown also that in the case of 
experiments by the method of mixtures, errors 
resulting from the usual methods of correcting 
for thermal leakage can be avoided by a suitable 
experimental procedure. 

The problem discussed by King and Grover 
is that of a calorimeter in which is immersed a 
charge, such as a bomb. It is assumed that at a 
time fo =0 the charge is at a uniform temperature 
©», and the calorimeter is at a lower uniform 
temperature @)=0. It is further assumed that 
the temperatures © and @ of charge and calo- 
rimeter, respectively, vary with time in a manner 
which is described by the differential equations 


d@/dt= —a(O—6) (33)? 


* Equations taken from reference 1 are given the same 
numbers (expressed in Arabic numerals) as in that paper. 
The numbers of other equations given in this paper are 
expressed in roman numerals. 
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and 


dé/dt=b(Q@—6) —c(@—98,), (34) 


where a, 6, and c are constants of the apparatus, 
and 6, is the convergence temperature of the 
calorimetric system, that is, the temperature 
which the system would attain in an infinite 
time if the rate of stirring and all external 
conditions remained constant. 

Elimination of © from Eqs. (33) and (34) 
vields the second-order differential equation 


d*0/dt?+2 pd0/dt +q0=q6x., (35) 


where 2p=(a+6+c) and g=ac. The second- 
order differential equation in © is of the same 
form as (35). The solutions of these differential 
equations may be written in the form 


6=6.+Ae-'+ Be-@r-0t (36) 
O=6.+ Ce + De?!" (1) 


and 


where «= p—(p*—q)!, and A, B, C, and D are 
constants of integration. Insertion of the nu- 
merical data given by King and Grover into 
Eq. (36) leads to the following expression for 
the time-temperature relation for the hypo- 
thetical experiment 


8= 1.7082 +1.010029e—°-0 


—2.718229e-°-**, (IT) 


It can be shown that the temperature rise of 
the calorimeter corrected for thermal leakage is 
given by 
" —€A —(2p—6€)B+c(09— 6x) 
0..°—~@-—— Eh OR 9 SIR 2 
a+b 
The value of (6..*— 6) calculated from Eq. (III), 
using the values of the constants given by King 
and Grover is 2.69302°.8 
It should be noted that although the mathe- 
matical treatment of the problem under con- 
sideration by King and Grover leads to the 
correct value of the temperature 6,* which the 
calorimeter would attain if there were no thermal 
leakage, this temperature can be calculated by 
means of Eq. (III), or the equivalent equation 


(III) 





3 King and Grover give the value 2.6931° for (0..*—o). 
The difference may be due to the use of rounded values 
of some of the constants in the present paper. The value 
given in the text is consistent with Eq. (II) and the 
values of (a+5) and ¢ given in reference 1. 
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given by King and Grover, only when the time- 
temperature relation is of the form given by 
Eq. (36). Experimental time-temperature curves 
usually cannot be represented by equations of 
this form. After the ignition of a sample of 
combustible material in a bomb, or after a 
heated charge is dropped into a calorimeter, 
there is usually a short interval of time during 
which the rate of temperature rise is much 
smaller than the rate which would be calculated 
by means of an equation of the form of (36). 
This is the result of various kinds of lag, espe- 
cially stirring lag. The form of experimental 
time-temperature curves after the temperature 
has begun to rise rapidly indicates also that the 
rate of transfer of heat between calorimeter and 
charge is not a linear function of (0—6@) as was 
assumed in the derivation of Eq. (36). It will 
be shown later that the value of the corrected 
temperature rise for some types of calorimetric 
experiments would differ from that calculated 
from Eq. (III) even if the time-temperature 
curve were of the form represented by Eq. (36). 

The calculation of the thermal leakage cor- 
rection for a bomb-calorimetric experiment will 
now be considered. The time-temperature rela- 
tion for such an experiment may be represented 
qualitatively by the curve ABDFG of Fig. 1. 
The parts AB and FG of this curve represent 
data obtained in the ‘‘rating periods’’ for use in 
calculating the thermal leakage. Even though 
the rates of change of © and @ may not be in 
accord with Eqs. (33) and (34) when (Q—@) is 
large, it will be assumed that these equations 
are valid during the rating periods when (0—6@) 
is small. 

During the rating periods Eqs. (36) and (I) 
reduce to 


0=0,+Ae-* (IV) 
and 

O=6.+Ce“. (V) 
Imposing the conditions (see Fig. 1) that 


6=6p3(=60) when t=¢g, and @=6r when t=t%,p, 
Eq. (IV) leads to the following equations for 
the initial and final rating curves 


@= 6.+(@2— 6, )es(tB—9 | 
6=0,.4+ (Or _ 0,.)estF—-, 


(VI) 
(VII) 


It is easily shown that the temperature of the 
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bomb will be represented by 


a 
© =6.+ -(O0n— 6, )et'*# ) (VILL) 
a—e 
during the initial rating period, and by 
a 
©=6.+ (Op — 6, Jet t?-? (IX) 
a-e 


during the final rating period. 

The observed temperature rise of the calo- 
rimeter will be taken as (@¢—6,). It follows 
from ‘Eqs. (VI), (VID, (VIID, and (LX) that if 
# and ©’ are the temperatures of calorimeter 
in the 
initial rating period, and 6” and ©” are the 
corresponding temperatures at any time ¢” in 
the final rating period, then 


and bomb, respectively, at any time (’ 


a 
Q/ =— 
G-~<¢ 


ee” — (@’’—6’). (X) 


In particular, the temperature rise of the bomb 
corresponding to the observed temperature rise, 
(0~e—6,), of the calorimeter will be given by 


a 
Op—Og= (O07 — Oz). 


a-e 


(XI) 


If the heat capacities of calorimeter alone and 
bomb alone are K, and Ko, respectively, the 
quantity of heat corresponding to the tempera- 
ture rise (@¢—6,) of the calorimeter, and the 
temperature rise (Or—Q,) of the bomb will be 
given by 


Ose =Ki(0r—62)+K2(Or—Oz) 
a 
= (Kit- K:) (0-0). (X11) 
a—e 


The heat gained by the calorimeter by thermal 


leakage between the times fg and ty will be 


given by 
tr 
aQ=—cKif (0—6,)dt. 
‘B 


This quantity may be taken as equivalent to an 
increase A@ in the temperature of the calorimeter, 
and an increase AQ = Aéa/(a—e) in the tempera- 
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ture of the bomb, where 


cK, f 
K,+Koea (a—e) tB 


The constant preceding the integral sign in 
(XIII) evidently may be considered as a cooling 
coefficient which could be used to calculate the 
increase in temperature A@ of the calorimeter 
due to thermal leakage on the basis that the 


A= — (0—6,)dt. (XIII) 

















corrected temperature rise, (@¢—@—A8@), of the 
_———— 
= 
- 
4 
c 
rr) 
a 
= 
is 
ro 
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Fic. 1. Diagram to illustrate Dickinson’s method of 


calculating the corrected temperature rise. 


calorimeter corresponds to a temperature rise of 
the bomb equal to 


a 
——(0~—60p,—A8). 


a-eé 


This cooling coefficient can be shown to be equal 
to the constant «. From Eqs. (IV) and (V) it 
follows that during the rating periods the rate 
of change of the temperature of the calorimeter 
is given by 


d0/dt = — «(@—0). (XIV) 


In other words the constant ¢ is the cooling 
coefficient for the calorimeter which would be 
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calculated from the data of the rating periods 
on the basis of Newton’s law of cooling as 
expressed by Eq. (XIV). Therefore the use of 
these data for calculating the correction to the 
observed temperature rise would lead to a 
corrected temperature rise, (@¢—6,—Aé6), for the 
calorimeter which would correspond to a temper- 
ature rise of the bomb equal to 


a 
—§(0r—0p—A6). 
a-e 


The fact that the corrected temperature rise 
of the calorimeter calculated from the rating 
period data is not equal to the corresponding 
temperature rise of the bomb is of no consequence 
in bomb-calorimetric measurements for the 
following reasons. In order to obtain accurate 
results by means of a bomb calorimeter the 
energy equivalent, or effective heat capacity of 
the calorimetric must be determined 
experimentally by supplying a measured quan- 
tity of energy to the system, either electrically 
or by burning a standard material of known 
heating value, and observing the _ resulting 
temperature rise. If the corrected temperature 
rise in the calibration experiments is calculated 
from the data of the rating periods on the basis 
of Eq. (XIV), then the value obtained for the 
energy equivalent will be that corresponding to 
a temperature rise of the bomb equal to a/(a—e) 
times that of the calorimeter. The experimental 
value of the energy equivalent is therefore not 
the true heat capacity of the calorimetric 
system. It is evident, however, that it is the value 
appropriate for use in calculating the results of 
measurements of heat of combustion when the 
thermal leakage correction for such measure- 
ments is computed as described above. 

The fact that the effect of lag of parts of a 
calorimeter is completely eliminated by an 
experimental calibration has been pointed out 
by White.* The preceding discussion in the 
present paper is a somewhat more detailed 
demonstration of this fact for the special case in 
which the lagging part of the calorimeter is a 
bomb. 


system, 


*W. P. White, The Modern Calorimeter (Reinhold 
Publishing Corporation, New York, 1928), p. 88. 
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It has been shown by Dickinson® that the 
correction for heat of stirring in a bomb-calori- 
metric experiment is included in the correction 
for thermal leakage, when this correction is 
calculated from the data of the rating periods in 
the manner described above. The effect of the 
production of heat at a constant rate by stirring 
is merely to raise the convergence temperature 
slightly above the temperature at which the rate 
of heat transfer between calorimeter and sur- 
roundings would be zero. The heat of stirring 
is also taken care of in the calculation of (0,,.* — 60) 
by means of Eq. (III) if the value of @, used in 
this equation is the convergence temperature of 
the stirred calorimeter. The temperature 6,,* is 
then the temperature which the system would 
attain if the rate of loss of heat were at all times 
equal to the rate of production of heat by stirring. 

As an example of the correction of bomb- 
calorimetric data for thermal leakage, Dickinson's 
method? of calculating the corrected temperature 
rise will now be applied to a_ hypothetical 
bomb-calorimetric experiment for which the 
time-temperature relation is that represented by 
Eq. (II). The time-temperature relation during 
the initial rating period for such an experiment 
would be given by 


6=1.7082(1 —e—9-95*), (XV) 
This equation was obtained by substituting the 
numerical values of 6z, 6., and ¢« given by King 
and Grover into Eq. (VI). The corresponding 
relation for the final rating period is represented 
by Eq. (II) with the last term omitted. 
Dickinson’s method of finding the corrected 
temperature rise consists in extrapolating the 
rating curves AB and FG (Fig. 1) to a time f, 
such that the areas BCD and DEF are equal. 
The difference between the temperatures 6¢ and 
§¢ obtained in this way is the corrected tempera- 
ture rise. The value of t, for the experiment under 
consideration was calculated from the relation 


f Atarevart f Be-2r-9'dt=0 (XVI) 


and found to be 2.01 minutes. The values of 6¢ 
and 6, calculated by inserting {=/,=2.01 into 


5H. C. Dickinson, Bull. Bur. Stand. 11, 189 (1914). 
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the equations for the rating curves are, respec- 
tively, 0.01708° and 2.70813°. The difference 
between these temperatures is 2.69105°, which 
is lower by 0.073 percent than the value 2.69302° 
calculated from Eq. (IIT) for (6.* — 0). 

The reason for the discrepancy is that the two 
values of temperature rise of the calorimeter 
correspond to different values for the tempera- 
ture rise of the bomb. Thus the temperature rise 
(6.*—6) for the calorimeter corresponds to a 
temperature rise of (@,*—@Q,) for the bomb, 
where © ,, in accordance with Eq. (VIII), is 
given by 

Op = (a0z — €6;.) ‘(a —e). 


Hence the temperature rise for the bomb in 
this case is 
ab p— €). 


0,* —On=0,* — (XVII) 


a—e 
On the other hand the temperature difference 
(@e—86-) for the calorimeter, in accordance with 
Eq. (X), corresponds to a temperature difference 
of 
a 
Or — Oc =——(0r— 6c) 
a—e 


(XVIII) 


for the bomb. 

It is easily shown by actually performing the 
calculation that the quantity of heat required 
to produce the temperature rise (6,*— 4) 
= 2.69302° in the calorimeter, and the corre- 
sponding temperature rise (6,.*— © ,) =2.71376° 
in the bomb is exactly the same as the quantity 
of heat required to produce the temperature 
increases (@¢—8c¢)=2.69105° and (Og—Qc¢) 
= 2.72373° in calorimeter and bomb, respectively. 
Hence either (6,*—6)) or (@g—@c) may be 
considered as the corrected temperature rise for 
the bomb-calorimetric experiment under con- 
sideration if used with the proper value for the 
energy equivalent of the system. From Eq. 
(XVIII) and the discussion given in the text 
following Eq. (XIV) it will be seen that the 
value of the energy equivalent appropriate for 
use with the temperature rise (@¢—4@c) is the 
value which would be obtained as the result of 
an experimental calibration of the. calorimeter 
in which the thermal leakage correction was 
calculated on the basis of Newton’s law of 
cooling in the form (XIV). 
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TABLE I. Values of temperature and corresponding time in 
a hypothetical bomb-calorimetric experiment. 








Time Temperature Differences | Time Temperature  Differentes 
minutes degrees degrees minutes degrees degrees 
—10 —0.08758 12 2.6527 
— 9 — 07863 900895 | 14 2.64747 
= .00892 — 0.00372 
— 8 — .06971 15 2.64375 
. < .00886 — .00409 
— 7 — 06085 16 2.63966 
.00883 pn i. — .00429 
— 6 — .05202 a 17 2.63537 
= .00878 — .00441 
— 5 — .04324 18 2.63096 
.00873 — .00447 
— 4 — 03451 19 2.62649 
.00869 — .00450 
— 3 — .02582 20 2.62199 
. .00864 — .00451 
— 2 — 01718 21 2.61748 
.00862 vs — .00451 
— 1 — .00856 00856 22 2.61297 — '00449 
0 .00000 : P 23 2.60848 ra ‘00448 
0.5 5987 24 2.60400 = 00446 
1 1.0645 25 2.59954 ee ‘00445 
2 1.7082 26 =2.59509 : 
ns pa — .00442 
+ 2.3304 27 2.59067 — ‘00439 
6 2.5591 28 2.58628 a ‘00438 
8 2.6288 29 2.58190 pe ‘00436 
10 2.6507 30 =.2.57754 . . 


The value of the energy equivalent of a bomb 
calorimeter appropriate for use with the temper- 
ature rise (6*,—6 ) depends to a slight extent 
upon the magnitude of this temperature rise, 
apart from any dependence of heat capacity 
upon temperature. In general it cannot be as 
conveniently determined experimentally as the 
energy equivalent corresponding to the tempera- 
ture rise (@g—6c). In the case of an experiment 
by the method of mixtures, the value of the 
energy equivalent which is equal to the effective 
heat capacity of the calorimeter alone, multiplied 
by (@.*—6) will yield the heat given up by the 
charge in cooling from its initial temperature to 
6,*. This energy equivalent can be determined 
experimentally either by electrical heating of the 
calorimeter alone, or by the use of a charge of 
known heat capacity. However, in general it 
will not be possible to use Eq. (II1) to evaluate 
(6,*—6 ) either for bomb-calorimetric experi- 
ments or for experiments by the method of 
mixtures, because of the fact that experimental 
time-temperature relations are not of the form 
assumed in the derivation of Eq. (III). The 
calculation of the thermal leakage correction 
for an experiment by the method of mixtures 
will be considered later in this paper. 

In practice the time-temperature relation for 
a calorimetric experiment is not given in terms 
of equations, but in terms of temperatures 
measured at definite times. Ordinarily, therefore, 
the method used above of calculating the 
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“ABLE II. Calculation of corrected temperature rise for a 
hypothetical bomb-calorimetric experiment. 


Initial period Final period 


64 = —0.08758° OF 


= 2.62199° 
6B = 0.00000 @6G= 2.57754 
Aé@= 0.08758 . Aé@ = —0.04445 
Aé atl c A@ a p 
) = 0.008758 deg. min.~! ) = —0.004445 deg. min.~! 
M/ p Mi/s 
ép = —0.04324° @s= 2.59954° 
Aé . 
4( 2 ) = —0.013203 deg. min.~! 
A@= 2.64278° 
Aé . 
A( =) Aé = —0.0049959 min.~! 
2.01 ce 17.99 ‘ 
iq-t3=-——— = 1.005 min. trp —tr= a = 8.995 min. 
1.005 X0.008758 = 0.00880° —8.995 K —0.004445 = 0.03998° 
6p= 0.00000 Or = 2.62199 
690= 0.00880 6re= 2.66197 
6p = —0.04324 @s= 2.59954 
A@= 0.05204 Aé = —0.06243 
0.05204 KX —0.0049959 +0.06243 K —0.0049959 
= —0.000260) ‘deg. = —0.000312 deg. min.~! 
A@ : Aé 
= 0.008758 min.~! = —0.004445 
M P ‘ Al Ss 
A@ : Aa 
[v ) 0.008498 deg. min.~! ) = —().004757 deg. min.~! 


At /R 
—0.004757 X —17.99 =0.08558° 
6p =2.62199 


0.008498 X 2.01 =0.01708° 
@B =0.00000 


6c =0.01708' 60 =2.70757° 


Corrected temperature rise = 2.70757 —0.01708 =2.69049° 


temperatures @¢ and 6g from equations for the 
rating curves would not be applied. It will be of 
interest to consider the calculation of the 
corrected temperature rise for a bomb-calori- 
metric experiment, using numerical data rather 
than equations. In Table I are given values of 
temperature calculated by means of Eqs. (II) 
and (XV). It will be assumed that these temper- 
atures were obtained from observations in a 
bomb-calorimetric experiment, and they will be 
used to calculate the corrected temperature rise 
of the calorimeter by Dickinson's method. The 
period from t= —10 to t=0 will be taken as the 
initial rating period. Since the equation for the 
final rating period is assumed not to be known, 
it would be necessary for the experimenter to 
judge from the data when the final rate of 
temperature change has been attained. If 
temperatures were read each minute to 0.0001° 
it would be evident to the experimenter that the 
final rate had not been attained before t=18, 
and he might be expected to take the period 
from t= 20 to t=30 as the final rating period. 
Ordinarily the value of t, would be obtained 
from numerical data by plotting such data on 
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large-size coordinate paper and determining 
areas by some such method as counting squares. 
In the present case, however, we are interested 
in comparing the results of extrapolating the 
rating curves to t, by means of the equations for 
these curves, with the results of extrapolating 
the curves by the use of numerical data, without 
introducing the effect of possible errors in the 
determination of ¢,. Therefore the value t,=2.01 
which was calculated by means of Eq. (XVI) 
will be used here. 

Since the present writer, using Dickinson's 
method, has obtained a distinctly different 
value for the corrected temperature rise from 
that reported by King and Grover as having 
been obtained by the same method, the details 
of the present calculation are given in Table IT. 
The . following description of the method of 
calculation will be understood by reference to 
Table II and Fig. 1. 

The average rate of change of temperature, 
A6/At, in each of the rating periods is taken as 
corresponding to the ‘‘observed’’ temperature 
at the mid-time of that period. From the values 
of the two rates and the corresponding tempera- 
tures the average change in rate per degree 
change in temperature, A(A@/At)/A@(= —e), is 
calculated. The approximate temperature, 69, 
which would have been attained at the time tg 
(Fig. 1) midway between fz and ¢, if the’ rating 
period had continued to the time fe, is then 
calculated by linear extrapolation of the curve 
AB. From the value of A@/At)p, the difference 
between the temperatures @g and 6p (Fig. 1), 
and the value of A(A@/At)/Aé@, the value of the 
rate (A@/At)g corresponding to the temperature 
6 is calculated, and this rate is multiplied by 
(t,—tg) to obtain (@¢—6g). Similarly, the rate 
corresponding to the temperature 6” (Fig. 1) is 
calculated and multiplied by (tp—t,) to obtain 
(0z—O6,). These data, together with the values 
of 6g and 6, give the corrected temperature rise 
(6z—O6¢). 

The value of the corrected temperature rise 
obtained in the manner just described is 2.69049°, 
which is lower by 0.021 percent than the corre- 
sponding value 2.69105° calculated by means of 
the equations for the rating curves. Three 
factors contribute to this discrepancy: (1) 
Although the experimenter would probably think 
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that the final rate had been attained at time 
t=20, the calculated value of the last term of 


Eq. (II) for t=20 is 0.00012°, so that the 
observed value of the temperature for t= 20 is 
lower by this amount than the temperature 
corresponding to the final rating curve. (2) The 
linear extrapolation of the final rating curve 
gives too low a value for @z. (3) Although dé/dt 
in the rating periods varies linearly with temper- 
ature, A@/At does not vary linearly with the 
temperature at the midtimes of the intervals for 
which the values of A@/At are calculated.® All 
three of these factors enter into the calculations 
in such a manner as to lower the resulting value 
of @¢. The corresponding error in @¢ arising from 
factors (2) and (3) is negligible because the 
extrapolation of the initial rating curve is over 
such a short interval. 

The errors arising from the factors (2) and (3) 
mentioned above result from the fact that the 
method of extrapolating the final rating curve by 
the use of the numerical data does not take full 
account of the curvature of this curve. The 
small error from this source could be avoided by 
so conducting the experiment that the final 
temperature of the calorimeter is very nearly 
equal to the convergence temperature. The final 
rate would then be very small, and the change 
in rate between the temperatures 6g and @s 
(Fig. 1) would be negligible. 

The hypothetical experiment under considera- 
tion differs in several important respects from 
an actual bomb-calorimetric experiment with 
apparatus suitable for measurements of high 
precision. The differences are such as to make 
the errors in the calculation of the thermal 
leakage correction from rating period data larger 
for the hypothetical experiment than for an 
actual experiment. The curves of Fig. 2 represent 
the time-temperature relation for the hypo- 
thetical experiment, and a_ time-temperature 
relation obtained by multiplying the observed 
temperatures in a routine bomb-calorimetric 
experiment by the ratio of the corrected tem- 
perature rise, 2.69105°, for the hypothetical 
experiment to the corrected temperature rise, 
2.9519°, calculated from the data of the actual 


* Except in the special case that the intervals are all 


equal. The three intervals involved in the calculation of 6 


are tg—ta4=tg—tr=10, and tre—t,=17.99. 


134 





experiment. The “observed” curve of Fig. 2 
may therefore be considered as representing 
the data of an actual bomb-calorimetric ex- 
periment for which the corrected temperature 
rise is the same as for the hypothetical experi- 
ment. The final rate in the actual experiment was 
about 0.0001° per minute, so that the change in 
rate corresponding to the temperature difference 
(@x—O@s) (Fig. 1) is negligible. 

It will be seen that the observed curve (Fig. 2) 
approaches a steady final rate much 
rapidly than the hypothetical curve. Even 
though the observed curve rises more slowly at 
the start, its later rise is so rapid that the final 
rate is reached at t=10, while for the hypo- 
thetical curve the final rate is not reached until 
about ‘=20. This more rapid approach of the 
actual curve to the final rate, and the consequent 
smaller interval over which the observed final 
rating curve must be extrapolated to obtain 6z, 
result in a smaller error in this temperature 
arising from the factor (1) mentioned above, 
even if the error in 6 is the same as in the 
hypothetical experiment. The value of ¢ for the 
actual experiment is less than half that assumed 
for the hypothetical experiment. The thermal 
leakage correction for the actual experiment is 
less than for the hypothetical experiment because 
of the smaller value of e, because of the shorter 
time of the actual experiment, and because 

(@—6@.) is smaller during a relatively large 
part of the actual experiment. This correction 
amounts to less that 0.3 percent for the actual 
experiment, while for the hypothetical experi- 
ment it is over 2 percent. 

Although the above discussion has been 
restricted to Dickinson’s method of calculating 
the thermal leakage correction, it is evident that 
the proper use of any method based on Newton’s 
law of cooling can lead to only one result. 

The significance of the constants c and € occur- 
ring in Eqs. (34) and (36), respectively, and the 
relation between them will now be discussed 
briefly. It is evident from Eq. (34) that the con- 
stant c is the cooling coefficient of the calorimeter 
alone, that is, not including a bomb or other 
charge. From (VI) and (VII) it is evident that 
the cooling coefficient of a calorimetric system 
which includes a bomb or other charge is not ¢ 
but e. From the definition of ¢ it follows that if 
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(e’— ec) is negligible in comparison with ac, then 
e=ca/(a+b). In other words, the ratio ¢/c is 
approximately equal to the ratio of the heat 
capacity of the calorimeter alone to the heat 
capacity of calorimeter plus bomb. Thus, from 
the values of a, 6, and c given by King and 
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Fic. 2. Comparison of the hypothetical time-temperature 
curve represented by Eq. (II) with a time-temperature 
curve observed in a routine bomb-calorimetric experiment. 


Grover the value of ca/(a+) is calculated to 
be 0.00501 as compared with their value 0.00500 
for e. The fact that € is not exactly equal to 
ca/(a+b) is due to the lag of the temperature 
of the charge or bomb behind that of the 
calorimeter. 

The discussion of the thermal leakage cor- 
rection for the bomb-calorimetric 
experiment may be summarized as follows. The 
data of the rating periods, which correspond to 
the cooling coefficient ¢ rather than to c, are the 
appropriate data to calculating the 
thermal leakage correction on the basis of 
Newton’s law of cooling. The results of such 
calculations will be affected by the lag of the 
temperature of the bomb behind that of the 


case of a 


use in 
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calorimeter, but errors due to this lag will be 
completely eliminated by the use of an experi- 
mentally determined value of the energy equiva- 
lent of the calorimetric system. Errors in the 
calculation of the thermal leakage correction can 
be made negligible by suitable experimental 
procedure, especially by so conducting the 
experiments that the final temperature of the 
calorimeter will be approximately equal to the 
convergence temperature of the system. 

The above conclusions were reached on the 
basis of certain assumptions regarding the 
conditions existing during a bomb-calorimetric 
experiment, which were introduced in order to 
simplify the mathematical treatment. The con- 
clusions have been shown to be valid, therefore, 
only for the type of ideal experiment defined by 
the assumptions. King and Grover, in their 
discussion, have introduced essentially the same 
simplifying assumptions, and therefore their 
conclusions are necessarily subject to the same 
restrictions. The treatment given in this paper 
is not a complete justification of the usual 
methods of correcting for thermal leakage, but 
it does show the fallacy of the conclusion of 
King and Grover that the usual methods are 
not applicable to the type of ideal bomb- 
calorimetric experiment which has been discussed. 

A rigorous mathematical treatment, of the 
problem of correcting the data of actual bomb- 
calorimetric experiments for thermal leakage 
will not be attempted here. However, it appears 
from the following discussion that the errors 
involved in the usual methods of correcting such 
data are probably negligible. It follows from 
the mathematical analysis given in this paper 
that, in the case of the ideal type of experiment 
considered, the effect of the difference between 
the temperatures of and_ calorimeter 
during the rating periods is fully taken into 
account by an experimental calibration of the 
calorimeter. It is evident that this would be 
true also for actual calorimetric experiments 
provided the time-temperature curves of calibra- 
tion experiments and combustion experiments 
are identical. This condition is practically satis- 
fied when the calibration is by means of a 
standard substance of known heat of combustion, 
but is not satisfied when the calibration is by 
electrical heating. It is conceivable, therefore, 


bomb 
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that when a bomb calorimeter is calibrated by 
electrical heating the difference in temperature 
between bomb and calorimeter during the rating 
periods might not be taken into account com- 
pletely, although it is evident that this difference 
would be at least approximately taken into 


account. For the calorimeter used in obtaining 
the data represented by the observed time- 
temperature curve of Fig. 2 it has been calcu- 
lated that the difference between the tempera- 
tures of bomb and calorimeter affects the 
observed temperature rise by about 0.03 percent. 
It seems highly probable that the error resulting 
from the difference in form of the time-tempera- 
ture curves of electrical calibration experiments 
and combustion experiments is a small fraction 
of this amount. 

This supported by results 
obtained by Dickinson’ in two sets of electrical 


conclusion is 


calibration experiments one set with a heater 
which was not in contact with the bomb and the 
other with a heater which was in good thermal 
contact with the bomb. The time-temperature 
curves for the first set of experiments were nearly 
linear during the heating period, while those for 
the second set were somewhat similar in form 
to the time-temperature curve for a combustion 
experiment. The results of the two sets of 
experiments were in practically perfect agree- 
ment. This indicates that the difference in form 
of the time-temperature curves had no measur- 
able effect on the corrected temperature rise 
calculated by the usual methods. 

The method described above for calculating 
the thermal 
calorimetric experiment is not applicable to an 
experiment by the method of mixtures, because 
in this case the data of the initial and final rating 
periods correspond to different heat capacities 
of the calorimetric system, and because of the 
fact that the effect of the lag of the charge in 
such an experiment is not eliminated by an 
experimental calibration of the calorimeter. A 
method of calculation appropriate for experi- 


leakage correction for a bomb 


ments by the method of mixtures will now be 
described. 

It is assumed that the charge is dropped into 
the calorimeter at time tg (when @=6,), and 
that the mass and initial temperature of the 
charge are such that the final temperature of 
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the calorimeter will be so nearly equal to 6, that 
the lag of the charge will not cause its tempera- 
ture to differ by a significant amount from that 
of the calorimeter during the final rating period. 
The final rating period is assumed to begin at a 
time t=fp. 

The cooling coefficient of the calorimeter alone 
will be represented by c. The rate of increase in 
temperature of the calorimeter alone due to 
thermal leakage will then be represented by 


dé /dt= —c(@—9;). (XTX) 


If K, represents the heat capacity of the calo- 
rimeter alone, the expression 


dQ dt= —cK,(@—94,) (XX) 
will represent the rate at which heat is added to 
the calorimeter by thermal leakage. This rate is 
assumed not to depend on whether the charge is 
or is not in the calorimeter. 

From Eq. (XX) it follows that the heat 
added to the calorimeter by thermal leakage 
between the times fg and fp is given by 

tr 
AQ= = Ki f (0—6,.)dt. (XX1) 
tB 
That part of the temperature rise (6”— 6s) which 
is due to thermal leakage is therefore given by 


AQ tr 
A§=——— --af (6—6,)dt, 
K,+K,. tB 


where Kez represents the heat capacity of the 
charge and 


(XXIT) 


K, 
a = c— 
Ki+Ksz 
is the coefficient of 
calorimeter plus charge. 

The constant a may be calculated from the 
data of the rating periods. Thus if (A@/A?) p» and 
(A@/At)s are the average rates of change of 
temperature during the initial and final rating 
periods, respectively, we may write 


Aé . 
( -) = —c(Op:— ,) 
At7 p: 


Ad Ki 
(—) ike ltl, 
At/ s Ki+ Kz. 


(approximate) cooling 


JOURNAL OF APPLIED PHYSICS 











whence 
K, 
( — 
K,+ Ky, 


Ky, (~) (=) 
K,+Ke\ At/ p’ At? s. (XXIII) 


Op —Os 


The constant —a is thus equal to the value of 
A(A@/At)/A@ for a calorimeter having a_ heat 
capacity of (K,+Ke), and thermal 
represented by Eq. (XIX). 

The thermal leakage correction to the observed 
temperature rise can be calculated by inserting 
the value of a into Eq. (XXII), and evaluating 
the integral by any convenient method. As the 
object of an experiment by the method of mix- 
tures is usually the determination of Ke, it will 
be necessary to use an approximate value of this 
constant in calculating a, and more than one 
approximation may be required to obtain the 
desired accuracy. An alternative, and perhaps 
better procedure, would be to calculate the 
quantity of heat represented by Q=K,(@r—96,) 


leakage 


and to calculate the thermal leakage cor- 
rection to this quantity by means of Eq. 
(XXI), using the value of c calculated by 


means of Eq. (XXIII). An error in the value 
used for K2 would then be multiplied by (A@/Az) s, 
which is small, rather than by (A@/Af) p-, wich is 
large, and the error would therefore have a 
much smaller effect on the final result. 

It should be emphasized that the above 
procedures for correcting the data of an experi- 
ment by the method of mixtures are based on 
the assumption that the final rate is so small 
that the lag of the charge does not affect the 
final temperature by a significant amount. The 
effect of the lag in such experiments is not 
eliminated by an experimental calibration of 
the These procedures are 
subject to the restriction that the rate of pro- 
duction of heat by stirring does not change when 
the charge is dropped into the calorimeter. 

It appears from the discussion in this paper 
that the important contribution of King and 


calorimeter. also 
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Grover consists in pointing out that the lag of 
the bomb in a determination of heat of com- 
bustion, or of a solid body introduced into a 
calorimeter in an experiment by the method of 
mixtures, affects the observed temperature rise 
of the calorimeter, and the corrected temperature 
rise calculated by the usual methods. It is 
possible that failure to appreciate this fact has 
been a source of error in bomb-calorimetric 
measurements in cases where the energy equiva- 
lent of the system was calculated from the 
masses and specific heats of the various compo- 
nent -parts of the calorimeter, and in measure- 
ments by the method of mixtures made under 
conditions such that the rate of temperature 
change in the final period was large. 

The effect of the lag of the bomb in a bomb- 
calorimetric experiment is practically completely 
eliminated by an experimental calibration of the 
calorimeter. In this case, any method of calcu- 
lation based on Newton’s law of cooling, if 
properly applied, should yield the correct result. 
Small errors in the calculation of the thermal 
leakage correction can be reduced by arranging 
the experiment so that the final temperature is 
very near the convergence temperature, and the 
final rate, therefore, very small. If the effects of 
evaporation are to be considered, it is preferable 
that this final temperature be below the con- 
vergence temperature rather than above it, so 
that any heat transfer by evaporation will be 
always in the same direction. Since experimental 
calibration of the calorimeter has been the usual 
practice in precise measurements for twenty-five 
years or more, the results of modern bomb- 
calorimetric experiments are not affected’ by 
the errors discussed by King and Grover. 

In the method of mixtures, in which a solid 
body is introduced into the calorimeter at the 
end of the initial rating period, the effect of the 
lag of the body can be practically eliminated by 
so arranging the experiment that the final 
temperature is at or very slightly below the 
convergence temperature. In this case also the 
usual methods of calculation based on Newton’s 
law of cooling should yield the correct result. 


7See F. D. Rossini, Chem. Rev. 18, 243-244 (1936). 
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Placement Register for Physicists 


Hike American Institute of Physics has main- 

tained a relatively rudimentary form of 
placement service since 1937 which has func- 
tioned simply as an agency for bringing qualified 
applicants desirous of finding positions in touch 
with employers seeking men. It has, during its 
existence, made approximately 175 such con- 
tacts, some of them calling for the employment 
of several men. 

Recently it has become apparent that, due to 
rapid changes of the status of applicants, to a 
shortage of applicants, and to the ever increasing 
‘ demands made on it, this system was not per- 
forming its function as effectively as it might. 
Early in 1940 it was deemed desirable to con- 
template changes in it and in April of that year 
a study was made, with the cooperation of the 
American Chemical Society, of the Placement 
Service used by that organization. 

As a result, the placement service is discon- 
tinuing its present system of operation and plans 
to put into operation a system somewhat similar 
to. that used by the American Chemical Society. 
Under this plan, tentatively called a Placement 
Register, a listing of physicists interested in new 
employment will be made available to prospec- 
tive employers at certain important meetings of 


the Founder Societies. Data on applicants will 
be less complete than that asked for under the 
old system and particular attention will be given 
to applicants in attendance at the meeting and 
who are available for interview by the employer 
or his representative. A separate listing of appli- 
cants not at the meeting will also be available to 
the employers. If it seems desirable, the register 
files will be retained for a relatively short period 
after the meeting for which they were made up 
and names from these files may be referred to in- 
terested employers unable to attend the meeting. 

Applicants will register at the meeting as near 
its opening as possible or may pre-register prior 
to this time on special forms to be supplied by 
the Institute. 

It is expected that this plan will be tried out 
at the spring meeting of the American Physical 
Society at Baltimore at the end of April. It 
should be noted that the venture is in the nature 
of an experiment which will succeed only with 
the full cooperation of all physicists and em- 
ployers of physicists. 

Further regarding methods of 
registration and details of operation will appear 
in succeeding issues of this journal, this being in 
the nature of a preliminary notice. 
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New Instrument Booklets 








The latest issue of The Laboratory, publication of the 
Fisher Scientific Company and Eimer & Amend, features an 
article entitled ‘‘The Pharmaceutical Chemist,’’ and con- 
tains numerous other articles describing the uses of Aerosol 
and such new apparatus as the Snelling-Fisher Hotspotter, 
a device for providing concentrated variable heat; the 
Fisher-Serfass Electronic Relay; the Fisher-Engel Heat of 
Vaporization apparatus; and a d.c.-a.c. inverter to convert 
direct current to alternating. Further information may be 
obtained from Eimer & Amend, 633-635 Greenwich Street, 
New York, New York, or from the Fisher Scientific Com- 
pany, 711-723 Forbes Street, Pittsburgh, Pennsylvania. 


The General Radio Experimenter, October, 1941, issue, 
contains articles on, “‘A vacuum-tube-driven tuning-fork 
oscillator” and ‘Impedance bridges assembled from labora- 
tory parts—Part IV.”’ The November issue contains 
“Impedance bridges assembled from laboratory parts— 
Part V,” as well as an article entitled, “‘A redesign of the 
vacuum tube bridge.”’ Part VI of the series on ‘‘Inductance 
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Measurements” appears in the December issue of The 
General Radio Experimenter, with a second article entitled 
“Using the Cathode-Ray Oscillograph in Frequency Com- 
parisons.”’ Published by the General Radio Company, 
30 State Street, Cambridge A, Massachusetts. 


The Ohmite News, December, 1941, describes the use of 
Ohmite rheostats in the new WABC transmitter at 
Columbia Island in Long Island Sound, and gives additional 
information about other resistors and rheostats made by 
the Ohmite Manufacturing Company, 4835-41 Flournoy 
Street, Chicago, Illinois. 


Harrison and 
Peoria Streets, Chicago, Illinois, has just released a new 
condensed which Bulletin 
No. Z4000 


The Wheelco Instruments Company, 


catalog supersedes the old 


The latest issue of Mechanical Topics, a publication of 
the International Nickel Company, Inc., 67 Wall Street, 
New York, New York, presents new developments in the 
use of nickel alloys of interest to metal users who are con- 


’ cerned with war needs today. 
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